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ESTIMATION OF BACTERIAL DENSITIES BY MEANS OF 
THE “MOST PROBABLE NUMBER’* 


WiiuramM G. CocHRAN 


School of Hygiene and Public Health 
The Johns Hopkins University 


Presented at a joint session of the Engineering, Laboratory and Statistics 
Sections of the American Public Health Association, Biometrics Section of the 
American Statistical Associations and the Biometric Society, New Y. ork, 
October 27, 1949. 


INTRODUCTION 


HIS PAPER attempts to give a simple account of the concept of the 

“most probable number” (m.p.n.) of organisms in the dilution 
method. The concept is quite old, going back to MeCrady (4) in 1915, 
and has been discussed by various writers from time to time, so that 
little of what I shall present is new. In addition, some advice is given on 
the planning of dilution series. 

The dilution method is a means for estimating, without any direct 
count, the density of organisms in a liquid. It is used principally for 
obtaining bacterial densities in water and milk. The method consists 
in taking samples from the liquid, incubating each sample in a suitable 
culture medium, and observing whether any growth of the organism 
has taken place. The estimation of density is based on an ingenious 
application of the theory of probability to certain assumptions. Tor a 
biologist, it is more important to be clear about these assumptions than 
about the details of the mathematics, which are rather intricate. 


ASSUMPTIONS 


There are two principal assumptions. In statistical language, the 
first is that the organisms are distributed randomly throughout the liquid. 
This means that an organism is equally likely to be found in any part of 
the liquid, and that there is no tendency for pairs or groups of organisms 
either to cluster together or to repel one another. In practice this implies 
that the liquid is thoroughly mixed, and if the volume of liquid is not too 
great some shaking device is usually employed for this purpose. 


*Paper 254 from the Department of Biostatistics, 
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The second assumption is that each sample from the liquid, when 
incubated in the culture medium, is certain to exhibit growth whenever 
the sample contains one or more organisms. If the culture medium is 
poor, or if there are factors which inhibit growth, or if the presence of 
more than one organism is necessary to initiate growth, the m.p.n. gives 
an underestimate of the true density. 


MATHEMATICAL ANALYSIS 


In the mathematical analysis we relate the probability that there will 
be no growth in a sample to the density of organisms in the original hquid. 
Suppose that the liquid contains V ml., the sample contains v ml., and 
that there are actually b organisms in the liquid. By the second assump- 
tion, there will be no growth if and only if the sample contains no organ- 
isms. We will calculate the probability that none of these b organisms is 
in the sample. 

Consider a single organism. By the first assumption, the probability 
that it lies in the sample is simply the ratio of the volume of the sample 
to that of the liquid, i.e.v/V. The probability that it is not in the sample 
is therefore (1 — v/V). Since there is assumed to be no kind of attract- 
tion or repulsion between organisms, these two probabilities hold for 
any organism, irrespective of the positions of the other organisms. 
(Strictly, this requires the additional assumption that the space occupied 
by an organism is negligible relative to v.) Consequently, by the multi- 
plication theorem in probability, the probability that none of the b 
organisms is in the sample is 


p= (1 =»v/V). 


When v/V is small, this is closely approximated by 


—vb/V 


De 


where e, about 2.7, is the base of natural logarithms. Finally, since 
b/V is the density 6 of organisms per ml., we have 


| = eae 
where p is the probability that the sample is sterile. 


THE CASE OF A SINGLE DILUTION 


If n samples, each of volume v, are taken, and if s of these are found 
to be sterile, the proportion s/n of sterile samples is an estimate of p. 
Hence we obtain an estimate d of the density 6 by the equation 
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This gives 


i . In () == 2205 log () (1) 


nr v 


where /n and log stand for logarithms to base e and to base 10 re- 
spectively. 

The estimate d is the ‘‘most probable number”’ of organisms per ml. 
The derivation given here does not reveal why this name has been 
ascribed to the estimate. In fact, the concept of m.p.n. is scarcely needed 
for this simple case. We will, however, reexamine the analysis so as to 
introduce the concept, which becomes useful in the more complex situa- 
tion where several dilutions are used. 

If p is the probability that a sample is sterile, the probability that s 
out of n samples are sterile is given by the binomial distribution as 


ni 
ely. —~s)! 


0 2) ae (2) 
Since p = e”’, this expression may be written 


mt , 
CE s)! ie ans 5 a (3) 
If we have obtained s sterile samples out of n, this formula enables us to 
plot the probability of this event against the true density 6. Such curves 
always have a single maximum. 

A curve of this type suggests a method for estimating 6 which is 
plausible on intuitive grounds. For if we are considering two possible 
values of 6, it seems reasonable to prefer the one which gives a higher 
probability to the result that was actually observed. This argument, 
carried to its conclusion, leads to a choice of the value of 6 for which 
the probability of obtaining the observed result is greatest. It is this 
value of 6 that has been called the ‘“‘most probable number” of organisms. 
It can be shown mathematically that this is the value of 6 for which 
p = s/n. Consequently the m.p.n. is the same as the estimate previously 
given. 

In practice, more than one dilution is usually needed. The reason is 
that the precision of the m.p.n. is very poor when the volume v in the 
sample is such that the samples are likely to be all fertile or all sterile. 
When all are fertile, the maximum on the probability curve (3) occurs 
when 6 is infinite, so that the estimated density is infinite. When all are 
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sterile the estimated density is zero, as may also be verified from equation 
(1). Thus a single dilution is successful only if v happens to be chosen 
so that some samples are sterile and some are fertile. Such a choice of 
v can be made only if the density 6 is known fairly closely in advance. 

If we possess this knowledge, it is best to select v so that the expected 
number of organisms per sample lies somewhere between 1 and 2. For 
this choice the expected percentage of sterile samples will lie between 
15% and 35%. In default of this knowledge, the practice is to use 
several dilutions (i.e. several different values of v) in the hope that at 
least one of them will give some sterile and some fertile samples. 


THREE DILUTIONS 


The case of three dilutions serves to illustrate the general problem. 
Let the suffix 7 indicate the dilution. For the 7th dilution the volume of 
the sample is v; , and s; out of n,; samples are found to be sterile. How 
do we estimate 6 from these results? 

From equation (1) we can obtain a separate estimate for each dilu- 
tion: 1.e. 


i ee log e) 


te 


However, the best way to combine the three estimates d; into a single 
value is not obvious. Since, as we have seen, some dilutions give very 
poor estimates, it is not satisfactory to take the arithmetic mean. 

One solution is provided by the m.p.n. concept, which extends easily 
to this situation. Following the approach used in the previous section, 
we first write down the probability of obtaining the observed results for 
any hypothetical value of the true density 6. The observed results are 
that s,; samples out of n, are sterile at the first dilution, s, out of m2 at 
the second, and s; out of nm; at the third. The probability that these 
three events should all happen is the product of three terms, each like 
expression (3) in the previous section. As before, the graph of this 
probability against 6 shows a single maximum. The value of 6 at this 
maximum is taken as the m.p.n. 

The value of the m.p.n. cannot be written down explicitly. The 
equation which it satisfies is as follows. 


_ —vid —_ era <= = 
SV =5 Se + $303 = (m = sie ™ + (M2 = 82)v,¢"*" ae (ns 83)V3€ 


—vid —va 
== 1—e 


Methods for solving this equation by trial and error have been given by 
several writers: e.g. Halvorson and Ziegler (3), Barkworth and Irwin (1) 
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and Finney (2). In laboratories where the numbers of samples n; and 
the dilution ratios are standardized, it is convenient to have a table which 
gives the m.p.n. for all sets of results that are likely to occur. <A table 
is provided in “Standard methods for the examination of water and 
sewage” (5), for dilution series in which 5 samples are taken at each 
dilution and there are three 10-fold dilutions. A more extensive table, 
for dilution ratios of 2, 4, and 10 and any number of levels (except two 
levels with a 10-fold dilution) is given by Fisher and Yates (6). This is 
not a table of the m.p.n., but of a different estimate which seems to be 
just about as precise for series of the size usually conducted in practice. 
This estimate is derived from the total numbers X and Y of fertile and 
sterile samples. The quantities x = X/n, y = Y/n are entered in the 
table, from which an estimate of log d is obtained. 


CRITIQUE OF THE M.P.N. 


We have seen that the m.p.n. is an estimate of the density of organ- 
isms. Considered more generally, it is a procedure for obtaining estimates, 
since the same argument could be applied to other statistical problems. 
The only justification which I have mentioned for the procedure is 
that it seems intuitively reasonable. From a reading of the literature 
I am not certain as to the reasons which led early investigators to select 
this estimate, though either the intuitive approach or an appeal to a 
theory of inverse probability may have been responsible. 

During the past 25 years the problem of making estimates from data 
has received much attention from statisticians. Today, most statis- 
ticians would, I believe, reject an appeal to intuition or to the theory of 
inverse Probability as a reliable procedure for constructing estimates, 
since both have been found on occasion to be untrustworthy. They 
might also object to the name “most probable number,” on the grounds 
that the adjective ‘probable’ in that phrase has a different meaning 
from the one given to it in the theory of probability. The estimate is 
“most probable” only in the roundabout sense that it gives the highest 
probability to the observed results. But they would not reject the m.p.n. 
procedure itself, which has come to be regarded as a remarkably reliable 
tool of very wide utility. At the risk of a slight digression it is interesting 
to indicate the reasons for the reputation which the method has acquired. 

The modern approach is to appraise any method of estimation by 
results. For the m.p.n. this is done, ideally, by conducting a large num- 
ber of dilution series with given v’s and n’s, in circumstances where the 
true density is known. For each series the density is estimated by the 
m.p.n., so that we accumulate a large number of observations on the 
amounts by which the m.p.n. is in error. These observations can be 
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summarized conveniently by plotting the frequency distribution of the 
m.p.n. about the true density. If this frequency distribution groups 
very closely about the true density, we know that the estimates are 
usually good. Such a set of experiments would be difficult and expensive 
to conduct, but if we assume that the mathematical analysis which has 
been applied to the dilution method is valid, we can work out the 
frequency distribution by purely mathematical methods. 

As the numbers of samples 7; become large, the frequency distribu- 
tion of such an estimate (m.p.n. or other) usually tends to assume a 
certain limiting form—the normal distribution. An important general 
result has been established about these limiting distributions (7), to the 
effect that the limiting distribution of the m.p.n. has the smallest 
standard deviation that can be achieved by any method of estimation. 
Roughly speaking, this means that the m.p.n. gives on the average at 
least as precise estimates as any other method used on the same data. 
There is no point in seeking further for a more precise estimate. The 
theorem cannot be proved in general when the numbers of samples are 
small, but experience suggests that the m.p.n. technique is among the 
best methods of estimation in this case also. Consequently the m.p.n. 
method is now generally used in a great variety of problems of statistical 
estimation, though it more frequently goes by the name of the ‘“‘method 
of maximum likelihood.” 


THE PLANNING OF DILUTION SERIES 


In preparation for an estimation by the dilution method, three de- 
cisions must be made: (i) what range is to be covered: i.e. what are to 
be the highest and lowest sample volumes; (ii) what dilution factor is 
to be used; and (iii) how many samples should be taken for each dilution. 

Specific decisions must depend on a knowledge of the limits within 
which the true density is likely to lie and on the precision desired in the 
estimate. The way in which precision is to be measured needs some 
comment. Suppose that the true density is thought to lie somewhere 
between say 2 and 400 organisms per ml. No matter where the true 
density should happen to be within this range, we want to plan the series 
so that the estimate will have a specified “precision.” This might be 
taken to mean that the standard error of the estimated density should be 
say 30 organisms. But this does not seem a reasonable definition of 
“equal precision,” because although an estimate of 360 + 30 organisms 
seems satisfactorily precise, an estimate of 5 + 30 organisms seems very 
imprecise. Instead, we take “equal precision” to imply that the standard 
error bears a constant ratio to the true density, in other words that the 
coefficient of variation of the estimated density is constant. A further 
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potent reason for adopting this concept is that in a well-designed series 
the m.p.n. estimates do have approximately the property that the coeffi- 
cient of variation is independent of the true density. Thus in a sense we 
are making a virtue of necessity. 

The following remarks are intended as a rough guide in the planning 
of dilution series. They were derived from investigations of the precision 
of the m.p.n. 

HIGHEST AND LOWEST SAMPLE VOLUMES 


These are determined by the range of densities with which we expect 
to have to cope. With a single dilution it was mentioned that for the 
best results the expected number of organisms in the sample volume v 
should lie between 1 and 2. It follows that in a series of dilutions the 
expected number of organisms in the highest sample volume vy, should 
be at least 1, otherwise there is a risk that all samples will be sterile. 
Similarly the expected number of organisms in the lowest sample volume 
vz Should not exceed 2, to avoid the risk that all samples will be fertile. 
This line of reasoning would lead to the rule that a dilution series is 
capable of estimating any density that lies between 1/v, and 2/v, . 

This rule is satisfactory if a substantial number of samples, say 
20 or more, are being taken at each dilution. With very small numbers 
of samples per dilution, which are typical in certain lines of work, the 
rule is not quite stringent enough, in that it allows too much risk that 
all samples may be fertile. Suppose that we have three 10-fold dilutions, 
with sample volumes 0.01, 0.1 and 1 ml. This series should be able to 
estimate any true density between 1 and 200 organisms per ml. If, 
however, the density happens to be 200 per ml., so that the expected 
number of organisms per sample in the lowest sample volume is 2, then 
the probability of a sterile sample at this dilution is e *, or 0.135. The 
probability of a fertile sample is 0.865. If only four samples are used per 
dilution, the probability that all four are fertile is (0.865)*, or 0.56. At 
the two higher concentrations, all samples are practically certain to be 
fertile. Thus the worker runs about a 50-50 chance that all his samples 
will be fertile, which usually necessitates repetition of the series. On the 
other hand, with 20 samples per dilution, the probability that all are 
fertile is (0.865)”, or only about 0.05. 

Thus in small experiments it is safer to reduce the upper density 
value from 2/v, to 1/v, . In practice, we use this rule by first guessing 
two limits 6, and 6, between which we are fairly certain that the true 
density lies. The sample volumes are then chosen to satisfy the rules 
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For example, if we are confident that the density lies between 10 and 750 
per ml., the highest sample volume should be at least 1/10, or 0.1 ml. 
The lowest sample volume should not be more than 1/750 ml. The 
three 10-fold dilutions 1/10, 1/100 and 1/1000 ml., or the four 5-fold 
dilutions 1/10, 1/50, 1/250 and 1/1250, would amply cover this range 
of densities. 


THE DILUTION RATIO 


As regards the selection of a dilution ratio, there are two relevant 
results. If the total number of samples in the whole series is kept fixed, 
the average precision is practically the same for any dilution ratio 
between 2 and 10. The advantage of a low dilution ratio, which requires 
more work, is that the precision is more nearly constant throughout the 
range of densities between 1/vq and 1/v,; . These points may be illus- 
trated by a comparison between the dilution ratios 2 and 10, in series 
designed to cover the same range of densities and to use the same total 
number of samples, 72. The details for the two series are as follows. 


Dilution No. of samples Volumes of samples (ml.) 
ratio per dilution 
2 9 .01, .02, .04, .08, .16, 
.32, .64, 1.28 
10 24 {Oth cil), WO 


The two series should cover a range of densities from 1/vyz to 1/v, , or 
from about 1 to 100 organisms per ml. The dilution ratio 2 requires 
eight dilutions, with 9 samples per dilution, whereas the dilution ratio 
10 requires only 3 dilutions and allows 24 samples per dilution. 

In Figure | the standard error of the m.p.n., expressed as a percent 
of the true density, is plotted against the true density (on a log scale). 
With both dilution ratios the standard error per cent is fairly constant 
for any true density between | and 100 organisms per ml. Outside these 
limits the standard error begins to rise steeply, except that with the 10- 
fold series, which has 24 samples per dilution, the rise is postponed until 
6 = 200, for reasons given in the previous section. Inside the limits the 
standard error shows a periodic fluctuation which is noticeable with the 
10-fold dilution but negligible for the 2-fold. With a 5-fold dilution 
(not shown), this periodic effect would be just perceptible. It is present 
with the 10-fold series because practically all the information is con- 
tributed by a single dilution. When the true density is about 1.5 or 15 
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or 150, so that one of the dilutions has about 1.5 organisms per sample, 
there is a trough, with peaks in the intervening densities where no sample 
has a density close to this value. With the 2-fold series, several dilutions 
contribute information and the periodic effect is smoothed out. On the 
whole, the 2-fold dilution gives a slightly lower standard error over the 
range from 1 to 100 organisms per ml., the difference being about 7 per 
cent. For these reasons a low dilution ratio is preferable if the extra 
work involved can be accomplished easily. 


60 
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FIGURE I. COMPARISON OF DILUTION RATIOS 2 AND 10 


The curves in Figure | were calculated by assuming that the formula 
which holds for the standard error in the limiting distribution, appropri- 
ate for very large samples, could be applied to this example in which the 
total number of samples is 72. Some unpublished work by Dr. I. J. 
Bross on the distribution of the m.p.n. in small samples indicates that the 
standard errors are higher than those obtained in this way from the 
limiting distribution. Further, the periodicity with the 10-fold dilution 
does not follow the course predicted for it. However, the two principal 
conclusions from Figure 1 still appear to hold in small samples, namely 
that the standard error is more stable with a low dilution ratio, and also 
tends to be slightly lower. * 


*This work was carried out under contract with the Office of Naval Research. 
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STANDARD ERROR OF THE M.P.N. 


In many types of investigation there may be only a few samples for 
each dilution. In this event the distribution of the estimated density d 
is very skew, and to attach a standard error to d is misleading. The 
distribution of log d is more nearly symmetrical, and it is recommended 
that tests of significance and the construction of confidence limits be 
performed from log d rather than from d. If there are n samples per 
dilution (assumed the same in all dilutions), the standard error of logio d 


may be taken as 
Wea eo 
n 


where a is the dilution ratio. This formula can be used for any density 
which lies between 1/vz and 1/v, , and for any dilution ratio of 5 or less. 
For a dilution ratio of 10, a more conservative factor of 0.58 is preferable 
to 0.55, to allow for the contingency that the estimation may have been 
made at a point where the standard error has one of its peaks. Thus 
for dilution ratio 10 the formula becomes simply 0.58/+/n. Note that 
the formula does not explicitly involve the number of dilutions used. 

To test the significance of the difference between two estimated 
densities, made from independent series, we compute 


log d; — log d; 
0.55 ee ay es log a, 
1 


No 


and refer to the normal probability tables. 

The construction of confidence limits may be illustrated by assuming 
that we have three 10-fold dilutions, with 5 samples per dilution. The 
standard error of log d is 0.58/+/5, or 0.259, so that the 95 per cent 
confidence limits for log d are (log d + 0.518). It follows that to get the 
upper confidence limit for d, we must multiply d by antilog (0.518) or 
3.3, and to get the lower confidence limit we must divide d by 3.3. 

For the common dilution ratios, 2, 4, 5, and 10, Table I shows the 
standard error of log d for any number of samples per dilution between 
land 10. The table also gives the factor by which the estimated density 
must be multiphed and divided in order to obtain upper and lower 95» 
per cent confidence limits respectively. In the example presented by 
Fisher and Yates (6), the number of rope spore organisms per gram of 
potato flour was estimated to be 760. The dilution ratio was 2 and there 
were 5 tubes per dilution. From Table I, the factor for n = 5, a = 2 is 
1.86. Hence the upper confidence limit is 760 X 1.86 or 1414, while the 
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TABLE I 
STANDARD ERROR OF LOG d AND FACTOR FOR CONFIDENCE LIMITS 


: Factor for 95% 
No. of Sree sc 2) confidence limits 
samples 
per dil. Dilution ratio (a) Dilution ratio (a) 
n 2 4 5 10 2 4 5 10 
1 O01 | .427 . 460 . 580 4.00 7.14 8.32 | 14.45 
2 SPAN | Gage tO .325 .410 2.67 4.00 4.47 6.61 
3 .174 246 265 330 223 eo OM tones 4.68 
4 .150 214 . 230 290 2.00 2.68 | 2.88 | 3.80 
5 185 Eton 206 . 259 L867, 2:41 2.58 | 3.30 
6 123 | 2174 188 287 1776) 2523 2.38 | 2.98 
7 114} .161 174 .219 1.69 2.10 2.23 2.74: 
8 ay) clot | “2i6sa £205 |) 1.64,089.004] 2.19 1 9187 
9 . 100 .142 .153 .193 1.58 1.92 2.02 2.43 
10 .095 185 .145 .183 il FHS) 1.86 1595) 21732 


lower limit is 760/1.86 or 409. This factor clearly fulfills the same 
general purpose as would a standard error, if it had been appropriate to 
attach one to d. 

The table makes it evident that the dilution method is of low pre- 
cision, as is to be expected from a method that does not use direct 
counts. Large numbers of samples must be taken at each dilution if a 
really precise result is wanted. Further, the table is likely to over- 
estimate the accuracy of the method, since it is derived on the assumption 
that the mathematical analysis corresponds exactly to the practical 
situation. With a large volume of liquid that cannot be mixed, the 
distribution of organisms may be far from homogeneous. The method 
will determine the density in that part of the liquid from which the 
initial sample was taken. This might be very different from the average 
density over the whole liquid, and this source of error could be more 
important than the error in the dilution method itself. 


SUMMARY OF STEPS IN PLANNING 


The decisions to be made involve a choice of the dilution ratio, a, 
the number of dilutions and the actual sample volume in each dilution, 
and finally the number of samples n to be used at each dilution. The 
steps may be set out as follows. 

1. Decide on the limits 6, and 6, within which the true density 


appears certain to lie. 
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2. Calculate the lowest and highest sample volumes by means of 

the relations 
ai wane 
CS ae On 

3. Select a dilution ratio. A low ratio is preferable whenever 
feasible. 

4. The number of dilutions and the actual volumes for each dilution 
may now be chosen so as to satisfy the requirements that the highest 
sample volume must not be less than vy, and the lowest must not 
exceed vz, . 

5. The precision to be expected for any specified number 7 of samples 
per dilution may be appraised from Table I, if the number of samples 
per dilution is less than 10, or from the formula for S.E. cog a) . Choose 
the number of samples in the light of the precision that is desirable 
and the amount of work that it is practicable to do. 
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STATISTICAL ANALYSIS IN SANITARY ENGINEERING 
LABORATORY STUDIES* 


Rour Evtassen, Sc.D. 
Professor of Sanitary Engineering 
Massachusetts Institute of Technology 
Cambridge 39, Mass. 


HE APPROACH of the sanitary engineer to statistics is that of one who 

is seeking an essential tool for the design of laboratory experiments 
and the interpretation of their results, and for the planning of the 
organized collection of field data and the analysis of the significance of 
this data. The principal distinction between statistical operations on 
laboratory and field data is that of using small numbers of samples in the 
laboratory, usually with a few controlled variables, as opposed to the 
use of very large samples in the field, many times with an appreciable 
number of uncontrolled variables. This paper will discuss the statistical 
methods which the writer and his colleagues have used in the planning 
and execution of laboratory studies. 

A sanitary engineering research laboratory is concerned with experi- 
mental work involving physical, chemical and biological principles as 
applied to the treatment of water, sewage and industrial wastes, and the 
contamination of milk, foods, atmospheres and bodies of water. Many 
of the phenomena under investigation are biochemical in nature. The 
vagaries of biochemical reactions, together with the heterogeneity of 
the media employed, frequently lead to experimental results which may 
follow a general trend but vary with individual observations. Under 
these conditions it is obvious that the sanitary engineer and his associates 
in the fields of chemistry and biology have a great need for the applica- 
tion of statistics in the planning of experiments and evaluation of their 
results. 

In addition to a thorough knowledge of their own field of research, 
sanitary scientists ought to be conversant with the basic concepts of 
statistics. They should possess the ability to choose the most effective 
statistical procedure for the analysis of the data at hand, to scrutinize 


*Presented before the Engineering, Laboratory and Statistics Sections at 77th Annual Meeting of 
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the fundamental assumptions underlying the particular statistical pro- 
cedure chosen, and to recognize by testing to determine whether these 
assumptions are fulfilled by the situation surrounding the data. Fre- 
quently they may have to call upon the services of mathematical 
statisticians to assist in this work. Statistics then becomes a scientific 
tool for testing hypotheses and estimating the significance of experi- 
mental results obtained in the sanitary science research laboratory. 

In the planning of experiments in new fields of sanitary engineering 
research, it is frequently necessary to develop new analytical techniques. 
Such a situation arose in connection with a study of the effect of chlorin- 
ated hydrocarbons on the suppression of hydrogen sulfide production in 
sewage (1). The problem involved the sweeping out of hydrogen sulfide 
from an anaerobic culture medium by means of an inert gas and trapping 
the hydrogen sulfide for quantitative analysis with the complete exclusion 
of air. 

Repeated trials of the technique, together with necessary revisions, 
were observed by means of statistical tests, particularly standard devia- 
tions and their coefficients of variation. Starting with coefficients as 
high as 30 per cent, controls were effected until the results in Table I 
were obtained. 


TABLE I 


SULFIDES PRODUCED IN ANAEROBIC FERMENTATION OF SEWAGE 
AFTER 42 HOURS AT 37°C. 


Run No. HS Total Sulfides 

ppm ppm 

OR Betas) <ecab hy aie sity teens et ok re 210 239 

Qe Rr ne hates Ses ae tae ery eee 221 215 

See he ames tact AS he ae emcee, AE OL 218 240 

ARO Nae, Ae ls, eae gem ee eee rman: 228 22, 

ite i A TE a en a Lae Cae Lx wea 220 230 

\ehe 7. SAN Mesa emey ene 5 Sie MrCI ec Boa 8 mr 227 226 

ig”, Sas PT aaa Pa Ae AP Me Na 5.” 223 2G 

CA SN ec ee MS hed, Pkt BA tages Bs 1 224 228 

Oe eae Lae ea eck, aa So ee eee ees eS 192 208 

Mean ase ir) 1: ti gee NASB nes rapt) ot mee 218 227 

Standards Deviations (cele ant enn Natal 10.2 

Coefficient of variation (%) ........ fay. 4.5 


These statistical results indicate that two-thirds of the time this 
analytical procedure can be expected to yield quantitative measures 
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which will be within 5 per cent of the mean value which would be ob- 
tained if a number of runs were made. Furthermore, 95 per cent of the 
time the results would be within 10 per cent of the mean. Following 
this, the research work could then be carried out with the assurance 
that experimental errors in analytical techniques had been reduced to a 
satisfactory level. 

A different type of correlation analysis was involved in the develop- 
ment of a rapid method for the determination of dissolved oxygen in sea 
water. When making pollution studies in harbors and estuaries, it is 
preferable that all chemical analyses be made on the sample boat. The 
Winkler test (2) does not lend itself too readily to such a situation. In 
1935 Gilcreas (3) introduced a colorimetric method using Amidol for 
dissolved oxygen determination in ordinary waters. However, applica- 
tion to the analysis of sea water was not feasible because the presence 
of sodium chloride interfered with the test. In an attempt to adapt the 
method to sea water Lieber (4) conducted a laboratory investigation as a 
master’s thesis under the author’s direction at the College of Engineering 
of New York University. The results of many tests with Amidol (diami- 
dophenol hydrochloride) were compared on a statistical basis with those 
using the Winkler Method in waters having a wide range of salinity. 

By the method of least squares the data was fitted to a linear equation 
of the type y = a+ bz. This could be used in graphical form to convert 
the results of field observations with Amidol (y) to equivalent Winkler 
values (x). From the lines of regression for each salinity value the coeffi- 
cients of correlation were obtained. These values are indicated in 
Table II 


TABLE II 


STATISTICAL RESULTS IN COMPARING AMIDOL AND WINKLER METHODS 
FOR DISSOLVED OXYGEN DETERMINATION IN SEA WATER 


Line of Regression StandardError | Coefficient of 
Percent Vso Oe of Estimate Correlation 
Sea Water 
a b Sy r-% 
0 —0.223 1.108 0.176 99.25 
20 —0.338 OU 0.230 99.65 
40 —0.222 1.325 0.262 99.50 
60 —0.123 1.341 0.348 98.99 
80 —0.505 1.473 0.367 99.35 
100 —0.386 1.536 0.258 99.65 
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Applying the method of least squares, or the simpler method of 
graphical analysis, it would be desirable to express the parameters a 
and b as linear functions of percent salinity. These relationships have 
been evaluated approximately by the graphical method: 


a = —0.250 — .001 (% salinity) 
b = 1.100 + 0.0045 (% salinity) 


For a complete analysis, these values of a and b would have to be tested 
for correlation with observed data and the constants in the above equa- 
tions revised by trial and error to the degree of accuracy required. 

A remarkably high degree of correlation was evident in the fit of the 
equations determined by least squares and the observed data. For 
example, the Sy for 60 percent salinity was 0.348. This indicates that 
the method will yield results within 0.35 ppm about 68 per cent of the 
time or within 0.70 ppm about 95 per cent of the time. Thus, the Amidol 
method may be relied upon to yield results with a good degree of re- 
liability, at the same time providing a highly useful tool for pollution 
studies. 

In bacteriological studies the fitting of curves to scatter diagrams 
is particularly useful for the interpretation of results. On a sponsored 
research project at New York University, Krieger (5) studied the effect 
of various chlorine contact periods on the rate of killing of coliform 
organisms. The observed data indicated that the phenomena could best 
be interpreted on the basis of a curve of non-linear regression in the 
form of 


1 
oes a+ bx ; 
TABLE IIT 


CURVILINEAR REGRESSION CHARACTERISTICS OF M.P.N. VS. 
CHLORINE RESIDUAL EXPERIMENTS 


Square of Square of 1 & 
Contact Time Standard Standard SES a+ be 
Minutes Error Deviation 
(Sy)? (o,)? a b 
5 0.192 0.240 (0). 0.32 
10 0.649 0.899 OR2 iI 0.74 
15 0.672 1.180 | 0.28 1.04 
20 ets? poe 0.65 1.02 
30 iL aay 1.640 1.91 0.47 
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The M.P.N. per 100 ml. were represented by y and the chlorine residual 
(R) after specified contact periods by x in the general equation: 


c k 
M.P.N. = exp (- ae Ih 

By the method of least squares, curves were developed for the various 
contact periods employed in the experiments. The statistical character- 
istics of the curves are shown in Table IIT. 

Plotting the values corresponding to the derived equations gave the 
family of curves shown in Fig. 1. By means of curves of a similar type 
derived for any particular sewage, state departments of health, designing 
engineers and plant operators can be guided in the design of chlorine 
contact chambers and the proper chlorine dosages to achieve the degree 
of disinfection of the sewage required for the receiving stream. 

Laboratory studies of sewage oxidation and stream degradation re- 
quire extensive use of the biochemical oxygen demand (B.O.D.) test. 
In spite of the vagaries of biochemical reactions, particularly when used 
for quantitative measures, the results can be interpreted successfully by . 
conducting a sufficient number of experiments and subjecting the results 
to statistical analysis. 

It has been generally agreed that the first stage biochemical oxidation 
of sewage is a unimolecular phenomenon proceeding in accordance with 
the relation y = L(1 — 10°“). L represents the ultimate first stage 
B.O.D. (ppm), y is the B.O.D. after time ¢ in days, and k is the reaction 
rate constant. For research workers the determination of the param- 
eters k and Z is extremely important. Various statistical procedures 
may be applied to the observed data. The most accurate method would 
be by least squares, but for an equation of this type the computations 
would prove too cumbersome for large amounts of data. The most 
convenient method, and one which has proven very accurate, is that 
suggested by H. A. Thomas (6) and known as the Method of Moments. 

The statistical hypothesis employed is that the observed data can 
be fitted to a unimolecular curve by assuming that the zeroth moment 
and the first moment for the data equal the same moments for the fitted 
theoretical curve. The type of curve and the equations of moments 
are shown in Fig. 2. Evaluation of the summations for all of the data 
may be facilitated by means of graphs prepared by Professor Thomas for 
specific sequences of days of B.O.D. measurements. By the summation 
of 2 y; and » y,t; from the observed data, the parameters k and L 
may be read from the graphs. 

This method has proved to be a useful research tool in an extended 
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XE 
y= L100) 


O = Observed 8.0.0. (ppm) 
at time * 


L = Ultimate B.0.0. (ppm) 


Tice G.aly/s 


HYPOTHESES 


Zeroth Moments yD Sen eT 


i=0 


First Moments S ty. = > ty 
1=0 1=0 


FIGURE 2. METHOD OF MOMENTS FOR FITTING UNIMOLECULAR CURVE TO 
OBSERVED B.O.D. DATA 


series of experiments involving a study of the effect of industrial wastes 
on the oxidation of sewage and organic matter in streams. 

The parameters k and L can be made the basis for further statistical 
operations in evaluating the significance of results and the difference 
between means resulting from radioactivity. 

In the course of the radioactivity studies at the College of Engineering 
in New York University, experiments were made to determine whether 
sewage taken from a specific manhole at the same time on Tuesday and 
Wednesday could be considered as coming from the same population. 
Since the ultimate effect of radioactivity would be measured by its 
effect on the rate factor, k, statistical tests were made on the k values to 
determine whether there was any significance in the difference between 
the means. 

An hypothesis was set up that the sewage taken on Tuesday and 
Wednesday would yield no significant difference in k values. As shown | 
in Table IV, ¢ was determined from the computed values of k. Using 
Student’s ¢ distribution, with six degrees of freedom, P was found to be 
equal to .005. This signifies that one could expect as great, or greater, 
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difference between the k’s for sewage sampled on Tuesday and Wednes- 
day only 5 times in 1000 due to chance alone. Therefore, the hypothesis 


was rejected. 


TABLE IV 
t-TEST FOR DIFFERENCE BETWEEN k VALUES OF SEWAGE SAMPLED ON 


TUESDAY 


AND WEDNESDAY 


k Values ‘. 4 
d = kw = kp d = ol (d = d)? 
Tues. Wed 
. 139 .155 + .016 — .015 000225 
.128 w74: + .046 + .015 000225 
124 .158 + .034 + .003 000009 
. 168 . 154 — .014 — .045 002025 
.129 ies + .044 + .013 000169 
3 Be pla + .050 + .019 .000361 
los 2, + .039 + .008 .000064 
2di= 215 x(d — d)? = .003078 
a= (O81 
ys we, “\2 
o, = ee Oe Reuss = .02265 
i = Al 6 
Sinbrtinwordin coe = EDS ages 
v/n V/7 
d— mm _ 0 : 
t= ( De Ua 3.6 Hypothesis: m = 0 


00s 


or 


From Student’s ¢ Distribution, 


Degrees of Freedom = 6 


P= .005 


On the basis of these statistical results, the experiments which fol- 


lowed had to be planned to be 


conducted only on sewage taken at a 


certain hour on a specific day for each successive run. 

Further planning of the radioactivity studies involved the determina- 
tion of the number of successive experiments to be conducted under 
similar conditions in order to arrive at mean values of k within the limits 
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of experimental error. On the basis of experience with the B.O.D. test, 
it was decided that a standard error of the mean of less than 0.008 
(5% of mean k) would be satisfactory. Ten runs were made, with values 
of k obtained as noted in Table V. Computations indicated that a 


TABLE V 
NUMBER OF RUNS NECESSARY TO ASSURE MEAN k VALUES WITH A 
STANDARD ERROR OF LESS THAN .008 


k bk (k — k)? 
1240 — 0252 0006 
1275 — .0217 0005 
1295 —.0197 0004 
1375 = O11; 0002 
.1390 — 0102 0001 
1475 —.0017 0000 
1510 + 0018 .0000 

1680 +.0188 0004 
1705 + .0213 0005 
1975 + .0483 0023 

1.4920 0050 


k = 0.1492 x(k — k)? = .005 


_ «(ak — ky? 005 _ a 
7 ee = 4/9 = 2.4 X 10 
oe ee =e hie 


» > /N~ 4/10 


standard error of 0.0076 could be expected with ten separate experi- 
ments. Therefore, subsequent experimental work was planned on the 
basis of ten B.O.D. runs for each concentration of radioisotope in the 
sewage. 

Many other examples might be presented to show the need for and 
use of statistical methods in sanitary engineering laboratory studies. 
The author has attempted to show some procedures which have been 
employed in laboratory work conducted under his supervision. As 
research workers become more familiar with statistical methods of 
analysis, they will recognize many places where statistics can be of great 
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assistance in the planning, execution and interpretation of laboratory 
work. 

Its most valuable place will be found in the interrelation of the oft- 
times scattered, but nevertheless consistent, data which frequently is 
obtained in the study of biochemical reactions such as the disinfection 
of water and sewage, the oxidation of organic matter in sewage treatment 
processes and in the self-purification of streams, and in anaerobic fer- 
mentation processes for sludges resulting from the treatment of sewage 
and industrial wastes. 
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THE APPLICATION OF STATISTICAL TECHNIQUES 
TO SEWAGE TREATMENT PROCESSES* 


JoHN W. FERTIG 
Columbia University 
and 
Austin N. HELLER 
The Barrett Division, Allied Chemical & Dye Corporation, New York 


Re OF THE inherent variability in the characteristics of sewage 
and because of the shortcomings of many of the primary analytical 
methods used, it is often impossible to determine the specific effectiveness 
of a given treatment process by examining a single or even a few samples 
of sewage. It is often necessary to perform many experiments. Because 
the resultant data may vary considerably, it becomes necessary to apply 
statistical methods for their evaluation. 

There are many phases of sewage treatment processes which can 
profit from the use of statistical methods. We have chosen to illustrate 
some of the techniques in terms of the problem of disinfection of raw 
sewage and sewage effluents. 

The results of an extensive series of experiments designed to measure 
some of the chlorination aspects have been reported by Eliassen, Heller 
and Krieger (1). It is not our present purpose to review this work, but 
rather to take extracts from the original experimental data to illustrate 
some rather fundamental statistical techniques. We shall also draw on 
hitherto unpublished data referring to the effect of mechanical mixing 
on M.P.N. of coliform bacteria of raw sewage. 

We should like to emphasize that throughout the presentation the 
stress shall be on statistical methods rather than on these data. We do 
not purport to believe that the results quoted here will have any great 
bearing on the operation of sewage treatment plants. We do believe, 
however, that some of the considerations we raise will have merit in the 
design of experiments in permitting a more accurate evaluation of treat- 


ment processes. 
DETERMINATION OF SANITARY QUALITY 


Sanitary quality is often measured by the count of coliform organ- 
isms per ml., as estimated by the “most probable number” (M.P.N.) 


*Presented at the Annual Meeting of the American Public Health Association in New York 
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technique. The M.P.N. of the chlorinated sewage was determined by the 
partially confirmed test using brilliant green lactose bile broth as the 
confirmatory medium. For the unchlorinated sewage, direct planting 
was made into this medium. All M.P.N. determinations were made on 
the basis of four ten-fold dilutions using three tubes for each dilution. 

Although four dilutions were used, the M.P.N:’s in most instances 
were determined completely (to two significant figures) by the three 
critical dilutions. Consequently most of the four-dilution codes were 
translated into M.P.N. by Hoskins’ three-dilution tables (2). In general, 
it is not difficult to select the three critical dilutions. Thus, if the code 
(for positive tubes) is 3 3 2 0, one selects 3 2 0 to enter the tables. In 
some codes, it is not so clear which are the three critical dilutions. Thus, 
given 3 3 1 1, the proper choice is 3 1 1, which gives the same result as 
the complete four-dilution code. Occasionally the M.P.N. read from the 
complete four-dilution code differs in the second significant figure from 
the value appropriate to the three-dilution code. Thus, 3 1 1 0 gives an 
M.P.N. of 74, whereas 3 1 1 gives 75 and 1 1 0 gives 73. 

Where the code 0 0 0 0 occurred, the convention was followed of 
assigning the lowest value corresponding to one positive tube. Where 
the code 3 3 3 38 occurred, it was assigned a value corresponding to no 
positive tubes at the next higher dilution. Other conventions would 
result in assigning other values to these “indeterminate” codes. How- 
ever, with a well-chosen system of dilutions, such codes will not arise 
very often. Furthermore, if there are a large proportion of such in- 
determinate codes, the data cannot well be analyzed as a series of 
measurements. 


ACCURACY OF M.P.N.’S BASED ON THREE TUBES AT EACH OF 
THREE TEN-FOLD DILUTIONS 

For a given true density of organisms, the M.P.N. obtained from sam- 
ples may vary considerably through the errors of chance. This variation 
was studied for the case of three tubes at each of three dilutions. The 
probability distribution of M.P.N. corresponding to a number of true 
densities was calculated directly according to the fundamental formula 
(2). Relatively few of the 64 possible codes occur with any appreciable 
frequency for any given true density, so that the calculations were not 
unduly laborious. 

The distribution is discrete and very irregular. The terms are not 
spaced uniformly and do not increase regularly as the true density is 
approached from either higher or lower sample values. As Halvorson 
and Ziegler (3) pointed out for the case of 5 or more tubes, the mean 
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exceeds the true density and the standard deviation increases with the 
true density. 

The distribution of log M.P.N. is more regular, however. While the 
distribution is discrete and the terms do not increase regularly as the 
center is approached, it is more symmetrical than the distribution of 
M.P.N. The mean of the logs is slightly larger than the log of the true 
density, but much less so than in the case of M.P.N. The standard 
deviation of log M.P.N. is remarkably constant, varying from 0.29 to 
0.33 for the range of true densities tried (0.1 to 5.0). The value of the 
standard deviation is in fact very close to the value of 0.30 given by 
Fisher’s method of determining the standard error of maximum likeli- 
hood estimates (4). This close agreement is surprising since Fisher’s 
method is strictly applicable only for the case where the number of tubes 
at each dilution is large. 

The distribution of M.P.N.’s based on three tubes at each of four 
dilutions was also investigated for several true densities. It is only 
shghtly different from that described above for three dilutions. 

From the above discussion of the distribution it follows that an 
M.P.N. based on three tubes at each of three (or four) dilutions may 
differ considerably from the true value. If we may regard the distribu- 
tion of log M.P.N. as something like a normal curve with three standard 
deviations as the limit of variation, then an observed M.P.N. may 
apparently vary from 1/10 to 10 times the true value due to the errors 
of chance. Two observed M.P.N.’s would have to differ enormously to 
indicate a significant difference. With a series of experiments, however, 
one compares the average M.P.N.’s and these have much greater stability 
than a single M.P.N. 


DAILY VARIATION OF OBSERVED M.P.N.’S 


In the latter part of 1946, a series of experiments was run on Bronx 
sewage with chlorine demands ranging from 1.2 to 3.6 p.p.m. The 
coliform count was determined under various degrees of chlorination and 
mixing of the sewage and chlorine solution (1). All of the M.P.N. deter- 
minations were made in triplicate. In the table below are given the 
data for nine different days under conditions of no chlorination and no 
mixing. 

It is immediately apparent that the days vary widely from each other, 
so that the above set of 27 values cannot be considered as homogeneous. 
On each day, however, the triplicates should exhibit merely random 
variation. One may compute a measure of this random variation by 
averaging the nine daily standard deviations, each of which has two 
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TABLE 1 


TRIPLICATE VALUES OF LOG M.P.N. ON NINE DIFFERENT DAYS, 
UNCHLORINATED SEWAGE, NO MIXING 


Day Log M.P.N. (=2) Total Average 
1 2.97, 3.08, 2.97 9.02 3.01 
2 4.04, 4.38, 3.63 12.05 4,02 
3 1.56, 2.15, 2.386 6.07 2.02 
4 2.63, 3.36, 3.36 9.35 3.12 
5 2.97, 3.18, 2.63 8.78 2.93 
6 DA, Gin (By Be 10.35 3.45 
7 3.63, 3.36, 3.36 10.35 3.45 
8 3.36, 3.36, 3.63 10.35 3.45 
9 2.63, 3.36, 2.36 8.35 2.78 

Total 84.67 3.14 


degrees of freedom. This average has eighteen degrees of freedom and 
may conveniently be computed as 


; E » _ (9.02)? + (12.05)? + +++ + (8.35)? 
—— POM = 3 


where Sa” is the sum of the squares of the individual values (274.8275). 
¢ = 1/7 0.1014 = 0.82 


This standard deviation agrees remarkably well with the measure of 
random variation deduced theoretically (0.29-0.33). 

Utilizing the triplicates over a number of other situations—100% 
chlorination, mixing, etc.—we obtain an estimate of random variation of 
0.33, based on 108 degrees of freedom. 

From the data above we may also obtain a measure of the variation 
between the daily averages based on eight degrees of freedom. The 
variation between daily averages (multiplied by 3) is readily obtained as 


ES + (12.05)? + -++ + (8.35)? _ (84.67)? 
3 27 


| + 18 = 0.1014 


| + 8 = 0.9355 


A comparison of this variation with a measure of random variation 
will test whether the days do in fact have significantly different M.P.N.’s 


on the average. The comparison is made by the variance ratio or F 
test (5) 


F = 0.9355/0.1014 = 9.2 
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For 8 degrees of freedom in the numerator and 18 in the denominator, 
an F of 3.7 is exceeded only 1% of the time when days are the same. 
Consequently the observed day to day variation appears to be highly 
significant. Of course, it is well known that the sewage characteristics 
throughout the course of a day are just as variable or more so than those 
from day to day. 

Because of the significant variation from one type of sewage to 
another arising from sampling different days, times or places, it is evi- 
dently necessary that any comparison of treatments be so arranged that 
each treatment is represented with each type of sewage sample. 


COMPARISON OF TWO DIFFERENT TREATMENTS 


On 10 different days during the early part of 1947, a series of experi- 
ments was run on Bronx sewage in order to compare the bactericidal 
effect of various methods of mixing the chlorine solution and the sewage. 
For our purpose we shall take the data comparing a rapid mix method 
(initial rapid mixing of 15 sec.) and a no mix method at 100% chlorina- 
tion. The samples for determining the M.P.N. were taken after a contact 
period of 10 minutes with the chlorine solution. 


TABLE 2 


COMPARISON OF RAPID MIXING WITH NO MIXING 
OF SEWAGE CHLORINATED 100% 


No Mix — 
Day Rapid Mix No Mix Rapid Mix (=A) 

1 3.04 3.04 0.00 
2 2.38 Babies +1.00 
3 2.38 3.38 +1.00 
4 1.97 2.38 +0.41 
5 1.36 3.04 +1.68 
6 0.97 1.36 +0.39 
76 1.97 2.18 +0.21 
8 1.63 2.66 +1.03 
9 0.15 2.66 +2.51 
10 1.63 3.04 +1.41 
Total 17.48 27712 +9, 64 
Average 1 2.71 +0.96 

(= bie, AY) 


*M_P.N. of 2400 arbitrarily assigned to code 3 3 3 3. 
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In this particular case the superiority of one method over the other 
is evident immediately since for the 9 days where a difference is present, 
the “no mix” has the higher value. The probability that 9 differences 
should have the same sign in the absence of a real factor is 1/256. 

Although most of the differences are quite large, only a few of them 
would be statistically significant per se. The difference between two log 
M.P.N. values is subject to random variation measured by a standard 
deviation of 


o = 0.33+/2 =.0.47 


Only three differences are as large as three times this standard deviation. 

When testing the significance of the difference between two averages, 
it is customary to compare the difference with the standard error of the 
difference. This latter quantity is usually evaluated from the standard 
errors of the two averages. This is not the proper procedure in this case, 
however. Since the days differ significantly from each other and each 
method is represented on each day, one must work directly with the 
daily differences in the two methods. 

The standard deviation of the 10 differences is computed as 


ee Aaa Ga = WARE Oil ROD = 8 


The standard error of the average difference is then computed as 
Oye = 0) \/ 10 =" 0:73) \/ 10 == 0.23 


The average difference of 0.96 between the two methods is thus very 
significant since it is about 4 standard errors distant from zero. Such an 
occurrence is very improbable according to the normal curve, if zero is 
indeed the true difference.* 

The standard deviation of the ten daily differences (0.73) arises from 
the fact that the two methods do not portray the same difference on 
each day. In fact, the obtained value should theoretically equal 
0.33+/2 = 0.47. The difference is just barely significant. This indi- 
cates that the difference between methods varies from one type of 
sewage to another. To generalize concerning the superiority of one 
method over another, it is evidently necessary to make the comparison 
over a range of different sewages. It is perhaps not necessary to make 
the experiments on different days, but only with different sewages. 

On the basis of a larger series of experiments, Eliassen et al (1) report 


*The use of the t-test instead of the normal curve would rarely alter the decision 
as to significance. 
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an even larger difference in average M.P.N. between Rapid Mixing and 
No Mixing than that reported here. 


EFFECT OF PROLONGED RAPID MIXING ON M.P.N. OF UNCHLORINATED SEWAGE 


In none of the experiments reported by Eliassen et al (1) was there a 
significant difference between the M.P.N. for rapidly mixed and unmixed 
non-chlorinated sewage. In an unreported series of experiments on 
Bronx sewage designed to test the effect of rapid mixing for various 
periods of time, there is a small but statistically significant increase in the 
average M.P.N. as the time of mixing increases from 0 sec. to 120 sec. 
The average log M.P.N. increased from 3.66 to 4.14, corresponding to an 
increase in M.P.N. from 4600 to 14,000. 

The type of data is similar to that presented in the preceding section 
of this paper, except that more than two averages are being compared. 
The technique of statistical analysis used was an extension of that dis- 
cussed above, known as the analysis of variance (5). 


EFFECT OF CONTACT PERIOD OF CHLORINE SOLUTION ON BACTERIOLOGICAL KILL 


A series of experiments was performed during 1946 to study the rela- 
tion between bacteriological kill and the period during which the chlorine 
had been in contact with the sewage. The sewage used in these studies 
was collected mostly from Ridgewood, N. J., and from the Bronx, the 
chlorine demand ranging from 1 to 6 p.p.m. 

There was a rather high proportion of indeterminate results in this 
study, i.e., results in which the tubes were either all positive or all nega- 
tive. This was due to a poor choice of dilutions, usually at the shorter 
contact times. Consequently, the data cannot be analyzed as a series 
of measurements as in the previous sections of this paper. For this 
reason, the comparison of contact periods will be made in terms of the 
proportion of samples with M.P.N. less than some specified number. 
A value frequently cited is 3 per ml. and this will be used as a dividing 
point for our purpose. 

Taking the data referring to 100% chlorination, we have 27 samples 
at each contact time. 

. To decide if the proportion of samples with M.P.N. less than 3 varies 
significantly with contact time, we compute chi-square 


v= z| =m] = 19.4 


where 0 refers respectively to the six observed frequencies and T refers to 
the respective theoretical frequencies deduced on the basis of the over- 
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TABLE 3 
DISTRIBUTION OF 81 SAMPLES OF SEWAGE CHLORINATED 100% 
ACCORDING TO CONTACT PERIOD AND VALUE OF M.P.N. 


Samples with M.P.N. 


—— 


Contact Time (min.) Less than 3 3 or more Total 
3 u 20 27 
10 13 14 i 
30 23 4 27 
Total 43 38 81 


all proportion of samples under 3, i.e., 43/81. Thus, at each contact 
time, the theoretical frequency under 3 is 


43/81 X 27 = 143 


and consequently the theoretical frequency of 3 or more is 12.7. 

There are two degrees of freedom for x” since three proportions are 
being compared. From tables (5) of x° we find that a value of 9.2 is 
exceeded only 1% of the time when contact periods are the same. 
Consequently, our value of 19.4 is very significant. 

For 120% chlorination we find the same degree of significance, while 
for 40% and 70% the values are not significant, although the proportions 
vary in the same way with contact time. 

In each case the proportion at 10 min. is slightly higher than that 
at 3 min., but it cannot be considered significant, even if all the data are 
considered collectively. 

The x” method of treating the data illustrated above appears to be 
very arbitrary in that some division points must be used. Actually 
several division points were tried without appreciably altering the re- 
sults. It should be emphasized, however, that if the data can be treated 
as measurements, such treatment is preferable to the y° method. 


USE OF STATISTICAL METHODS IN STUDY OF OTHER TREATMENT PROCESSES 


Some of the methods outlined above should be particularly useful in 
the evaluation of secondary sewage treatment processes, and modifica- 
tions thereof. In addition, the effect of industrial wastes on biological 
processes, both aerobic and anaerobic, can probably be more accurately 
and readily established if requisite experiments are designed to permit 
statistical evaluation. 
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SUMMARY 


Various statistical techniques useful in the study of sewage treatment 
processes have been illustrated in terms of data on M.P.N. of coliform 
bacteria, particularly as influenced by various phases of chlorination. 

(1) Most probable numbers determined from 3 tubes at each of 4 
dilutions are subject to large chance fluctuations. Their logarithm 
exhibits a rather constant variability of approximately 0.30, regardless 
of the true density of organisms. 

(2) The M.P.N. of sewage studied here varies significantly from day 
to day, as shown by comparing the variation between days with random 
variation. 

(3) Rapid mixing of sewage chlorinated 100% resulted in a much 
lower M.P.N. on the average than when no mixing is employed, as shown 
by comparing the average difference with the standard error of the 
average difference. The difference varies from day to day, suggesting 
that a comparison of methods should be made over a variety of sewages. 

(4) Rapid mixing of unchlorinated sewage resulted in a slight but 
statistically significant increase in the average M.P.N. as the mixing time 
increased, as shown by an application of the analysis of variance. 

(5) At 100% and 120% chlorination the proportion of samples with 
M.P.N. less than 3 per ml. increases significantly with the contact period 
of the chlorine solution, as shown by the x’ test. While the results were 
in the same direction at 40% and 70% chlorination, significance could 
not be established. 
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FIDUCIAL INTERVALS FOR VARIANCE COMPONENTS 


Irwin Bross 


Department of Biostatistics 
The Johns Hopkins University* 


NATURE OF THE PROBLEM 


Rees workers frequently have occasion to estimate variance 
components, especially in genetic and sampling problems. For 
example, a 1940 Iowa AAA corn acreage study where two sections were 
selected from each of 1617 townships has the following analysis of 
variance of the corn acreage per section: 


Source of Variation df. Mean Square] Expected Mean Square 
Between Townships 1616 6511.9 ao? + 2o;? 
Within Townships 1617 1954.3 o 


An unbiased estimate of o; is 


EOS 1h. 105408 
oF = 2 


= 2278.8. 


A symbolic representation of the above table would be 


Source of Variation d.f. Mean Square| Expected Mean Square 
1 Mm Vi oy, = a2 + acy? 
2 Ne V2 a? 


Now the cautious research worker would like to have fiducial or 


*This work was done at North Carolina State College, Raleigh, North Carolina in partial satis- 
faction of the requirements for the degree of Master of Science in Experimental Statisties, 
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confidence limits to place on this estimate. The standard error of 3? 
may be estimated by 
a Vn, Ny 


and when n,; and n, are large, the distribution of a? tends to be normal 
so that fiducial limits can be set in the usual way (1). When n, and n, 
are small, however, the distribution of the estimate departs considerably 
from normal, so it would be desirable to have a more refined procedure 
for small n, and nz . 

Satterthwaite (2) has examined the distribution of a,v, + av. and 
suggests that it be approximated by a type III curve with effective 
degrees of freedom 


= (aw, + A>) 
as a7 Pe ae 
Qi , x5 


Ny Ne 


Ne 


He warns that if, as in the case of a; , the a’s are negative, then caution 
must be exercised in using the approximation. This is evident since 
negative estimates of o; can occur which are incompatible with a type 


IIT curve. 
This paper will present an approximate fiducial interval (which 
appears to be useful for n, as low as six) which is based on R. A. Fisher’s 


approach to the problem. 


APPROXIMATE FIDUCIAL LIMITS 


R. A. Fisher (3) approaches the problem from a different viewpoint. 
Let 


2 2 
Ue Sy is T= Of aa Og 
2 
Vy = O71 
a. oF a ea. 
Vo F2 
then 
uP U 


ee ats rl 
Now let x; and 2, be variables distributed as chi-square with n, and nz 
d.f. respectively. Then: 


2 2 
me L104 ee X02 fh VN, abs VoN2 
"My 7 Ne eZ) 


and if we take 
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Dri VoNy Fn, _ Ne 
1; Tt Gi ab 1 Lo 
u (3s) iM == jl 


then 


us il Fn, ne | z Io} 
PD < ly) = Pgh ( i3/ an : 


may be found by integrating out «, and a with the region of integration 
determined by the equation in the parenthesis. If we let 


Wa boMld Te 
FF - DL, +2 
Le 


7] 


the region of integration is: 


Oto © 
25%,<0 
so that 
(1.01) PUL < In) = | fled) dra {| fles) de 


which may be regarded as the basic equation. 
Now if we can find an Z, such that P(L < L,) = a we may obtain 
our fiducial limits immediately since: 


P(L < Ly) 


I 


(2 < In) = P(T < Leu) = Place] < Lo; — v,)} 


Pig Lea) 


Direct integration is not very easy, while tabulation requires four 
entries, (”, , %2 , F, a), and involves interpolation, so that it seems more 
practical to look for an approximate solution. 


The basic equation (1.01) may be solved directly in three cases. 
If n. becomes infinite, then 


No A Fn, 
ee tt yr sey, 


so that the second integral does not involve x, and 


Pe a ( | es as)( i ie dz) ts 
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Hence, if x4 is the value in the chi-square table, (n, , a), then 
z= Xa 
F 
In MAGS x 


This may also be written (in terms of F-table entries) as 


7mm 
1 
eae 
where F,, is the entry in the F-table for d.f. (n, , ©) 
A second special case is for F —o where 


limz = = 


Foo In 
and hence L, = 1/F,. Finally if F = F, , then LZ, = 0 is the solution 
of (1.01). Forif ZL, = 0, 


P(L <0) = py (Fats =f) 1x of = pp Pats 2 ma} 
Zig aid dl == Ml Wy 


iy be 
NoX 
ee) 
Niko 


But (n22,)/(m#2) is distributed as tabular Ff’, F, , and hence 


P(L < 0) = P(F, > F.) = a. 
Now consider the function 


F 
tage 
Yr ee | ieee ce 
ee 
ihe 
T= | eae EL 
Then 
i] 
bina 
Inna 1G, = 
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PF Fo 

ei = ie =e : 

eT Ue ee O 
Fy 


so that L is exact for the limiting cases. 

Other functions beside L have this property and have been investi- 
gated. However L gave the closest agreement with quadrature solu- 
tions for the functions investigated. 

Some typical values are given in Table I which may serve to indicate 
the degree of approximation for the case np = 12 and a = .05. 


TABLE I 
Ste 15 2.0 3.0 5.0 
mM N 
8 
1 02 5.13 5.06 5.05 
3 B17 5.31 5.23 5.10 
8 5.21 5.21 5.19 Bk 


This means, for example, that ifn, = 3 and F/F, = 2, the research 
worker would actually be working at the 5.31% level instead of the 5% 
level if he used the proposed approximation. Discrepancies of this order 
would not ordinarily lead the research worker very far astray. 


APPLICATION OF THE APPROXIMATE FIDUCIAL LIMITS 


The approximate fiducial interval can be obtained by multiplying 
3, by L to obtain the lower limit, and by L to obtain the upper limit 
(central interval). Suppose, for example, that the 5% lower limit is 
desired. Three quantities are required: 


F The ratio v,/v, obtained from the data. 
F os The entry in a 5% F-table for d.f. n, and nz . 
F'os The corresponding entry for d.f. n; and ©. 


If F < Fos, the lower limit will be taken as zero, since o; is a positive 
quantity. If F > Fo; , then 
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The corresponding upper limit 


may be found by using the fact that 


1 


F s(t , M2) = Ea neeen) 
-05 23 1 


i.e., the 95% value of F may be found by entering the 5% table with 
the degrees of freedom interchanged and taking the reciprocal of that 
quantity. Similarly 


1 


Bie = ==) ) 7g ees = 
nt F o5(~, nN) 


The interval between the two points Lo,” and Lo,” constitutes an approxi- 
mate 90% fiducial interval. 


NUMERICAL EXAMPLE 


For purposes of comparison, the numerical example chosen is one 
where current methods might be expected to apply. A portion of an 
analysis of variance table (C/21 * NC7) from H. F. Robinson’s data 
on hybrid corn (5) will be used. 


Source of variation Olin M.S. E(M.S.) 
Females in males in blocks 184 . 0087 oe + 20%? 
Pooled error 234 .0033 oe 


The variance component, a; , is the additional variance among 
paternal half sibs due to female differences. Evidently 


fe — BES = 0027 
0087 _ 


(0025 Ze 


jPes 


Tye al 26 J ee Cea (interpolated). 
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Whence 
2.64 
1.26 _ 1.098 y 
Be Pa Clam, ity es 
ame eG 
and Lo, = .00202. Similarly, 
Times 1 Pot 
Le Tie (Garey ey 
ieee 1 evel 
Bos = ice co), 194) 0) 
Thee aoe = jl a3] 
2.64 1.20 = Mb 
so that 


Le? = .00354. 


A comparison of the corresponding approximate intervals gives: 


Interval Used 


Lower 5% limit 


Upper 5% limit 


iclhiicia ee = ep oe, bee ee .0020 .0035 
INUGy Maer Wal ee cee Makes nen Gtomch 6-"c .0019 .0035 
Chi-square (Satterthwaite) .... . .0021 .0037 


— 


so that, as would be expected due to the large number of degrees of 
freedom, the three methods lead to similar results. ; 

In small samples the methods do not always lead to similar results 
as the following data on Plankton hauls (6) indicates. The analysis of 
variance of vertical hauls (Table III) 


Source of Variation df. M.S. E(M.8.) 


Eland 2 8) i, ee ee A ee 6 1011 CGH + 6a, 


He eee Ns x say 35 .0208 oH? 
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leads to an estimate of oy of .0134. The interval estimates are given 
below. The Confidence method will be discussed briefly in the next 
section. 


Interval Used Lower 5% limit} Upper 5% limit 
rca. @ ust oe See el tS .00420 .0580 
GonhiGences: Aen er us. ses ee kere, 8 .00347 .0604 
(CLC SGTIN GED Ree San ae near ee Meee es .00556 .0812 
|GET .* Os Arne aan ee A a ee .00000 0294 


The lower limit of the normal approximation is inconsistent with 
the F-test which is significant at the 5% point. 
The analysis of variance of Haul 327 


Source of Variation d.f. M.S. E(M.8.) 
ET S11 PE DE sone eral ot voy eos tor nd 59 9 . 1926 oGH + box? 
(Cavayay 24 1S hr een ee ems ee fie 45 .0970 oon? 


gives .0159 as an estimate of oj. The intervals are 


Interval Used Lower 5% limit | Upper 5% limit 
irceiver alliemae eu Meee keh ec ey as tial See .00000 .0690 
Contdencesmew wd ee af iste Sf ee .00000 .0740 
Gig Sep aLe mar pees riers) sie es oe < o .00532 . 3090 
INGEN! A ORs sco, ees Oe eee .00000 .0414 


In this case F is not quite significant at the 5% level so the Chi 
Square result is spurious. A non-zero lower limit will occur with the 
Fiducial method if and only if the F-test is significant so that contra- 
dictions of this nature are impossible. Moreover, as has been indicated 
earlier, it may be used even when the degrees of freedom are small. 


AN APPROXIMATE CONFIDENCE INTERVAL 


A rough confidence interval may readily be obtained by using the 
fact that when o? = 0, F is distributed as ordinary F times o;/o, . Hence 


144 BIOMETRICS, JUNE 1950 


Ptr Si, = a 
02 


By manipulating the quantities within the parenthesis, it follows at 
once that an exact lower confidence limit is 


ie )2 

(- : a 

By substituting a sample value for the parameter, a rough confidence 
interval may be obtained, 


F 
Ry we 2 
F =, 1 Op - 
This result greatly resembles the fiducial approximation. For the first 
numerical example previously considered, the confidence interval ap- 
proximation leads to the interval .0018—.0039. It would appear from 
considerations of expected values (4) that this rough confidence interval 
tends to be over-conservative. It is interesting to note that, despite 
theoretical divergences, confidence and fiducial interval methods tend 
to lead to approximately the same results here. 
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DETERMINING SCALES AND THE USE OF 
TRANSFORMATIONS IN STUDIES ON 
WEIGHT PER LOCULE OF TOMATO FRUIT 


LeRoy Powers 


Principal Geneticist, Bureau of Plant Industry, Soils, and Agricultural Engineering, 
Agricultural Research Administration, United States Department of Agriculture 


ia BIOLOGICAL RESEARCH the scale on which the variates have been 
measured may not be in harmony with the nature of the action of the 
biological processes bringing about the expression of the character under 
study, as was pointed out by Wright in 1926 and Fisher, Immer, and 
Tedin in 1932. In such an event the data will not follow the normal 
probability integral and, if statistical methods of analyzing the data 
assuming normal distribution are to be used, some transformation of 
the data causing it to follow the normal probability integral is essential. 
Mather 1949 suggests that the use of special scales employing some metric 
or metrics that would be in harmony with the nature of the forces 
bringing about the expression of the character might furnish a solution. 
The use of a conventional scale, and if necessary some transformation, 
is preferred from the standpoint that the equipment for making the 
measurements in the case of the conventional scale is more readily 
available and usually does not require any special training for its use. 
Then the problem resolves itself into determining the scale that is in 
harmony with the nature of the processes bringing about the expression 
of the character undergoing investigation and if necessary employing 
a transformation that will cause the data to follow the normal probability 
integral. 

This paper is concerned with determining scales and the use of trans- 
formations in studies on weight per locule of tomato fruit. The locule of 
the tomato fruit is the chamber containing the seed, its juices, and the 
placental tissue. Weight per locule was calculated by dividing average 
weight per fruit of each plant by average number of locules. Hence, 
the weight per locule also includes other tissues of the tomato fruit, such 
as the fleshy part which is composed of the locular walls. The data are 
extensive, involving a large number of plants of four hybrids and covering 
a period of three years. 
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In the present study the genetic design of the experiment included 
Pe, By to P, fy. > 7 Bato Per, ands, LoInace populations and the 
field design of the experiment was a randomized complete block with 20 
replications. The experiments conducted in 1938 had 24 plants planted 
per plot, and two plots of each of the segregating populations (B, to P, , 
F, , and B, to P,) were grown per replication to a single plot of any non- 
segretating population (parents and F,). In 1939 and 1940 only 12 
plants were planted per plot, and one plot of each population was grown 
per replication. 


EXPERIMENTAL RESULTS 


In analyzing data to determine the appropriate scale or scales all 
constants must be calculated from the individual plant data, and not 
from frequency distributions. Likewise, transformations should be ap- 
plied to the individual plant data. Another precaution involves the use 
of chi-square for testing goodness of fit between the obtained and theo- 
retical frequency distributions. In case the theoretical frequency 
distribution is calculated from the mean and standard error of the 
individual plant data from which the obtained frequency distribution 
is constructed, then chi-square calculations must be based on the differ- 
ences between the theoretical and obtained frequency distributions and 
not upon a common frequency distribution such as is used when the 
theoretical frequency distribution is calculated from a provisional genetic 
hypothesis based on the data from the nonsegregating populations. 


Scale Primarily Logarithmic 


The data from the Danmark (Lycopersicon esculentum Mill.) X Red 
Currant (L. pimpinellifoliwm Mill.) tomato hybrid furnish an example of 
a scale primarily logarithmic. The first step in determining the scale 
capable of describing the nature of the action of the processes responsible 
for variation of a character is to compare the obtained means with those 
calculated on the basis of an arithmetic progression and those calculated 
on the basis of a geometric progression. 

The obtained means and standard errors, theoretical arithmetic and 
geometric means, grand total variances, and total number of individuals 
in the study on weight per locule of tomato fruit are given in table 1. 
In every instance the theoretical arithmetic means are larger than those 
obtained, whereas, the theoretical geometric means are not significantly 
different from those obtained. These results indicate that the effects of 
the forces responsible for the variability of weight per locule of the fruit 
in the populations of the Danmark X Red Current hybrid are multi- 
plicative. Then, the appropriate scale is logarithmic. 
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TABLE 1 


OBTAINED MEANS AND STANDARD ERRORS, THEORETICAL ARITHMETIC AND 

GEOMETRIC MEANS, GRAND TOTAL ARITHMETIC VARIANCES, AND TOTAL NUMBER 

OF INDIVIDUALS FOR WEIGHT PER LOCULE OF FRUIT; DANMARK X RED CURRANT 
TOMATO HYBRID GROWN IN 1938* 


Mean 
Grand Number 
Population Theoretical total of 
Obtained arithmetic | individ- 
Arith- Geo- variance uals 


metic metric 


Grams Grams | Grams 
Red Currant 0.45 + 0.017 0.017588 420 
B,; to Red Currant 97+ .045 2.92 0.99 . 183857 932 
ia ley see he Ge 2.33 + .130 5.40 2.16 . 759823 475 
Iiirs BO > <i ata Ox 2.12+ .105 5.40 2.16 1.332184 932 
B, to Danmark 4.824 .253 7.88 Ane 5.213842 928 
Danmark 10.386 + .581 18. 286069 457 


*In this table and all following tables which involve the Danmark X Red Currant hybrid, D. 
signifies Danmark and R., C., Red Currant. 


As a consequence of these findings, the individual plant data for 
weight per locule of the fruit of these hybrid and parental populations 
were transformed to logarithms. The obtained means and standard 
errors, theoretical logarithmic means, and grand total logarithmic vari- 
ances of weight per locule are listed in table 2. After transformation, 
the theoretical arithmetic formulae are used in calculating the theoretical 
means of the logarithms. The reason for so doing is that multiplicative 
effects are additive on the logarithmic scale. As was to be expected the 
theoretical means are not significantly different from those obtained. 
Since the two parental means were used in calculating the theoretical 
means of the other four populations given in tables 1 and 2, these figures 
furnish only indirect evidence as to the nature of the action of the forces 
causing environmental variability of the parental plants. 

To obtain further information concerning the scale appropriate for 
describing the nature of the action of the forces causing the variability 
noted for weight per locule of tomato fruit, the data were classified into 
frequency distributions. The obtained frequency distributions for weight 
per locule for the two parental and F’, populations are shown in table 3. 
In this table the upper class limits are given for both the arithmetic and 
logarithmic scales. If the appropriate scale is logarithmic then the 
frequency distributions of table 3 should follow the normal probability 
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TABLE 2 
OBTAINED LOGARITHMIC MEANS AND STANDARD ERRORS, THEORETICAL LOGA- 
RITHMIC MEANS, AND GRAND TOTAL LOGARITHMIC VARIANCES FOR WEIGHT PER 
LOCULE OF FRUIT; DANMARK X RED CURRANT TOMATO HYBRID GROWN IN 1938 


Mean Grand total 

Population logarithmic 
Obtained Theoretical variance 

logarithmic logarithmic 

Red Currant —0.364833 + 0.018357 0.018692 
B, to Red Currant —0.051210 + 0.014673 —0.029018 | 0.033374 
18) 9 ID), S€ IR, C. 0.334631 + 0.026734 0.306296 0.031875 
15. 1D, S€ Uke (Ce 0.272647 + 0.014645 0.306296 | 0.045940 
B, to Danmark 0.635670 + 0.017059 0.641611 0.043278 
Danmark 0.976926 + 0.026607 0.035156 


integral. The theoretical frequency distributions based on the assump- 
tion that the data do follow the normal probability integral can be 
calculated from the obtained means and the standard errors of a single 
determination, which in turn are calculated from the grand total 
variances. 

The grand total variances for the populations must be used in 
estimating the standard errors of a single determination to be employed 
in calculating theoretical frequency distributions, because the obtained 
frequency distributions include the variability due to replications as 
well as all other variability within a population. The arithmetic grand 
total variances are given in table 1 and the logarithmic in table 2. 
The details of calculating the theoretical frequency distributions will 
be given later in connection with the same calculations for the segregat- 
ing populations. Chi-square values were calculated to determine whether 
the deviations noted between the obtained frequency distributions and 
the theoretical frequency distributions were greater than expected on 
the basis of random sampling. 

The degrees of freedom, chi-square values, and P values for testing 
goodness of fit between the obtained frequency distributions and those 
calculated on the basis of the arithmetic and logarithmic scales are given 
in table 4. Since all of the values for weight per locule of the fruit for 
the Red Currant parent fell in the first class of the frequency distribution, 
there are no values for this parent. All of the chi-square values for the 
theoretical frequency distributions calculated on the arithmetic scale 
are too large to be attributable to probable errors of random sampling. 
The P value of the chi-square calculated from the F, obtained and 
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TABLE 4 


DEGREES OF FREEDOM, CHI-SQUARE, AND P VALUES FOR TESTING GOODNESS OF 

FIT BETWEEN THE OBTAINED FREQUENCY DISTRIBUTIONS AND THOSE CALCU- 

LATED ON THE BASES OF THE ARITHMETIC AND LOGARITHMIC SCALES FOR 

WEIGHT PER LOCULE OF FRUIT FOR THE Pi, F1 AND P: POPULATIONS; DANMARK X 
RED CURRANT TOMATO HYBRID GROWN IN 1938 


Degrees of freedom Chi-square P lies between 
Population 
Arith- | Loga- | Arith- | Loga- | Arithmetic | Logarithmic 
metic | rithmic | metic | rithmic 
Red Currant 0 0 — —_ — — 
1a 5 ID, OS TR OL 3 4 11.792 | 14.142 | 0.01 & — | 0.01 & — 
Danmark 16 16 48.393 | 29.546 | 0.01 & — | 0.05 & 0.02 


logarithmic frequency distributions is less than 0.01, and P for the 
corresponding calculations for the Danmark population is less than 0.05 
but larger than 0.02. Another population of Danmark was grown in this 
same experiment, but in connection with the Danmark X Johannisfeuer 
cross. Chi-square calculated from the obtained frequency distribution 
and the theoretical logarithmic frequency distribution had a P value 
lying between 0.50 and 0.30, showing that the proper scale must be 
logarithmic for the data from this parent. All of the chi-square values 
listed in table 4 can be calculated from the data given in tables 1, 2, and 3. 

If the genetic variability is due to multiplicative effects of the forces 


TABLE 5 
MEANS, GRAND TOTAL VARIANCES, STANDARD ERRORS OF A SINGLE DETERMINA 
OF FRUIT; SEGREGATING POPULATIONS OF THE DANMARK X RED CURRANT 
TRANSFORMED TO LOGARITHMS 


Frequency Distribution 


Population Mean Variance | Standard Upper limit of class 
error 

— ) a) nN oD 
ron = oS rn) 
=) re ro) nN a) 
a) & x 69 ro) 
i D> SY 19 = 
i oe) oe S x 
° —) S ° ( 
% % % % % 

Bito Red Currant .. . —0.051210 | 0.033374 | 0.182686 | 89.3 10.0 0.6 0.1 

Uppy ND) GUE CR i Gm 6 e 0.272647 | 0.045940 | 0.214336 | 32.6 39.3 17.9 (hae 22583 


‘Bio Danmi enka eens 0.635670 | 0.043278 | 0.208034 14" 113.2053 82052 5-9 
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causing it, as was indicated by a comparison between the obtained and 
theoretical means given in tables 1 and 2, and if such is true for the 
environmental variability also, then the obtained frequency distributions 
of the segregating populations should not deviate from their correspond- 
ing theoretical logarithmic frequency distributions further than can be 
explained by chance. The means, standard errors of a single determina- 
tion, and theoretical frequency distributions for weight per locule of 
fruit for segregating populations of the Danmark X Red Current hybrid 
are given in table 5. The theoretical frequency distributions given in 
table 5 were calculated from the means and standard errors of a single 
determination by use of Pearson’s (1930) tables of the normal proba- 
bility integral, namely, that portion of table II giving the area and 
ordinate in terms of the abscissa. In using this table of Pearson’s 
the value under the column heading “x” is obtained by subtracting 
any given class from the mean and dividing the remainder by the 
standard error. As an example, take the B, to Danmark population 
given in table 5 and the class having the column heading 0.397940. 
We have (0.634670 — 0.397940) divided by 0.208034, which results in 
an “x” value of 1.14. Looking this value up under column heading ‘‘a”’ 
of Pearson’s (1930) table II, page 3, we find the value 1/2(1 + a) to be 
0.873. Since the value of ‘‘x”’ is positive this must be subtracted from 
1.0 and multiplied by 100 to give the percent of the population having 
a logarithm equal to or less than 0.397940. The value obtained is 12.7 
percent. Since 1.4 percent of the population is expected to fall in the 
preceding class, the proportion of the population expected in the 
0.397940 class is 12.7 percent minus 1.4 percent, or 11.3 percent as 


TABLE 5—Continued 


TION, AND THEORETICAL FREQUENCY DISTRIBUTIONS FOR WEIGHT PER LOCULE 
TOMATO HYBRID GROWN IN 1938; INDIVIDUAL PLANT DATA FOR WEIGHT IN GRAMS 


Frequency Distribution 
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shown in table 5. For asecond example, take the class headed 0.977724. 
We have (0.635670 — 0.977724) divided by 0.208034, which gives an 
‘“o” value of —1.64. The value of 1/2(1 + a) for the “x” value of 1.64 
is 0.949 (Pearson’s table IT). Since the ‘x’ value is negative this value 
is not subtracted from 1.0 but is multiplied by 100 to give 94.9 percent, 
which is the percent of the population having a logarithm less than or 
equal to 0.977724. Since 92.1 percent of the population is expected to 
fall in classes having a lower value, 94.9— 92.1, or 2.8 percent of the 
population is expected to fall in the 0.977724 class. The values for all 
of the other theoretical frequency distribution tables, whether on a 
logarithmic or arithmetic scale, were calculated in a manner identical 
to the above given examples. 

The theoretical frequency distributions given in table 6 were calcu- 
lated by multiplying the percent expected in any class by the total 
number of individuals in the population. Thus the theoretical number 
for the 0.176091 class of the B, to Danmark population is 1.4 percent of 
928, or 13 individuals. With the exception of the frequency distribu- 
tion used in calculating the chi-square value of the B, to Red Currant 
population, in which classes 1 and 2 were combined, the classes were 
grouped so that any given theoretical frequency distribution did not 
have less than 10 individuals in any class. This grouping was started 
at the extremes of the frequency distributions as indicated in table 6. 
In case of the exception noted (B, to Red Currant population) the 
0.544068 and 0.653212 classes were grouped, even though after so doing 
the theoretical frequency distribution had only 7 individuals in the 
last class. 

An examination of table 6 reveals that when classes 1 and 2 are not 
combined the deviations between the obtained and theoretical frequency 
distributions for the B, to Red Currant population are greater than ex- 
pected, due to the probable errors of random sampling. This is shown 
by the high chi-square value (10.765) and its corresponding P value 
which is less than 0.01. The corresponding chi-square value for the B, 
to Danmark population is also somewhat large. In fact, for every 
population of table 6 the 0.176091 class (first class) had more individuals 
in the obtained frequency distributions than in the theoretical. The 
reverse is true for the 0.397940 class (second class). It seems that the 
data transformed to logarithms are not normally distributed for classes 
1 (0.176091 class) and 2 (0.397940 class). The chi-square values with 
their corresponding P values show that when the first and second classes 
are combined the deviations between the obtained and theoretical fre- 
quency distributions are no greater than expected due to chance fluctua- 
tions, as all of the P values lie between 0.30 and 0.20. 
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From the above analysis the following conclusions can be drawn. 
The data transformed to logarithms are not normally distributed for 
the first two classes of the frequency distribution, but are normally 
distributed after the first two classes have been combined. ‘This fact 
could not be detected by a study of the means alone. Hence, in determin- 
ing the scale, or scales, capable of describing the nature of the forces 
bringing about variation of a character, it is desirable to include in the 
analysis both the means and the frequency distributions. Primarily 
the data follow the logarithmic scale, and in making a genetic analysis 
of the data by the partitioning method (Powers, Locke, and Garrett, 
1950) the weights per locule of individual plants should be trans- 
formed to logarithms and the first two classes of the frequency distribu- 
tions combined. 


Scale primarily arithmetic 


The obtained means and standard errors, theoretical arithmetic and 
geometric means, and total number of individuals for weight per locule 
of fruit for the Johannisfeuer X Bonny Best tomato hybrid grown in 1939 
are given in table 7. None of the arithmetic means differ significantly 
from their comparable obtained means. However, without exception 
the theoretical arithmetic means are of smaller magnitude than their 
corresponding obtained means. Also, the geometric means are con- 
siderably smaller than their comparable obtained means, and for most 
populations significantly so. Clearly, the obtained means do not follow 
a geometric progression. Likewise, though not as decisive, the evidence 
indicates that the obtained means do not follow an arithmetic progression 
either. This does not rule out either the arithmetic or logarithmic scale, 
as partial dominance of greater weight per locule of the tomato fruit 
could be responsible for the rather poor agreements noted. It will be 
remembered that the formulae for calculating both the theoretical 
arithmetic and geometric means assume no dominance. 

In order to obtain further information as to the nature of the action 
of the biological processes bringing about the expression of the character 
weight per locule of fruit in the tomato hybrid populations of the cross 
Johannisfeuer X Bonny Best, an analysis was made of the frequency 
distributions. The degrees of freedom, and chi-square and P values for 
testing goodness of fit between the obtained frequency distributions and 
those calculated on the bases of the arithmetic and logarithmic scales 
are listed in table 8. With the exception of the B, to Johannisfeuer popu- 
lation the fit is good between the theoretical arithmetic frequency 
distributions and their corresponding obtained frequency distributions. 
P is considerably less than 0.01 for the chi-square value calculated from 
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TABLE 7 


LOCULE OF FRUIT; JOHANNISFEUER X BONNY BEST TOMATO HYBRID 
GROWN IN 1939* 


Mean 
Population Theoretical Number of 
Obtained individuals 
Arithmetic | Geometric 
Grams Grams Grams 
Johannisfeuer 6.20 + 0.272 224 
B, to Johannisfeuer 9.96 + ~.510 9.44 8.22 224 
1 sated he Ml BAM BS 13.31 + .469 12.68 10.90 227 
Eee se els eB: 13.01 + .546 12.68 10.90 222 
B, to Bonny Best 17.02 + .563 15.92 14,45 2211 
Bonny Best 19.15 + .487 213 


*In this table and all following tables wnich involve the Johannisfeuer * Bonny Best, hybrid J. 
signifies Johannisfeuer and B. B., Bonny Best. 


the theoretical arithmetic frequency distribution and obtained frequency 
distribution of the B, to Johannisfeuer population. The fit between the 
theoretical logarithmic frequency distribution and the obtained fre- 
quency distribution for this population is good, as the P value for the 
chi-square test lies between 0.70 and 0.50. With the exception of the 
Bonny Best and B, to Johannisfeuer populations the fits between the 
theoretical arithmetic frequency distributions and the comparable ob- 
tained frequency distributions are better than the fits between the 
theoretical logarithmic frequency distributions and their comparable 
obtained frequency distributions. In case of the Bonny Best population 
both chi-square values have a P value that lies between 0.20 and 0.10. 

The scale is primarily arithmetic, but it should be pointed out that 
the data are not as discriminatory as is desired, and more work on some 
populations is necessary to discriminate between scales. However, in 
making a genetical analysis of these data by the partitioning method the 
arithmetic scale is satisfactory for all populations except the B, to 
Johannisfeuer. To avoid metrical bias the weights per locule for the 
individual plants of this population must be transferred to logarithms. 
Just why the data for the B, to Johannisfeuer follow a logarithmic scale 
is not clear, especially since the fF, and Johannisfeuer populations give 
such a good fit to the arithmetic scale. The fact that the obtained fre- 
quency distributions are normal when the data are expressed on one or 
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TABLE 8 
DEGREES OF FREEDOM, CHI-SQUARE, AND P VALUES FOR TESTING GOODNESS OF 
FIT BETWEEN THE OBTAINED FREQUENCY DISTRIBUTIONS AND THOSE CALCU- 
LATED ON THE BASES OF THE ARITHMETIC AND LOGARITHMIC SCALES FOR 
WEIGHT PER LOCULE OF FRUIT; JOHANNISFEUER X BONNY BEST TOMATO HYBRID 
GROWN IN 1939 


Degrees of Chi- P lies between 
freedom square 
Population 
Arith-| Loga-| Arith-| Loga-| Arithmetic | Logarithmic 
metic | rith- | metic | rith- 
mic mic 
Johannisfeuer 6 7 | 6.010/32.484| 0.50 & 0.30 | 0.01 & — 
B, to Johannisfeuer 12 il BBE WHS) SB 7@ai| IL ais — .70 & 0.50 
Hi 5 do OS dats 18}. 11 11 | 8.956)12.679} .70& .50 .50 & .30 
1ilis5 do. DS 1835 18% 14 14 |11.723/80.456} .70& .50 OL & — 
B, to Bonny Best 14 14 |13.084]19.540} .70& .50 .20& .10 
Bonny Best 14 12 |19.349)17.245) .20& .10 .20 & .10 


the other, the arithmetic or the logarithmic scales, together with the fact 
that the obtained means are larger than either the theoretical arithmetic 
or logarithmic means proves phenotypic dominance of greater weight per 
locule of fruit in the Johannisfeuer X Bonny Best tomato hybrid, and 
also is convincing evidence in support of genetic dominance. 


Scales arithmetic and logarithmic 


The obtained means and standard errors, theoretical arithmetic and 
geometric means, and total number of individuals for weight per locule 
of fruit for Johannisfeuer X Red Currant tomato hybrid grown in 1939 
are given in table 9. In every instance the arithmetic means are larger 
than those obtained and the geometric smaller. By employing the ob- 
tained mean of the Red Currant parent and the obtained mean of the F, 
population of the Johannisfeuer X Red Currant hybrid, the theoretical 
geometric mean is 1.09 grams per locule, whereas the obtained mean is 
1.04 + 0.040 grams. Similarly, by employing the mean of the F, 
population of the Johannisfeuer X Red Currant hybrid and the mean of 
the Johannisfeuer parent, the theoretical arithmetic mean of the B, 
to Johannisfeuer population is 4.45 grams per locule, whereas the ob- 
tained mean is 4.48 + 0.139 grams. From these results one might expect 
the appropriate scale for the Red Currant and B, to Red Currant popu- 
lations to be logarithmic, and that for the B, to Johannisfeuer and 
Johannisfeuer populations to be arithmetic. 
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TABLE 9 
OBTAINED MEANS AND STANDARD ERRORS, THEORETICAL ARITHMETIC AND 
GEOMETRIC MEANS, AND TOTAL NUMBER OF INDIVIDUALS FOR WEIGHT PER 
LOCULE OF FRUIT; JOHANNISFEUER X RED CURRANT TOMATO HYBRID GROWN 


IN 1939* 
Mean 
Number of 
Population Theoretical individuals 
Obtained 
Arithmetic | Geometric 
Grams Grams Grams 
Red Currant 0.44 + 0.015 229 
B, to Red Currant 1.04 + .040 1.88 0.85 233 
Hingetelee <a we (Ce 2.70 + .058 3.02 1.65 225 
lite ds Seed Gs 2.12+ .106 3.32 1.65 225 
B, to Johannisfeuer 4.48 + .1389 4.76 3.20 230 
Johannisfeuer 6.20 2- .272 224 


*In this table and all following tables which involve the Johannisfeuer X Red Currant hybrid, 
R. C. signifies Red Currant and J., Johannisfeuer. 


The degrees of freedom, chi-square, and P values for testing goodness 
of fit between the obtained frequency distributions and those calculated 
on the bases of the arithmetic and logarithmic scales for weight per 
locule of fruit for the Johannisfeuer * Red Current hybrid grown 
in 1939 are given in table 10. Since all of the Red Currant plants fell into 
one class, this population does not provide any information as to the scale 
appropriate for describing the nature of the action of the processes 
bringing about the variability of the character weight per locule. For 
the B, to Red Currant population the chi-square values are those ex- 
pected, if the appropriate scale is arithmetic and not logarithmic, as 
was indicated by a study of the obtained and theoretical means. The 
chi-square values also indicate that the proper scale of measurement for 
the F, generation is arithmetic. The same is true of the Johannisfeuer 
population, but the B, to Johannisfeuer population data do not seem to 
be following the arithmetic scale, as was indicated by the analysis of 
the means. The chi-square values show that the data are following 
the logarithmic scale. The F, population data do not seem to be 
following either scale, as was to be expected since the B, to Red Currant 
followed the arithmetic scale and the B, to Johannisfeuer the logarithmic 
scale. 

This behavior of the B, to Red Currant and B, to Johannisfeuer 
populations, the former following the arithmetic scale and the latter 
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the logarithmic scale, needs to be considered further, as the reverse was 
expected from the analysis of the means. These findings raise the 
the question whether in case of the B, to Red Currant population the 
nature of the action of the processes causing the genetic variability is 
most adequately described by the logarithmic scale and the nature of 
the action of the processes causing the environmental variability is most 
adequately described by the arithmetic scale. If such were the case, 
the opposite would be true of the B, to Johannisfeuer population, that 
is, the genetic variability would be following the arithmetic scale and 
the environmental variability the logarithmic scale. Then, since both 
genetic and environmental variability would be present in these two 
populations the obtained frequency distribution would be a combination 
of the two and therefore might be expected not to give a good fit when 
tested against either the theoretical arithmetic or theoretical logarithmic 
frequency distributions. However, since the environmental variability 
makes up the greater proportion of the variability, as regards weight 
per locule, and the genetic variability a rather small proportion, the fit 
between the obtained frequency distribution and the theoretical arithme- 
tic frequency distribution might be good in case of the D, to Red Currant 
population and the fit between the obtained and the theoretical loga- 
rithmic frequency distribution good in case of the B, to Johannisfeuer 
population. This would mean that the data are not discriminatory as 
regards the genetic variability. In other words the environmental 
variability forms such a large proportion of the total variability as to 
obscure the scale that is followed by the genetic variability. 

The obtained means and standard errors, theoretical arithmetic and 
geometric means, and total number of individuals for weight per locule of 
fruit for Danmark X Johannisfeuer tomato hybrid grown in 1938, 1939, 
and 1940 are listed in table 11. For the 1938 data, with the exception of 
the B, to Danmark population, the obtained means are larger than either 
the theoretical arithmetic means or the theoretical geometric means. 
In case of the exception noted, the theoretical arithmetic mean is larger 
than the obtained mean, but not significantly so. For the 1939 data 
there are no exceptions, as in every case the obtained means are larger 
than the theoretical means, whether arithmetic or logarithmic. For 
the 1940 data, the differences between the theoretical arithmetic means 
and their respective obtained means for the B, to Johannisfeuer and F, 
populations are not significant, but for the 7, and B, to Danmark 
populations the theoretical arithmetic means are significantly larger than 
the obtained means. For the same year the logarithmic means are 
not materially different from the obtained means for the F, and B, 
to Danmark populations. However, the logarithmic means are sig- 
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TABLE 10 
DEGREES OF FREEDOM, CHI-SQUARE, AND P VALUES FOR TESTING GOODNESS OF 
FIT BETWEEN THE OBTAINED FREQUENCY DISTRIBUTIONS AND THOSE CALCU- 
LATED ON THE BASES OF THE ARITHMETIC AND LOGARITHMIC SCALES FOR 
WEIGHT PER LOCULE OF FRUIT; JOHANNISFEUER X RED CURRANT TOMATO 
HYBRID GROWN IN 1939 


Degrees of Chi- P lies between 
Population freedom square 
Arith-| Loga- | Arith-| Loga- 
metic | rith- | metic | rith- | Arithmetic | Logarithmic 
mic mic 

Red Currant 0 0 — — — — 
B, to Red Currant 1 1 0.193} 5.819} 0.70 & 0.50 | 0.02 & 0.01 
(Oh WA 4 awe OF z 2 1.665] 9.529) .50& .30 01 & — 
ry Jape bai C 3 4 | 8.692]17.856) .05& .02 01 & — 
B, to Johannisfeuer 5 6 |29.092) 7.067; .01 & — .50 & .30 
Johannisfeuer 6 fi 6.010/382.484; .50& .30 01 & — 


nificantly smaller than the obtained means for the B, to Johannisfeuer 
and F,, populations. Then, the analysis of the means indicates that both 
arithmetic and logarithmic scales are necessary to describe the nature of 
the action of the genetic processes differentiating the six populations of 
the Danmark X Johannisfeuer hybrid for the 3 years of the tests. 

The degrees of freedom, chi-square, and P values for testing goodness 
of fit between the obtained frequency distributions and those calculated 
on the bases of the arithmetic and logarithmic scales for weight per 
locule of fruit of the Danmark * Johannisfeuer tomato hybrid grown in 
1938, 1939, and 1940 are listed in table 12. For the 1938 data, the 
theoretical arithmetic frequency distribution and the obtained frequency 
distribution are in agreement for the Johannisfeuer population and the 
logarithmic frequency distribution is rejected. For the F, and Danmark 
populations the fit is good between the obtained and theoretical loga- 
rithmic frequency distributions, and the theoretical arithmetic frequency 
distributions are rejected. The obtained frequency distributions of the 
segregating populations (B, to Johannisfeuer, F, , and B, to Danmark) 
do not, in any case, give a good fit to either the theoretical arithmetic 
or theoretical logarithmic frequency distributions. In all cases the data 
are discriminatory. The results are those expected on the basis that 
the processes bringing about the variation of the character, weight per 
locule of fruit, are most adequately described for some genotypes by the 
arithmetic scale and for other genotypes by the logarithmic scale. In 
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TABLE 11 
OBTAINED MEANS AND STANDARD ERRORS, THEORETICAL ARITHMETIC AND 
GEOMETRIC MEANS, AND TOTAL NUMBER OF INDIVIDUALS FOR WEIGHT PER 
LOCULE OF FRUIT; DANMARK X JOHANNISFEUER TOMATO HYBRID GROWN IN 
1938, 1939, AND 1940* 


Mean 
Number of 
Year and Population Theoretical individuals 
Obtained 
Arithmetic | Geometric 
Grams Grams Grams 
1938 
Johannisfeuer 4.61 + 0.446 452 
B, to Johannisfeuer 6.72 + .425 5.94 5.58 928 
By Ds SK ae 7.96 + .419 (e206 6.76 469 
fi 5 ID. SK de 8.35 + .467 e206) 6.76 932 
B, to Danmark 8.32 + .399 8.59 8.19 921 
Danmark 9592722) 091 456 
1939 
Johannisfeuer 6.20 + .272 224 
B, to Johannisfeuer 9.45 + .426 8.28 AOS 230 
i, ID, SK dh 11.81 4+ .516 10.36 9.48 209 
fs 5 IDS S< di 10.92 + .382 10.36 9.48 215 
B, to Danmark 13.74 + .476 12.44 11.73 231 
Danmark 14.51 + .3860 228 
1940 
Johannisfeuer 7.34+ .062 220 
B, to Johannisfeuer 10.144 .146 9.80 9.08 224 
(9h IDs S€ dle 11.42 + .148 117), Paz/ 11.24 224 
185, IDs SK dh. 12.41 + .080 IPAS Til AL 223 
B, to Danmark 13.387 + .253 14.74 13.90 219 
Danmark 17.20 + .259 219 


*In this table and all following tables which involve the Danmark X Johannisfeuer hybrid, D. 
signifies Danmark and J., Johannisfeuer, 


such an event the frequency distributions of the nonsegregating genera- 
tions would be expected to show a good fit to a frequency distribution 
based on one or the other scale, whereas the obtained frequency distribu- 
tions of the segregating generations would not be expected to give a good 
fit when tested to either one or the other scale. As has been shown such 
was the situation. 

For 1939 the obtained frequency distributions of the Johannisfeuer 
and F’, populations gave a good fit when tested for the arithmetic scale 
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TABLE 12 
DEGREES OF FREEDOM, CHI-SQUARE, AND P VALUES FOR TESTING GOODNESS OF 
FIT BETWEEN THE OBTAINED FREQUENCY DISTRIBUTIONS AND THOSE CALCU- 
LATED ON THE BASES OF THE ARITHMETIC AND LOGARITHMIC SCALES FOR 
WEIGHT PER LOCULE OF FRUIT; DANMARK X JOHANNISFEUER TOMATO HYBRID 
GROWN IN 1938, 1939, AND 1940 


Degrees of Chi- P lies between 
freedom square 
Year and population 
Arith-| Loga-| Arith- | Loga- | Arithmetic | Logarith- 
metic | rith- | metic | rithmic mic 
mic 
1938 
Johannisfeuer 8 10 11.636) 44.326/0.20 & 0.10/0.01 & — 
B, to Johannisfeuer 13 15 |118.466| 35.936} .01 & — 01 & — 
Ii OE ee 14 13 67.582} 6.800) .01 & — .95 & 0.90 
lie. NDE ee Al 16 19 |186.745} 33.310} .01 & — .02 & .01 
B, to Danmark 16 17 77.482| 47.078} .01 & — .01 & — 
Danmark 16 16 52.572) 17.016) .01 & — .60 & .30 
1939 
Johannisfeuer 6 7 6.010) 32.484; .50 & .30) .01 & — 
B, to Johannisfeuer (a! fil 8.738] 6.244 .70 & .50} .90 & .80 
1D ea ke 10 ila 4.306} 27.528} .95 & .90| .01 & — 
Ji IDE ye AL 10 10 N70) SBR) AHO es QAO) nO es  Q4i0) 
B, to Danmark 10 12 28.285) 17.609) .01 & — Pies Ke) 
Danmark it 10 12.183] 6.056) .50 & .30/ .90& .80 
1940 
Johannisfeuer fi 7 7.931} 1.801) .50 & .30) .98 & .95 
B, to Johannisfeuer ies 11 55.484] 18.524) .01 & — LONG OD 
ED) ee 13 11 40.558} 9.951) .01 & — .70 & .50 
lites dD &< de 13 11 30.241) 6.013} .01 & — .90 & .80 
B, to Danmark 13 14 22.966] 12.876] .05 & .02) .70 & .50 
Danmark 15 14 | 26.493] 14.717) .05 & .02) .50& .30 


but not when tested for the logarithmic scale. The data for the Danmark 
population were not discriminatory as the obtained frequency distribu- 
tion gave a good fit when tested for either the arithmetic or logarithmic 
scale. The same situation was found for the B, to Johannisfeuer and F, 
populations. The data for the B, to Danmark population rejected the 
arithmetic scale but the fit between the obtained and theoretical fre- 
quency distributions was fair when the test was for the logarithmic 
scale. These results support the contention that for some genotypes 
the environmental variability is most adequately described by the arith- 
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metic scale and for other genotypes by the logarithmic scale. Since in 
some cases the data for 1939 are not discriminatory, a greater number 
of individuals per population are needed. 

For 1940 the obtained frequency distributions in all cases give a good 
fit to the theoretical logarithmic frequency distributions. However, the 
data for the Johannisfeuer population are not discriminatory, as the 
obtained frequency distribution for this population also gives a good fit 
to the theoretical arithmetic frequency distribution. Then, for 1940 
the scale capable of describing the action of the biological processes 
causing variability in weight per locule of the tomato fruit for the cross 
Danmark x Johannisfeuer is primarily, if not entirely, logarithmic. 

These data from the Danmark x Johannisfeuer cross are conclusive 
in proving that for some genotypes the appropriate scale for describing 
the nature of the action of the environmental forces causing variability 
in weight per locule of fruit is arithmetic and for other genotypes is 
logarithmic. Also, evidence in support of the contention that the scale is 
not the same for all years is convincing. In 1938 and 1940 the appropri- 
ate scale for plants of the /, genotype was logarithmic and for 1939 was 
arithmetic, and in all years the data for the Ff’, were discriminatory. 


CONCLUSIONS AND SUMMARY 


1. With the exception of the Danmark (Lycopersicon esculentum) X 
Red Currant (L. pimpinellifolium) tomato hybrid grown in 1938 the 
arithmetic and logarithmic scales have been sufficient to describe the 
nature of the action of the biological processes differentiating weight 
per locule of tomato fruit. 

2. In populations of Danmark < Red Currant hybrid, the environ- 
mental variability and genetic variability of all genotypes follow the 
same scale which is logarithmic. 

3. This is not true for the populations of any other hybrid, as the 
environmental variability of some genotypes was found to be arithmetic 
and that of others logarithmic. 

4. The same was true of the genetic variability for the Johannisfeuer 
(L. esculentum) X Red Currant hybrid, as the mean of the B, to Red Cur- 
rant was geometric and that of the B, to Johannisfeuer was arithmetic. 

5. When using the partitioning method of analyzing segregating 
populations, the environmental scale must be employed as the data for 
any given segregating population is partitioned on the basis of genotypes. 
By so doing the genetic variability is removed. 

6. In those cases where the environmental variability of all genotypes 
of a hybrid are not following the same scale, the transformations Gf any 
are necessary) must be based on the scales indicated by the environ- 


DETERMINING SCALES IN TOMATO FRUIT STUDIES 163 


mental variability of each individual genotype. 

7. For the populations of any given hybrid the genetic variability 
may be following one scale and the environmental variability another. 
This was found to be true for some genotypes of the Johannisfeuer X 
Red Currant and the Danmark X Johannisfeuer populations. 

8. Likewise, the scale is not necessarily the same for all years, as in 
1938 and 1940 the environmental variability of the F,; population of 
Danmark X Johannisfeuer hybrid followed the logarithmic scale and 
in 1939 the arithmetic scale. 

9. Red Currant, Danmark, and Johannisfeuer were crossed in every 
possible combination. For the Danmark x Red Currant hybrid the 
environmental and genetic variabilities were found to follow the loga- 
rithmic scale, primarily. Since Johannisfeuer was found to follow the 
arithmetic scale as regards environmental variability, for the Johannis- 
feuer X Red Currant and Danmark x Johannisfeuer populations the 
environmental variability of some genotypes would be expected to 
follow the arithmetic scale and that of others, the logarithmic scale. 
Such was found to be the case. Also, the same was true of the genetic 
variability. 

10. Within the range of the genotypes and environments encountered, 
both the genotype and the environment are factors in determining the 
scales appropriate for describing the nature of the action of the biological 
processes causing variability of weight per locule of tomato fruit. This 
raises the question whether the variability for any given genotype within 
one replication may be following one scale and within another replication 
another scale. The data indicate that such was not the case, as all 
replications for a given genotype seemed to be following the same scale. 

11. In some cases the data were not discriminatory, indicating that 
too few individuals were grown per population. It appears that each 
population should be composed of at least 400 individuals, and popula- 
tions of 900 or more individuals are preferred. 
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QUERY: 


QUERIES 


I have been carrying out an analysis of covariance of 


77 yield (y) and age (x) in a clonal cacao experiment of split-plot 
design, the two splits being for 6 clones (C) and for Buddings v 
Cuttings (7). Table 1 gives the yields and ages of buddings and of 
cuttings over the six clones and in the six blocks. 
As you will observe, clonal yields differ quite considerably and I 
have therefore partitioned the 5 d.f. for clones (C) orthogonally into:— 


TABLE 1 
YIELD (KILOGRAMS) AND SUM OF AGES (UNIT: 10 MONTHS)* OF 8 CACAO TREES IN 
EACH OF 72 PLOTS ARRANGED IN A SPLIT-PLOT EXPERIMENT 

ti = BUDDINGS, t2 = CUTTINGS 


Clone 
Block i 44 5 70 91 95 
th te ty ty ty te t ty th to ty te 
1 W 18 29 13 12 5 3 10 4 8 Tf 16 22 
X| 78 78 ae 0 as GF (is  2EN BS Oe 78 78 
2 Ww 30022 9 12, 8 8 6 3 7 19 25 
XG) eto Cou On | a Ome WS Bil 49 64 ihe: 0G 
3 Ww Up PAR 18 12 19 5 6 16 i 2 16 18 
|) WS AS | AS} 78 64 66 «67 al Gy (Omer 
4 W ily = A 225 26 12 6 9 10 il 4 19 29 
XG 7 Sins i 7 S383 COMO GS Ce 78 TA 
5 Ve |) BS Bil 1 12 13 8 9 3 12 9 16 18 
XG |e OmeecG 78 66 78 65 lee: fil Gil OMedal 
6 Wi) Bs aA iSO) 13 a iil aly ilil a 1 1B 
XG |e Samer io WSs GS |) We Ce 71 #68 74 66 78 78 
Total Wael ey) | Bb | |) Sr | asl ie BXGy || TKO TE 11. 
465 460 | 465 394 | 464 393 | 4388 341 | 413 388 | 466 452 
Total 
for 296 185 107 108 92 225 
Clone 925 859 857 779 801 918 


*Editor’s note: Querist reported yields in grams and ages in months. The rounded numbers contain 
all the essential information. Readers who do the computations will encounter the results of rounding. 
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(i) The mean of the 3 good clones (C'S 1,44,95) v the mean of the 3 
poor clones (JCS 55,70,91) (C’ — 1 df.) 
(i) Differences within the 3 good-clone group (G — 2 df.) 
(ii) Differences within the 3 poor-clone group (P — 2 df.) 


I have adopted a similar partition in clone X block interactions 
(Error 1), in C X T and in C X T X Blocks (Error 2), with the idea of 
testing C-components for significance against their corresponding C’ com- 
ponents X block interactions. 

On carrying out the tests of significance in the covariance analysis, 
Table 2, I find C, when compared with BC, to be significant at the 0.1% 
point but when comparing C’ with C’B, or G with GB, or P with PB, 
no significant effects emerge, although I should expect C’ to be even 
more highly significant than C! When C’ is compared with BC, how- 
ever, it is found to be highly significant, yet BC and BC’ do not differ 
significantly. 

TABLE 2 
ANALYSIS OF COVARIANCE IN SPLIT PLOT EXPERIMENT 


| 
Source of Variation D.F Sx? Sry Sy? 
1B tora si (WED ee ete a eee ea 5 18.79 13.63 14.19 
Clones (GO) we Sate aes oon aS 5 186.14 Ppl, 194 340.96 
(Gaoduys oot (Gia se eine Neu 1 124.79 185.47 275.68 
WithitinGood (Gi) San cike v0 ass 2, 27.34 34.70 64.16 
Withinpledor (2)e es eens) ae 2 34.01 0.95 il 1% 
Eprroneta(s Cm a) eon serene 25 47.31 12,24. 74.39 
OLB ENS Sek ee cpaihe. bass 5 2.09 2.06 7.74 
GBI wee ome sk eS 10 9.44 4.07 54.69 
IPED. OM 6 Vn eee 10 35.78 6.11 11.96 
Budshve Cubtines (Cl) blaee v 1 141.91 6.91 0.34 
(COE 8 po, WE iaat Io OR Om onC 5 72.43 12.48 PP) ANY 
(CPI ee Go OR ee ee re 1 18.03 14.56 Hil 75 
(COR 3. ote eee ee eee 2 24.71 5.09 2.28 
JOM 2 5 anes, ee re a oe sR 2 29.69 —7.17 8.44 
DOV) 2. trl es, Co. une A eee 30 109.35 34.09 54.07 


Am I wrong in attempting to compare treatment components such 
as C’, Gand P with their interactions with blocks, or should BC be used 
in tests of significance? If the former procedure is correct, how does 
it come about that none of the C-components attains significance, espe- 
cially as no significant difference exists between BC and BC’. 

Is it possible that mere linear correction is madequate in this par- 
ticular case and that a more valid analysis of covariance could involve 


the squares of ages? 
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First, a preliminary point not inquired about: if the clones 
ANSWER: _ were selected for trial because half were good and half 

poor, then the testing of C’ is justified by the design of the 
experiment; but if you designated the clones as good and poor because 
of their performance during the experiment, then I am not happy about 
the test. See this Journal, Vol. 1: page 26; Vol. 2: page 16; and Vol. 
5: page 99. Since the trees (or the experiment) are 8 years old, and since 
I have observed no mention of good and poor clones before 1945, I 
wonder if you are justified in the testing of C’. Actually, the yield of 
Clone 44 differs more from that of Clone 1 than it does from the yields of 
the “poor” clones. 

Granting that the design of the experiment included the comparison 
of the three good clones with the three poor, I think you are correct in 
testing C’, G and P against their discrepancies with blocks. Each 
discrepance is expected to estimate the real experimental error to which 
the corresponding effect is liable. BC and BC’ are not orthogonal, so 
that test gives little information. 

The arithmetical peculiarities which puzzle you are easily explained 
by plotting the mean yields against the ages, then drawing the error 
regression line. As an example, the line connecting the points repre- 
senting the means for the good and poor clones will be found to be 
sensibly parallel to the error regression C’B, the slopes being 185/125 = 
1.5 and 2.06/2.09 = 1.0. So a large fraction of the difference in yield, 
C’, is attributable to the age difference between the plantings. This 
explains the non-significance of C’. 

I doubt if curved regression is responsible for the irregularities in 
your data. Rather, I suspect that you have gained little advantage 
from the use of age as a covariate. There are two reasons. The first 
lies in the possible confusion of ages with yield among the 8 trees of 
each plot. In some plots, the average age is little more than half that in 
others. Either all the trees in these plots have been replaced by younger, 
or else some of the replacements have scarcely come into bearing. Unless 
the yield-growth curves of all these ages are straight and identical, the 
averages do not give an accurate measure of the yield which is associated 
with any specific age. As a preliminary, I should like to have a look 
at the yields of all the trees that have a common age. This would give 
me clearer information as to whether differences among the clones 
persist among trees of given age. 

This brings me to the second reason for questioning the use of age 
as a covariate. Perhaps age is itself an informative measure of the 
value of the clones. The data seem to indicate that some clones are 
adapted to the environment and others not. In both breeding and 
selection, livability as well as yield would seem to demand consideration. 
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QUERY: [| have seen the statement offered without proof that 
78 the most economical test of the difference between two unmatched 
groups is achieved when the numbers of observations in the two 
groups are equal. Thus, let us suppose that an experimenter wants to 
test whether two different methods of learning are equally good. He 
has two groups of subjects, A and B, unmatched in any way, and has 
both groups learn the same material. Group A learns by one method; 
B by the other. Assume that the two methods really do give different 
results. In general, will this experimenter be able to discover the dif- 
ference with the smallest total number of tests if the number of subjects 
in group A equals that in group B? If this is true, is there some source 
which gives the proof for this statement? 


Suppose we are trying to measure the difference between 
ANSWER: two quantities as precisely as possible with a fixed amount 

of effort. By “precisely” we mean that we should like 
the variance of our estimate of the difference to be as small as possible. 
But, if the measurements are independent, the variance of the difference 
of the means is the sum of the variances of the two means: 


i li) NA Cs ial A 
Further, we know that 


V(z) = V(a)/n, , 


Vy) = VYy)/n, , 
where V(x), V(y) are the variances of the z and y distributions. We 
therefore want to minimize 
Viz Vi 
(@) (y) 
N, Nn, 
subject to the condition that 
Nn, +n, = Nn 
where n is fixed. It is not difficult to show that the minimization 
occurs for 
Mad iby == o,/Cy ) 
that is, that the numbers of observations should be in the ratio of the 


standard deviations. In the special case where the standard deviations 
are equal, the observations should be equally divided. 

That the distribution of effort should depend on the variability 
of each group is clear if we consider the case where one variance is 
zero, for plainly a single observation will suffice for this group. 
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In practice, of course, we usually do not know the variances before- 
hand, and so cannot use the result directly. But also, if we do not know 
anything about the variances, it seems sensible to make the numbers 
equal. 

I am afraid I cannot give you a reference to a proof, though I have 
no doubt that the result is old. Somewhat similar problems have been 


considered in connection with sample surveys. 
C. P. Winsor 


QUERY: The problem upon which I am engaged is the deter- 
79 mination of the relative effects of various factors upon yield and 
size of orange fruits. Of course, yield and size are inversely 
correlated. Our main endeavor is to determine whether or not yield 
or size is correlated with the quantity of nutrients absorbed. For this 
purpose it has seemed probable that the use of standard partial re- 
gressions would be profitable and in fact this is turning out to be the case. 
However, a question has arisen which I shall appreciate having 
answered. I have assumed that a standard partial regression coefficient 
would have a value no more than 1. However, in one of our problems in 
multiple regression we obtain a standard partial regression coefficient of 
1.3. I have looked through a number of texts and through the data and 
computations for a considerable number of partial regression problems 
which I ran some years ago, and in none of these instances have I 
happened to find a value greater than 1. I will be grateful for a statement 
as to the theoretical possibility of obtaining a value greater than 1. 


There is no arithmetical limitation on the size of a standard 
ANSWER: partial regression coefficient. My experience has coincided 

with yours in that, for the usual type of linear regression, 
these coefficients tend to be small. If I get one larger than 2, I begin 
to look for trouble; but larger values can occur. 


CORRECTION IN QUERY 74, DECEMBER, 1949. Dr. Tukey 
writes that Professor Royal F. Bloom has pointed out an error in caleu- 
lating chi-square for the difference, 


| observed — expected | — 1/2. 


The difference should be 48.8, with the correct value of chi-square, 41.2. 
Since this is for 1 degree of freedom, its square root, 6.41, is roughly a 
normal deviate and compares adequately with the value, 6.13, obtained 
by the more direct method. 
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LURIA, 8. E. and R. DULBECCO (University of Indiana). 
87 Interpretation of the Formation of Active Bacterial Virus From 
Ultraviolet Inactivated Virus (Genetics 34: 93-125). 


Bacteriophage particles are inactivated at a logarithmic rate by 
ultraviolet ight. The In of the survival ratio gives the average number 
of inactivating hits (r) per particle. When two or more inactive particles 
of the same bacteriophage are adsorbed by a bacterial cell, there is a 
production of active bacteriophage in a fraction of these cells. This 
fraction increases with the average number z of particles adsorbed per 
bacterium and decreases with increasing r._ These results are interpreted 
as due to the production of lethal mutations by ultraviolet in discrete 
genetic units, of which each particle contains a constant number n. From 
the experimental values of r, x, and of the probability that a bacter1um 
liberates active phage, a relation is derived between these parameters, on 
the assumption that active virus is produced in a bactertum whenever 
the infecting particles as a group contain at least one copy of each of the 
n units in nonlethal form. The theoretical relation fits the experimental 
results rather closely and leads to a calculation of a minimum value for 
n for each bacteriophage. The limitations of this analysis and some 
systematic deviations from experiment were pointed out and discussed. 


NEWMAN, E. V.; M. MERRILL. (Johns Hopkins University). 
88 The Application of Equations Derived from Models, to ‘Central’ 
Circulatory Volume. 


The purposes of our studies are to (1) derive a theory which expresses 
the concentration change in the outflow fluid of a flow system such as the 
human heart and lungs after a single instantaneous injection of a known 
amount of an indicator substance such as T-1824, (2) test the theory by 
comparison to dilution curves obtained from mechanical models in which 
the flow and the volumes of the compartments are known, and (3) apply 
the theory to the analysis of dilution curves obtained with human 


subjects. 
169 
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An equation was derived which expresses the variation of concentra- 
tion with time as a function of the amount of dye injected, the rate of 
flow through the system, and the volumes of three chambers in series. 
The right heart, the lungs and the left heart are theoretically considered 
as the three separate successive volumes in which the dye is mixed and 
diluted. i 

The dilution curves obtained from a mechanical model are nearly 
identical with the theoretically derived curves. The equation gives the 
outflow-fluid concentration as the algebraic sum of three exponentials 
whose rate constants are made up of known constants consisting of the 
volumes in the system, the amount of dye injected, and the flow through 
the system. 

Comparison of the theoretically derived and mechanically produced 
dilution curves with human curves shows close similarity. The relation- 
ship of the constants derived from human curves to the volume of blood 
in the heart and lungs is discussed. These mathematical and mechanical 
models provide a basis for the rational interpretation of human dye- 
dilution curves. 

A device for rapid accurate collection of serial samples from flow 
systems has been constructed. 


89 DENSEN, P. M. (University of Pittsburgh). A Definition of the 
Group to be Followed (To appear in Human Biology). 


Follow-up studies in relation to morbidity of one kind or another 
have been appearing in the literature with imcreasing frequency as a 
direct consequence of the growing importance of the chronic diseases as 
causes of morbidity and mortality. A clear definition of the group to 
be followed is a necessary prerequisite to such studies because the 
objective is to obtain a set of facts which may serve as the basis for 
predictions about other similar cases. It is only to the degree that the 
study cases have been precisely defined that this similar group can be 
adequately described and the extent indicated to which generalization 
may be permitted. Several examples are given to illustrate this point. 

Definition from the standpoint of diagnosis is only a part of the 
total definition. Other factors besides diagnosis, such as age, sex, dura- 
tion of infection prior to diagnosis, etc., may play a very important part 
in the determination of the universe to which the findings may be gen- 
eralized. 

Often an investigator may have difficulty in finding material to 
work with of precisely the kind he needs and he may decide to take the 
best available. This is fair enough as long as he recognizes that the 


ABSTRACTS i 


universe to which generalization is made may not be the universe he 
originally had in mind. It is essential that the investigator specify the 
nature of the selection in order that others may make proper use of the 
findings. In this connection it must be recognized that a definition can 
be adequate only in relation to some specific objective. It is part of the 
job of the statistician to insist that a definition be made and that its 
implications be understood. It is not his job, however, to make the 
definition. 


HARRIS, T. E. (Rand Corporation), PAUL MEIER (Phila- 

90 delphia Tuberculosis and Health Association), and JOHN W. 
TUKEY (Princeton University). Timing of The Distribution of 
Events Between Observations.* 


When all that is known as to the time of occurrence of an event is that 
it had not occurred at one known time (‘last previous’’) and that it had 
occurred at another known time (‘‘first after’’), it is ordinarily useless 
to try to date the event more closely. But when we have information of 
this nature about many events, and when it is reasonable to think of 
the events as a sample from a distribution, it is possible to use the 
information to estimate the distribution. It is natural to express this 
distribution in terms of event-rates. The problem discussed here is how 
these event-rates may reasonably be estimated from such observations. 

Any attempt to infer the timing of events from sparse observations is 
subject to many pitfalls, some of which are discussed briefly. 

If a given body of data has avoided these pitfalls, a simple statistical 
model may apply, for which the best standard estimation procedures are 
available. The resulting maximum likelihood estimates are shown to 
correspond quite closely to those found by a very simple routine, namely: 


) 


(1) Each case between “last previous” and “first after” contribute 
one half to the number of cases exposed to risk, 

(2) Each case’s event is distributed over the interval between “last 
previous” and “first after’ according to the finally estimated 


chance of occurrence. 


(The appearance of the one-half and of these rules is a matter of mathe- 
matics, rather than the result of simplifying assumptions.) 

An iterative procedure of obtaining exact maximum likelihood esti- 
mates, and making rough (but usually adequate) significance tests is 
given, and the danger that maximum-likelihood-like estimates may fol- 
low the irregularities of the data too faithfully is discussed. 


*Prepared in connection with research sponsored by the Office of Naya] Research. 
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These procedures, and some modifications, are illustrated on a simple 
numerical example and on the “progression” of the cases of minimal 
tuberculosis. The latter example is based on the records of the Henry 
Phipps Institute as studied by Drs. A. L. Cochrane, H. W. Campbell 
and S. C. Stein. The present analysis confirms their conclusions, as 
obtained by simpler methods, and exhibits a marked tendency for the 
progression dates to avoid a period of 9 to 12 months after detection. 
The possible significance of this result is discussed. 


91 DORN, HAROLD F. (National Institute of Health). Methods of 
Analysis for Follow-Up Studies. 


Two practical difficulties arise in long time follow-up morbidity 
studies. (a) The impossibility, without excessive cost, of keeping each 
member of a group of cases under observation, resulting in only partially 
complete information, and (b) the necessity of combining, in a meaning- 
ful way, the experience of cases of varying durations. Three concepts 
essential in the analysis of data from such studies are defined: (a) cohort 
of cases, (b) person-time units of exposure and (¢) exposure to risk. 
Principles and procedures pertaining to the following problems are dis- 
cussed: (a) cases to be included, (b) starting date, (c) classification of 
status, (d) the length of the interval of tabulation, (e) time at which the 
status of an individual is to be determined, (f) handling of cases observed 
for part of an interval, and (g) handling of cases lost to follow up. 
Formulas for different types of rates are given and the arrangement and 
construction of a morbidity table is illustrated. 


WORCESTER, JANE (Harvard University) and STUART 8. 
STEVENSON (University of Pittsburgh). Malformations in the 
Boston Lying-In Hospital, 1930-1941. (To be published in 
Pediatrics). 


The gestational characteristics associated with the birth of 677 con- 
genitally malformed infants have been studied. The difficulties involved 
in the selection of control groups are discussed. 


93 YOUDEN, W. J. (National Bureau of Standards). Index for 
Rating Diagnostic Tests. 


Diagnostic tests have the task of correctly designating which are the 
diseased individuals in a population. The test may err in giving negative 
tests for diseased individuals (false negatives) and in giving positive 
tests for healthy individuals (false positives). An index, 
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~ (a+ de +4)’ 


is proposed as a measure of performance of a diagnostic test, where a 
is the number correctly diagnosed out of a + b diseased individuals, 
and d is the number correctly reported negative out of c + d controls. 
The index provides a ready means of comparing the merits of alternative 
diagnostic tests. 


94 GREENHOUSE, SAMUEL W. and NATHAN MANTEL (Na- 
tional Cancer Institute). The Evaluation of Diagnostic Tests. 


This paper presents statistical procedures for evaluating and com- 
paring diagnostic tests which yield a continuous range of scores for the 
known positive and negative individuals tested. These procedures are 
developed both assuming normality of the distribution of diagnostic 
test scores, and making no assumption about distribution forms. The 
evaluation procedures shown do not make use of any critical point which 
distinguishes between positives and negatives. This makes it possible 
for the statistician to make an evaluation unhampered by any prior 
judgements of the value of the critical point. 

For a single diagnostic test, the procedure for determining sample 
size necessary for evaluation is given. Criteria for selection of a critical 
point for differentiating between positives and negatives for a single test 
are also considered. 

In the comparison of two diagnostic tests, procedures are developed 
for the case where the same individuals are used for both tests, and also 
where different individuals are used for each of the tests. 


First InpIAN Recion BiomMeTRIC CONFERENCE 


Poona, January 1950 


95 DANDEKAR, V. M. Certain Modified Forms of Binomial and 
Poisson Distributions. 


The Binomial and Poisson distributions in which the probabilities at 
successive trials are independent have been slightly modified by intro- 
ducing the condition that the probability of occurrence is zero for a 
number of trials immediately following an occurrence. Appropriate 
distributions have been obtained and the occasions suitable for their 


application indicated. 


174 BIOMETRICS, JUNE 1950 


96 BANERJEE, BASUDEB and ANUKUL CHANDRA DAS. On 
the Response of Mimose-Pudica Leaflets to its Organ Extract. 


Ricca’s observation, that leaflets of Mimose-Pudica respond to a 
Chemical substance and it is present in its organ extract, led the botanists 
to isolate that Chemical substance. A purer sample of the extract will 
give response to a higher dilution. 

In this paper it has been established that for determining response, 
time at the beginning of reaction should be noted and time for its 
completion is not needed. The latter is constant for all concentrations 
and under all atmospheric conditions but the former varies. But whether 
the latter is constant for all degrees of purity is to be investigated further. 


97. RAO,S. RAJA. Normal Curve as an Approximation to Statistical 
Distributions. 


The probability integral tables of Statistical Distributions are used 
by Statisticians mostly in tests of significance. For this purpose all that 
we need to know is the value of the “Statistic” at specified percentage 
points like .5 per cent, 1 per cent, 2.5 per cent, 5 per cent, ete. 

In this note, combinations of 6, and 8, have been obtained such 
that 


/ F(a) dx | < 0.0025 where x, and x, are such that 


i i(e) de = iE o(x) dx = a (a given probability) and 


and F(x) is the Gram-charlier series of type A with the first five terms 
which can be taken as a good representation of moderately skew Dis- 
tributions. Under these conditions x, is a good approximation to 2, . 

Suitable graphs have been drawn to facilitate the use of this method 
in such approximations. Making use of this criterion, the size of the 
sample beyond which the sampling distributions can be considered as 
normal for purposes of tests of significance has been worked out for some 
important Statistical distributions like Student’s ¢, x’, ete. 


98 POTI, S. JANARDAN, Design and Analysis of Blood Group 
Sample Surveys. 


An attempt has been made to fix the sample size for detecting with 
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confidence any given differences in the gene frequencies of two groups of 
individuals. ‘Tables for some typical values of the gene frequencies have 
been calculated. 


99 RAO, C. RADHAKRISHNA, The Distribution of the Difference 
Between the D° Statistics Based on p and p + q Characters. 

In a paper published in Sankhya, Vol. 9, Part 4, the author derived 
the distribution of D>,, — D; when the population value A? is zero. 
In this paper the distribution is derived for nonnull A? and some illus- 
trations have been given when the statistic based on the difference 
D;+. — D; is useful in tests of significance. 


ANNUAL MEETING OF THE British REGION OF THE BIOMETRIC SOCIETY 
Tue Mepicat Socirry or Lonpon, Marcu 14, 1950 
The Medical School of London 
March 14, 1950 


100 GRIDGEMAN, N. T. The Graphical Calculation of the Results 
of Biological Assays with Graded Responses. 


The fiducial limits of error, at a given probability level, of a balanced 
assay (response linearly related to log dose) can be expressed as 


IIRC -)+ VC — i(RC+ DY 


where J = log dose-interval, Rk = standard-test response difference as a 
fraction of dose-interval response difference, and C = Fieller’s correction 
factor. A chart has been constructed enabling percentage limits of 
error (100 X antilogs of the above expression) to be read off. 


10] WOOD, E. C. The Estimation of Error in Certain Types of 
Biological Assays. 


In those biological assays in which one animal from each litter is 
assigned to each treatment-group, the ‘residual’ component of variance 
after removing the Treatment and Litter components is used for estimat- 
ing the error mean square. If, however, two or more litter-mates can be 
assigned to each treatment-group, as occasionally becomes possible, an 
estimate of the true error mean square is available independently of the 
Treatment X Litter interaction. It is then sometimes found that the 
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latter, or some of the orthogonal components into which it can be par- 
titioned, is significantly greater than the error. The meaning of this 
finding, and its effect on the calculation of the fiducial limits of the 
result, are discussed. 


102 FIELLER, E. C. The Problem of Combining the Results of 
Independent Assay. 


The paper proposes methods for combining the data of independent 
assays with a continuous response: 


(i) When the slope and residual variance are stable from one assay 
to another, 
and 


(ii) When the slope varies while the residual variance remains 
constant, 


and discusses the application of these methods to litter-mate assays. 


Joint MEETING OF THE INSTITUTE OF MATHEMATICAL STATISTICS AND THE HASTERN 
Norru AMERICAN REGION OF THE BIOMETRIC SOCIETY 


Chapel Hill, North Carolina 


March 17-18, 1950 


103 GHURYE, S. G. (University of North Carolina). A Method of 
Estimating the Parameters of an Autoregressive Time Series. 


The general autoregressive process of the second order is defined by 
the equations 


x, = X, + Nt 
and 


Ay aX 4-1 ate OX 15 == E79 


where x, is the value actually observed at time t, X, the corresponding 
theoretical value, e, the disturbance and », the superposed variation. 
The estimates of a; , a» given by Yule’s method are biased and incon- 
sistent if , is not identically zero, the permanent bias being a function 
of the unknown variance of ,.. The present: paper proposes a method of 
estimation which is unaffected by the presence of 7, , and seems to be 
better than any other known method; and this conjecture is supported 
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by the results of application to observational and artificial series. In 
this method, the estimates a; , a, are obtained by mimimizing 


- i 
> CN" em 2) 


N-k 2 
os (et Oye. Geby-2) (Sage Oi 2i44 a autias-d} 
r= 


where n is some number small in comparison with N (which is the number 
of observations). In the above expression, the usual approximation of 
substituting 


N-k 
(N = & —-2)y, for Se ae: 
t=3 


may be made for computational convenience. The method has been 
used for fitting autoregressive processes to the series of annual averages 
of Wolfer’s sunspot numbers and that of Myrdal’s Swedish cost of living 
index numbers. The method is applicable to higher order processes. 


104 HOEFFDING, WASSILY. (Institute of "Statistics, University 
of North Carolina). Most Powerful Rank Order Tests. 


chee A ay = pe De random. variables with 
a joint probability function P(S), and let P{X,, = Xan} = ONG A 
hii = 1,--- ,k). Let H, bea hypothesis which implies that P(S) is 
invariant under all permutations of X;,, °°: , Xm(@ = 1, °°: , k). 
Let r;;(7 = 1, --- , m:) be the ranks of X,,, °°: , Xin, - Under Hy the 
M = Iin, : rank permutations R = (mir, °°" ) Tint) °°" 9 Tet °° Ste) 
have the same probability Pik) = M 1 A test which: depends only 
on the permutations £ is called a rank order test (R.O.T.). A R.O.T. of 
size m/M which is most powerful (M.P.) against a simple alternative, 
P,(S), is determined by m permutations F for which P,(/) takes on its 
m largest values. 

For example, let the pairs (X1 , Yi), «°° ; (X, , Y,) be independ- 
ent and identically distributed. Let Ho state that X; , Y; are inde- 
pendent, and let H,(p) be the hypothesis that X,; , Y; have a bivariate 
normal distribution with correlation p. We may assume that. 

. < X, and consider the ranks 7; of the Y’s only. A R.O.T. which 
is uniformly M.P. against all H,(p) with p > 0 does not exist except 
for small n. The M.P.R.O.T. against small p > 0 is determined by 
the largest values of doi-1 (EZ;)(EZ,,), where EZ, is the expectation 
of the i-th order statistic in a sample of n from a standard normal 
distribution. The M.P. unbiased R.O.T. against small values of | p | 
‘5 based on the statistic >>; >»; (BZ:Z;)(EZ,.Z,;). The M.P.R.O.T. 
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against p close to 1 is obtained by expanding the probability of R in 
powers of {(1 — p)/(1 + »)}" 


COCHRAN, WILLIAM G. (Department of Biostatistics, Johns 
105 Hopkins University). The Comparison of Percentages in Matched 
Samples. 


In this paper the familiar x” test for comparing the percentages of 
successes in a number of independent samples is extended to the situation 
in which each member of any sample is matched in some way with a 
member of every other sample. This problem has been encountered in 
the fields of psychology, pharmacology, bacteriology, and sample survey 
design. A solution has been given by McNemar (1949) when there are 
only two samples. 

In the more general case, the data are arranged in a two-way table 
with r rows and ¢ columns, in which each column represents a sample 
and each row a matched group. The test criterion proposed is 


Cet PENCE Tet hy 

a SAS U;) 
where 7’; is the total number of successes in the j-th sample and wu; the 
total number of successes in the 7-th row. If the true probability of 
success is the same in all samples, the limiting distribution of Q, when the 
number of rows is large, is the x” distribution with (c — 1) degrees of 
freedom. The relation between this test and the ordinary y’ test, valid 
when samples are independent, is discussed. 

In small samples the exact distribution of Q can be constructed by 
regarding the row totals as fixed, and by assuming that on the null 
hypothesis every column is equally likely to obtain one of the successes 
ina row. This exact distribution is worked out for eight examples in 
order to test the accuracy of the x° approximation to the distribution of 
@ in small samples. The number of samples ranged from ¢ = 3 toc = 5. 
The average error in the estimation of a significance probability was 
about 14 per cent in the neighborhood of the 5 per cent level and about 
21 per cent in the neighborhood of the 1 per cent level. Correction for 
continuity did not improve the accuracy of the approximation, although 
it is recommended when there are only two samples. Another approxi- 
mation, obtained by scoring each success as “1” and each failure as “0” 
and performing an analysis of variance on the data, was also investigated. 
The F-test, corrected for continuity, performed about as well as the x” 
approximation (uncorrected), but is slightly more laborious. 


The problem of subdividing x’ into components for more detailed 
tests is briefly discussed. 


= 
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LUCAS, H. L. (North Carolina State College). A Method of 
106 Estimating Components of Variance in Disproportionate Sub- 
Class Data. 


By including sufficient effects in the forward solution of the Ab- 
breviated Doolittle method, components of variance may be estimated 
from disproportionate data. The procedure is very systematic, and thus, 
is adaptable to routine computational work. The computations will 
be described, and the utility of the method briefly discussed. 


ISAACSON, STANLEY L. (Columbia University). On the 
107 Theory of Unbiassed Tests of Simple Statistical Hypotheses 
Specifying the Values of Two Parameters (Preliminary Report). 


In the Neyman-Pearson theory of testing simple hypotheses, in the 
one-parameter case, a locally best unbiassed region is called “type A.” 
It is obtained by maximizing the curvature of the power curve at the 
point @ = 4% specified by the hypothesis, subject to the conditions of 
size and unbiassedness. For the two-parameter case, Neyman and 
Pearson considered “type C” regions (Stat. Res. Mem., vol. 2 (1938), 
pp. 36, ff.). The definition of these regions requires one to choose in 
advance a family of ellipses of constant power in an infinitesimal neigh- 
borhood of the point (6, , 92) = (6: , 62) specified by the hypothesis. 
The natural generalization of a “type A” region is a “type D”’ region, 
which maximizes the Gaussian curvature of the power surface at (6) , 62), 
subject to the conditions of size and unbiassedness. This definition does 
not require one to choose a family of ellipses in advance. This approach 
leads to a new problem in the calculus of variations. A sufficient con- 
dition is obtained which plays the role of the Neyman-Pearson funda- 
mental lemma in the “type A” case. An illustrative example is given. 
(Prepared under sponsorship of the Office of Naval Research. ) 


108 BOSE, RAJ CHANDRA. (Institute of Statistics, University of 
North Carolina). A Note on Orthogonal Arrays. 


Consider a matrix A = (a,;) with N rows and m columns, each 
element a,; standing for one of the s — 1 integers 0, 1, 2, --- , 8 — if, 
Let us take the partial matrix obtained by choosing any ¢ < m columns 
of A. Each row now consists of an ordered ¢-plet of numbers, and each 
element has one of s possible values, there are s’ possible é-plets. The 
matrix A may be called an orthogonal array (N, m, s, t) of size N,m 
constraints, s levels and strength 1, if by choosing any ¢ columns whatso- 
ever every possible ‘-plet occurs the same number of times. Clearly 
N = Xs‘ where d is an integer. Such arrays have been considered by 
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Rao and are useful for, various experimental designs. The existence of 
an orthogonal array (s°, m, s, 2) is equivalent to the existence of a set of 
orthogonal Latin squares of side s and m constraints (i.e. the number of 
Latin squares in the set ism — 2). The fundamental question that can 
be asked regarding orthogonal arrays is the following: What is the maxi- 
mum numbers of constraints for an orthogonal array, given N, s and ¢? 
Denote this number by f(N, s, ¢), then from known properties of Latin 
squares f(s’, s, 2) = s + 1, if sis a prime or a prime power, and a theorem 
by Mann states that f(s’, s,2) >r +1, ifs = pi’ --- pi’ where p,, --- , Dx 
are different primes, and r is the minimum of po’ , pi’, -:: p. % The 
following generalisation of Mann’s theorem is proved in this note. 


fWiN:2 pg Nx » $182 sti Shs tb) = Min F(N: » Si, i), TN Goh Sen cd ys 
ee NE Sent) 


FREEMAN, MURRAY F. and JOHN W. TUKEY (Princeton 
109 University). Transformations Related to the Angular and the 
Square Root.* 


The use of transformations to stabilize the variance of binomial or 
Poisson data is familiar (Anscombe, Bartlett, Curtiss, Eisenhart). The 
comparison of transformed binomial or Poisson data with percentage 
points of the normal distribution to make approximate significance tests 
or to set approximate confidence intervals is less familiar. Mosteller 
and Tukey have recently made a graphical application of a transforma- 
tion related to the square-root transformation for such purposes, where 
the use of ‘binomial probability paper’’ avoids all computation. We 
report here on an empirical study of a number of approximations, some 


intended for significance and confidence work, and others for variance 
stabilization. 


0 VERLINDEN, F. J. (North Carolina State College). Standard 
Inverse Matrices for Fitting Polynomials. 


For fitting polynomials of the type, y = box® + bw + box? + 
-++ + 6,2", with the x’s equally spaced, published tables of orthogonal 
polynomials may be used. This procedure does not yield the b’s directly, 
nor their variances or covariances, although such may be obtained by 
proper computations which are moderately tedius. In some types of 
statistical work, the b’s and their variances and covariances may be 
desired. These may of course be obtained directly by the method of 
least squares but the computational work is prodigious relative to that 


*Prepared in connection with research sponsored by the Office of Naval Research. 
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for the orthogonal polynomial approach. When the 2’s are equally 
spaced the elements of the variance-covariance matrix may be put in 
the simple form of sums of powers (including the zero power) of suc- 
cessive integers from zero to n (n equals one less than the number of 
observations). The elements of the inverses of matrices of this type have 
been worked out algebraically in terms of n for polynomials up to and 
including the quintic (m = 5). With these standard inverse matrices, 
the b’s and their variances and covariances may quickly be obtained 
once the elements are evaluated numerically. These elements have 
been evaluated numerically up to n = 20. 


RAFFERTY, J. A., M.D. (Dept. of Biometrics, School of Avia- 
111 tion Medicine, Randolph Field, Texas). Mathematical Models 
in Biology. 


From the point of view of a bio-medical research administrator, 
mathematical models will assume a greater role in biological research 
than heretofore. In anticipation of this trend, certain philosophical 
fmplications of models in biological theory and scientific theory in 
history are examined. A hierarchy of abstraction-levels in biology is 
delineated, and the role of mathematical models at these levels is illus- 
trated by examples from the literature. Proposals are made for a con- 
centration of mathematical effort on certain important biological prob- 
lems. Remarks are made on the capabilities and limitations of models 
in biology. 


112 BROSS, IRWIN. (Johns Hopkins University). Small Sample 
Performance of Biological Statistics. 


In this paper the dilution method for estimating bacterial density is 
investigated by an exact small sample method and also by an approxi- 
mate one. Methodologies and design of experiments are compared for 
various small sample cases. 


GREENBERG, B. G. and A. HUGHES BRYAN. (University 
113 of North Carolina). Methodology in the Study of Physical 
Measurements of School Children. 


In aseries of investigations to determine by small-sampling technique 
what physical differences, if any, occur between children of differing 
socio-economic backgrounds, several problems of methodology arose. 
A pilot study was undertaken to assure maximum efficiency at each step. 
It was found that the children could remain dressed (with the exception 
of boys’ bi-iliac measurement) without changing the magnitude of the 
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differences. The pilot study enabled us to decide how many observers 
to use, and how much duplication of measurements by them was neces- 
sary. Minimum sample sizes were estimated to indicate physical dif- 
ferences of predetermined magnitudes. It was found that the age group- 
ing 96-143 months was optimal from the standpoint of indicating 
physical differences between children of differing socio-economic levels. 
Boys and girls in the upper socio-economic levels were both taller and 
heavier for their age in this age group. There were no weight differences, 
however, when weight was adjusted for age and height. Measurement 
of the bi-iliac and transverse chest diameter provided little additional 
information on physical differences. The calf circumference, an indicator 
of muscle mass and subcutaneous fat, is suggested as being a sensitive 
supplementary index to indicate physical differences when age and height 
are adjusted. 


HOUSEHOLDER, A. S. (Oak Ridge National Laboratory). 
114 ne 
Tetrad Analysis in Yeast. 


In neurospora all four products of meiosis are recovered in the four 
spores of an ascus. In crosses AB X ab the asci are of three types, 
designated I, II or III according as all four, none, or two spores resemble 
parents. Frequencies of these types, P, P’ and P” are the observables. 
If there were no exchange P’’ would be zero; and one should have P’ = 0 
or 1/2 according to whether the loci were on the same or different 
chromosomes. 

Assuming only that no exchange occurs between sister chromatids 
and neglecting chromatid interference, one can calculate without further 
assumptions a frequency P” of exchanges between a single locus and its 
centromere from data on three or more genes taken in pairs by equations 


SF ="S0sso7 PU ae e OTe ts iS! 


where the subscript 0 refers to a centromere. Lindegren makes such 
calculations from his own data, by taking groups of three, but makes no 
effort to reconcile discrepancies. Neyman’s modified chi-square, how- 
ever, permits combining all observations in a set of equations that yields 
easily to a rapidly converging iterative solution. The equations are 
28; DS si(nis + mii) (nz + ni) = Dy s(n; + ni) Qn. = mij ) 


Aa ji 


where n;; is the number in class I and II combined for the loci ¢ and pe 
ni; the number in class III, and only those pairs (2, 7) are included which 
are found to be independent. 


The argument of A. R. G. Owen (Pr. R.S. B. 136, ’49), p. 67 can be 
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paraphrased for the present case and a suitable generating function 
PAX, w) is being sought providing a metric. The specific one proposed by 
Owen is ruled out since 


s = P(-3,w 


takes on a negative value for one locus, which is not possible with Owen’s 
function. 


RAPOPORT, ANATOL. (University of Chicago). Contribution 
115 to the Probabilistic Theory of Neural Nets. I. Randomization 
of Refractory Periods and of Stimulus Intervals. 


Agaregates of neurons are considered in which the frequency of 
occurrence of neurons with a specified value of the refractory period 
follows certain probability distributions. Input-output functions are 
derived from such aggregates. In particular, if input and output in- 
tensities are defined in terms of stimulus frequencies and firing fre- 
quencies per neuron respectively, it is shown that a rectangular distribu- 
tion of refractory periods leads to a logarithmic input-output curve. 
If input and output are defined in terms of the total number of stimuli 
and firings in the aggregate, it is shown how the “mobilization” picture 
leads to the logarithmic input-output curve. 

By randomizing the intervals between stimuli received by a single 
neuron and by introducing an inhibitory neuron a very simple “filter 
net” can be constructed whose output will be sensitive to a particular 
range of the input, and this range can be made arbitrarily small. 


LANDAHL, H. D. (University of Chicago). Theoretical and 
116 Experimental Aspects in the Removal of Air-Borne Matter by 
the Human Respiratory Tract. 


The principal factors governing the fate of a particle in the respitatory 
tract are impaction due to inertia, settling due to gravity and Brownian 
movements. For a given respiratory pattern, it is possible to calculate 
the probable fate of a particle from a knowledge of the geometry of the 
passages. These calculations have been carried out in such @ manner as 
to obtain the theoretical amounts of material deposited in various regions 
of the lungs as well as the relative amounts in various fractions of the 
expired air. Similarly, it is possible to estimate the probable fate of a 
particle which passes through the nasal passages. Experiments have 
been carried out to verify a number of these predictions. On the whole, 
the agreement, as illustrated in the slides, is fairly satisfactory when one 
considers the complexity of the calculations. 
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Wl WADLEY, F. M. (Navy Department). An Application of Bio- 
metrics to Zoological Classification. 


Statistical problems in taxonomy are discussed; attention must be 
paid to variation of individuals as well as of group means. Covariance 
analysis and the discriminant function technique are applied to multiple 
measurements in groups of molluscan fossils. 


MOSHMAN, JACK. (United States Atomic Energy Commis- 
118 sion). The Analysis of Hemotological Effects of Chronic Low- 
Level Irradiation. 


Several methods are investigated for analyzing the possible effects of 
chronic low-level irradiation upon the employees of the operating con- 
tractors of the US AEC. The effects investigated are those on the red 
blood count, hemoglobin, white blood count, lymphocytes and neutro- 
phils. The analysis includes measurements of significant differences 
among individuals, geographic sites and the exploration of various indices 
of exposure to radiation. A non-parametric determination of trend 
values for individuals which may be applied to mass data is considered. 


119 CURETON, EDWARD E. (University of Tennessee). Statistical 
Problems in Psychological Testing. 


Though great progress has been made in mathematical statistics in 
recent years, a number of the major statistical problems encountered 
in the development and use of psychological tests remain unsolved. 
Some of these problems are outlined, with particular reference to the 
mathematical models and assumptions implied by psychological theory, 
by the nature of the experimental data, and by the conditions under 
which the results and findings are to be applied. 


NICHOLSON, GEORGE E., JR. (Institute of Statistics, Uni- 
120 versity of North Carolina). Accuracy of a Linear Prediction 
Equation in a New Sample. 


The problem considered is as follows. Given two samples S, and S, of 
N, and N; observations on a p + 1 character random variable (y, x, , 
"++, @). Let Y, and Y; be the linear regression equation computed by 
the method of least squares from each sample. The effect of using Y, to 
predict the y’s in S, is considered. The ratio 


Fl D GoeeY a) a 


k 
S(y2 — Y.) 


ABSTRACTS 185 


is used as a measure of the predicting efficiency of Y, in S, relative to Y, 
when the X; are fixed for the usual regression model. The general 
multivariate case is also considered. 


BAHADUR, RAGHU RAJ. (Institute of Statistics, University 
121 of North Carolina). Smallest Average Confidence Sets for the 
Simultaneous Estimation of K Normal Means. 


Let v = (Gi; , -°* , Tint 3 *** 3 Ter, *** » Linz) Genote the combined 
sample point in samples of sizes n, , n2, +--+ , n, from normal populations 
™,,%2,°** , 7, respectively, 7; having mean py; and variance o;. Writing 
Km = (ui, M2, *** , He), denote the k dimensional Euclidean space of 
all points w» by R. Given any parameter point (u, o), where o = 
(0, , 02, °°* ,o,), and any set valued function f(v) defined for all sample 


points v and having subsets of F# as its values (which satisfies certain 
measurability hypotheses), let a(f/u, ¢) = probability of the statement 
“uef(v)” being false, and B(f/u, «) = expected Lebesgue measure of f(v). 
We consider the problem of constructing f(v) so as to make both a and 8 
‘as small as possible.” One of the results obtained is as follows. 

Given p, 0 < p < lI, let fr:t@ (0) Bit De ni[(z: — ws)/li)< 
f(p) - oi nils./1P}, where, = nz* DUN wii, 85 = ns" Dor’ wis — Bs)’, 
= (l,,l,, --- , i), the l,’s being given positive constants, and ¢(p) 
being determined by P(x: > ¢(p)-xw-.) = p, where x; , xw-» are inde- 
pendent chi-square variables with k, NV — k degrees of freedom (k < N = 
emt). Then (a) obviously a(f}:;) / u, CA) = p for all » and all c,o < 
¢ <o, and (b) if f(v) is any other function such that a(f / u,cd) < p 
for all » and all c, either (i) f(v) and fy-;:,)(v) differ by a set of measure zero 
for almost every 2, or (ii) supyer {B(f / u, cA)} >supyer {B(f:r~ / B, Or) 
for every c. 


KAWADA, YUKIYOSI. (Tokyo University of Literature and 
122 Science). Independence of Quadratic Forms in Normally Corre- 
lated Variables. 


An extension is given of theorems of Craig, Hotellmg and Matérn 
which includes the following theorem, proved by a new method: If 
two quadratic forms Q, , Q, in normally and independently distributed 
variates with zero means and unit variances satisfy the four conditions 
E(QiQi) = E(Qi)E(Q3), for 7,7 = 1, 2, then the product of the matrices 
of the two forms in either order is zero. 


VORA, SHANTILAL AMIDAS. (Institute of Statistics, Uni- 
123 versity of North Carolina). Bounds on the Distribution of Chi- 


Square. 
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Let 
j=k j=k . 
= > 0; — 1p)? /np,, x” = 2; +4 — Np)’ / Np; 
j=1 j=l 
where 


k k 
Vig aU) oo, = pe i0; > ype and N=n+k/2. 
1 1 


Bounds on the multinomial probability 7 in terms of x” are obtained. 
A triangular transformation of 


a, = 0; + 3 — Np,)/{Np. — Dts 
(i = 1,--- ,k — 1) to y; is applied so that 


t=k-1 


dx’ = pe, yi; ) 


where d is determined later by equating the coefficients of x *. Certain 
rectangles r(v) with (y; , +++ , Yx-1) aS a mid-point are non-overlapping 
and cover the entire space R,_, for v; = 0, +1, +2,:--. Ifx’? <e, 
then bounds on 7’ in terms of the integral of the (k — 1) dimensional nor- 
mal frequency function over the rectangle r(v) are obtained. Prob. 
{y’” < ec} isthe sum of T over x” < c, so the integral over the sum of rec- 
tangles whose mid-points lie within the hypersphere x” < cis considered. 
Two hyperspheres, one which contains the sum of those rectangles, and 
one which is contained in it are used for the bounds, giving 


Noe Hea (cs) —< Prob. {x”” << c} & Aa Fy -i(Gi), 


where F’,_;(«) is a chi-square distribution function with (k — 1) degrees 
of freedom and ), , dz , C1 , C2 are functions of c, n, k and p, , ++: , p, . 
Asn ©, both bounds tend to F,_,(¢c). Bounds of the same form are 
obtained for Prob. {x” < C}. Closer bounds for Prob. {x’ < C} are 
given in terms of a non-central chi-square distribution. 


1 HENDERSON, C. R. (Cornell University). Estimation of Genetic 
Parameters. 


Many applications of genetics and statistics to the improvement of 
plants and animals deal with experimental data for which the under- 
lying model is assumed to be 


Pp qd 
Yeu DL bees, Ue eee 
i=1 t=1 


where b; are unknown fixed parameters, x;, and z;, are observable 
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parameters, the wu; are a random sample from a multivariate normal 
distribution with means zero and covariance matrix | o:; || , and the 
€, are normally and independently distributed with means zero and 
variances og . If o;; = 0 when? # j and if «2 = o2,, the model is the 
one usually assumed when components of variance are estimated. 

Three different estimation problems are involved, (1) estimation of 
b; under the assumptions of the model, (2) estimation of u; and (3) 
estimation of ¢;; . The first two problems are not solved satisfactorily 
by the least squares procedure in which the wu; are regarded as fixed, but 
the maximum likelihood solution does lead to a satisfactory estimation 
procedure. 

Assuming that the ¢;; and o% are known, the joint maximum likeli- 
hood estimates of b; and wu; are the solution to the set of linear equations, 


> b,( ee Rey ih Oa) +e on u;( ye Piekea TO) 


‘= 


= > LiaYo/o, h= il 72° 5p 
Doe’ Cetin ete AG Made ceatual o) 
1=1 a t=1 a 


= > AhaYalon h=1,-*:,q 


Some important applications of this estimation procedure to genetic 
studies are described and certain computational short-cuts are suggested. 

The problem of estimating ¢;; has not been solved satisfactorily 
although under certain quite general assumptions the equations for the 
joint estimation of b; , u; , o;; , and o% can easily be written. The 
solution to the equations, however, is too difficult to make the procedure 
practical. Nevertheless unbiased estimates of o;; can be obtained by 
equating to their expected values the differences between certain reduc- 
tions in sums of squares computed by least squares and solving for the 
o;; - In general, the expectation of the reduction due to 


Cesare bee «ie < ) is EY), 
gh 


where d’” are the elements of the matrix which is the inverse of the 
(p + k)’ matrix of coefficients and the Y, are the right members of the 
least squares equations. 


COHEN, A. C., JR. (University of Georgia). Estimating the 
125 Mean and Standard Deviation of Normal Populations from 
Double Truncated Samples. 
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The method of maximum likelihood is employed to obtain estimates 
of the mean and standard deviation of a normally distributed population 
from double truncated random samples. Two cases are considered. 
In the first, the number of missing variates is assumed to be unknown, 
In the second, the number of missing (unmeasured) variates in each tail 
is known. Variances for the estimates involved in each case are obtained 
from the maximum likelihood information matrices. A numerical ex- 
ample is given to illustrate the practical application of the estimating 
equations obtained for each of the two cases considered. 


KALLIANPUR, GOPINATH. (Institute of Statistics, University 
126 of North Carolina). Minimax Estimates of Location and Scale 
Parameters. 


If the joint frequency function of the random variables X, , --- Xy 
contains only a scale parameter and is of the form 


re ey de 


then under mild restrictions the following theorem is proved. 
Theorem 1. If the loss function is of the form W[(a@a — a@)/a], the best 
or minimax estimate a,(a) minimizes the integral 


= Gi = & 1 xy ae =) 
I w( a )4 02, a ae 


and further 


Qo (ux Teen te Uy) = MOL (X; Pi Wy ots mall UD 


When both location and scale parameters are present and the joint 
frequency function is of the form 


1 es ta—#) 
Paige 2. ee Ne Sgr es ae a ommend 5 
a 


Qa Qa 


then (under conditions similar to those in Th. 1) one of the results 
obtained is 

Theorem 2. If the loss function is of the form W[(@ — @)/a], the best 
estimate (x) of @ minimizes the integral 


ote (RON Fors 6 Ly — 6 
i w( a ) Sof a 5 es) - ) a0 de 


and 
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These theorems have been applied to derive minimax estimates in 
the case of standard distributions. Finally, the problem of estimating 
the difference between the location parameters of two populations is 
briefly considered. The results obtained in this paper are a continuation 
of the line of approach suggested in Theorem 5 of Wald’s “Contributions 
to the Theory of Statistical Estimation and Testing Hypotheses.” Ann. 
Math. Stat., Vol. 10, 1939. 

The present work was carried out under ONR contract. 


ROY, S. N. Unstitute of Statistics, University of North Carolina). 
On Some Features of the Neyman-Pearson and the Wald 
Theories of Statistical Inference, Their Interrelations and Their 
Bearing on Some Usual Problems of Statistical Inference. 


127 


With two alternative hypotheses H, and H; it is shown that (i) the 
most powerful test of H, with respect to H, is automatically an unbiassed 
test in the sense that its power is never less than (and usually greater 
than) the level of significance a and (ii) there is also a least powerful 
test with its power not greater (usually less) than a. This means that 
all tests have powers lying in between, which gives a complete picture 
of the possible family of tests and provides a basis for defining efficiency 
of tests. 

With the first kind of error a is tied up a minimum second kind of 
error 8 (complementary to the maximum power P), and the level at 
which a is fixed depends upon some compromise between a and 8. 
This intuitive approach is formalised by the introduction of loss functions 
related to, and apriori probability weights for H, and H,, thus leading 
to the first stage in the Wald treatment of dichotomy with two solutions 
in the observation space corresponding respectively to minimum and 
maximum total risks. This is immediately generalised to the first stage 
in the Wald treatment of multichotomy with minimum and maximum 
total risk solutions. An important special case is discussed in which 
all the possible alternatives to a particular hypothesis are, by our test 
procedure, indistinguishable among themselves, thus effectively forming 
only one alternative to the hypothesis, which means a degenerate 
multichotomy. The bearing of this on most powerful tests on an average 
under the Neyman-Pearson theory is also discussed. 

The problem of testing of composite hypothesis which is usually 
treated in terms of the Neyman-Pearson theory is posed and treated in 
terms of the (first stage) Wald theory and an indication is given of ho 
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these notions could be applied to the usual problems of univariate and 
multivariate analysis. 


128 DAVIS, R. C. (U. S. Naval Ordnance Test Station, Inyokern, 
Calif.). Note on Uniformly Best Unbiased Estimates. 


For the estimation in an absolutely continuous probability distribu- 
tion of an unknown parameter which does not possess a sufficient sta- 
tistic, it is shown that no unbiased estimate for the unknown parameter 
exists which attains minimum variance uniformly over a parameter set 
of arbitrary nature. This result demonstrates the impossibility of ob- 
taining a generalized sufficient statistic first proposed by Bhattacharyya. 
Although not used in this note it is surmised that Barankin’s powerful 
results on locally best unbiased estimates can be applied to yield further 
results in this direction. 


ROBBINS, HERBERT E. (University of North Carolina). 
129 sates 
Competetive Estimation. 


Let 6 be a vector random variable with distribution function G(@) and 
let x be a vector random variable whose frequency function f(x; 6) 
depends on @. Two statisticians, A and B, are required to estimate 6 
from the value of «. If A’s estimate is closer to 6 he wins one dollar 
from B and vice versa; in case of a tie no money changes hands. It is 
shown that A should estimate @ by the function a(z) = median of 
posterior distribution of 6 given x; his expected gain will then be > 0 
whatever estimate B may use. If G(@) is not known to A he should 
estimate it from the series of values of @ which have been observed in 
previous trials. If these are not known, A should estimate G(6) from 
the values of « which have previously occurred; how this may be done is 
discussed elsewhere (see following Abstract). 

From the point of view of the theory of games, when G(@) is unknown 
we have a game in which the “rules”? are unknown and must be suc- 
cessively estimated from past experience. Other examples arise when- 
ever a game involves random devices whose probability distributions 
are not known to the players but must be inferred by statistical methods, 
in general from secondary variables which contain only part of the total 
information. The réle of statistical inference in such “long term’ games 
is fundamental. 


CHAND, UTTAM. (Department of Mathematics, Boston Uni- 
130 versity). The Effect of an Unknown ‘Location Disturbance’ on 
‘‘Student’s” t Based on a Linear Regression Model. 


ABSTRACTS 19] 


Consider y; , -+- yw. ,Y¥v, + 1, +++ yy, aset of observations ordered in 
time. If the y’s are normally and independently distributed according to 
N(a@ + B(t — 1), o) and we want to find out if the y’s have changed with 
time, we usually employ a “Student” ¢ type of statistic with N — 2 
degrees of freedom. If, as a consequence of the impact of a certain un- 
known political or economic change in the past on the y’s, the y’s actually 
constitute two independent, normal samples y; , --+ yy. , Yn, t 1,°°+ yy 
distributed according to N(m, , 0°), N(m 2, o°) respectively, a two-sample 
“Student” ¢ also based on N — 2 degrees of freedom would be the 
appropriate statistic to use for the hypothesis m, = m.. If, in fact, the 
latter situation describes the correct state of affairs, and the statistician 
employs the ‘‘Student”’ ¢ based on the regression model, he commits an 
error. The present paper investigates the nature of such an error in the 
light of the point of impact as determined by the magnitude of N, and 
the intensity of the impact as determined by the standardized 


‘distance’ = eee 
1 i 
o ~ eS 
N, N aa N, 


of this extraneous ‘shock’ on the ordered set of observations y. 


BRADLEY, RALPH A. (The University of North Carolina and 

131 McGill University). Corrections for Non-Normality for the Two- 
Sample t and the F Distributions Valid for High Significance 
Levels. 


The effects of non-normality of the parent population on common 
tests of significance have long been of concern in the application of 
statistical methods to experimental data. In this paper, the two-sample 
{ statistic is expressed as a simple multiple of the cotangent of an angle 
between two lines in a space of dimensionality one less than the total 
of the sample sizes; the F statistic for k samples is expressed as a multiple 
of the cotangent of an angle between a line and a plane of (k — 1) 
dimensions in a space, again, of dimensionality one less than the total 
of the sample sizes. The geometrical formulation is such as to suggest 
approximations to the distributions of these statistics valid for large 
values of the statistics, and these approximations are obtained. The 
approximations are shown to be exact in the special cases where the 
parent population is normal, and a method of evaluation of correction 
factors is given for a wide class of parent populations. The approxima- 
tion procedures are valid for the distributions under both null and 
non-null hypotheses. 


192 BIOMETRICS, JUNE 1950 


HANNAN, JAMES F. (Institute of Statistics, University of 
132 North Carolina). Some Tests Based on the Empirical Distribution 
Function (Preliminary Report). 


Let X = (X,, X.,--: , X,) be an independent sample of n where X; 
has the continuous ¢.d.f. F(x). Let S,(x) be the empirical distribution 
function. Acceptance regions of the type {X: S,(x) < $(x) for all x} 
are considered for different specifications of @ and their probabilities 
evaluated. The method of evaluation consists in identifying the regions 
with regions defined in terms of the order statistics of a sample of n from 
the uniform distribution on the interval (0, 1). The result obtained for 
o(«) = F(x) + c/n 0 < c integral < nis used to provide a direct proof 
of the Kolmogoroff result 

lim P{[n'” sup (S,(2) — F(a)) < 2] = 1—-—e”, 
while that obtained for (x) = F(x) + t,0 < ¢ < 1 gives the exact ¢.df. 
of the statistic sup, (S,(a) — F(@)). 


WALSH, DR. JOHN E. (The Rand Corporation). On a Gen- 
eralization of the Behrens-Fisher Problem. 


Let m + n independent observations be available where it is only 
known that a specified m of them are from continuous symmetrical 
populations with common median » while the remaining n are from 
continuous symmetrical populations with common median y. This is 
the generalization of the Behrens-Fisher problem investigated; some 
tests and confidence intervals for 4 — vy which are valid for the generalized 
situation are presented. For definiteness, suppose that n < m. The 
procedure used is to subdivide the m observations (common median sp) 
into n groups of nearly equal size and form the mean of the observations 
for each group. Pair the n means with remaining n observations and 
subtract the value of each observation from the value of the mean with 
which it is paired. The resulting n values represent independent obser- 
vations from populations with common median » — ». Tests and confi- 
dence intervals for 1 — y are obtained by applying the results of “Appli- 
cations of Some Significance Tests for the Median Which are Valid Under 
Very General Conditions” (Amer. Stat. Assoc. Jour., Vol. 44, 1949, pp. 
342-55) to these n values. To measure the “information” lost by using 
the generalized tests when one actually has two independent samples 
from normal populations, power efficiencies are computed with respect 
to: (a) Scheffe’s “best’’ ¢-test solution and (b) Most powerful solution 
when ratio of variances is known. Case (a) yields an upper bound while 
case (b) furnishes a lower bound for the actual efficiency. 
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SHRIKHANDE,S.S. Unstitute of Statistics, University of North 
134 Carolina). Construction of Partially Balanced Designs with Two 
Accuracies. 


Various methods of construction of Partially Balanced Designs first 
introduced by Bose and Nair (Sankhyd, 4 (1939), pp. 337-373) have 
been considered. Two of the methods given are generalisations of a 
Difference Theorem given by them. Another method is the inversion 
of an unreduced Balanced Incomplete Block Design with k = 2. Use 
has also been made of the existing Balanced Incomplete Block Design 
in another direction. A number of designs can also be obtained by 
methods of Finite Geometries and especially by omitting a number of 
treatments and certain blocks from the complete Lattice Designs. Use 
of curves and surfaces in Finite Geometries and the use of multifactorial 
designs given by Plackett and Burman (Biometrika, 33 (1946), pp. 
305-325) are also indicated. 


SHRIKHANDE, 8S. 8S. (Institute of Statistics, University of 
135 North Carolina). Designs for Two-Way Elimination of Heter- 
ogeneity. 

Use has been made of the existing Balanced and some Partially 
Balanced Designs for two-way elimination of heterogeneity with at most 
two accuracies. Particular cases of these designs were given by Youden 
(Contributions from Boyce Thompson Institute, EX (1937), pp. 317-326) 
and Bose and Kishen (Science and Culture (1939), pp. 136-137). The 
method depends upon interchanging the positions of various treatments 
in the different columns (blocks), if necessary, so as to satisfy certain 
conditions. 


SHRIKHANDE, S. S. (Institute of Statistics, University of 
136 ~ | ; 
North Carolina). Designs for Animal Feeding Experiments. 


In animal-feeding experiments change-over designs are generally pref- 
erable to continuous feeding experiments. In change-over designs both 
the direct and carry-over treatment effects are important. Use of 
Balanced and Partially Balanced Incomplete Block Designs toward this 
end has been considered. 


SANDELIUS, D. MARTIN. (Statistiska Institutionen, Uppsala 
Universitet). A Truncated Sequential Procedure for Interval 
Estimation with Applications to the Poisson and Negative Bi- 
nomial Distributions (Preliminary Report). 

Let z, ¥: , Yo, °°: be a sequence of random variables defined in 
(0, ©), and let be the smallest integer satisfying itt ys > tx, where 
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i > 0 isanon-random quantity. Define w, either as aoe y;/« or as the 
smallest integer exceeding )0'-, y;/a, k = 1,2, ---. Given the distribu- 
tion function F(a, 6) of x and, for any ¢, the conditional distribution of 
n with respect to 2, the distribution of wu, is obtained. The problem is 
to determine a confidence interval for 6 with confidence coefficient 1 — a 
on the basis of either an observation on uw, , if u, < t, or an observation 
on n, ifn < k — 1. The following procedure is proposed: If u, < ¢, 
choose 6,) and 6,, according to a rule satisfying Prob (A < @ < @11| uw < 
it) >1—a. Ifn < k — 1, choose @ and @, such that Prob (x < 
6<6,|n<k-—1) >1-—a. For continuous u, the following cases 
are discussed: A) a = @ with probability 1, and n has, for any ¢, a 
Poisson distribution with mean 46, B) x has a Gamma distribution with 
mean 6, and the conditional distribution of n with respect to x is, for 
any t, a Poisson distribution. Both cases may, for instance, be applied 
to bacterial counting. 


ROBBINS, HERBERT E. (University of North Carolina). A 
138 Generalization of the Method of Maximum Likelihood: Estimat- 
ing a Mixing Distribution (Preliminary Report). 


Let 6 be a vector random variable with distribution function G(é@) 
belonging to some class G, let x be a vector random variable whose 
frequency function f(x; @) depends on 6, and let g*(x) = Jf f(x; 0)dG(6) 
be the resulting frequency function of z. From a sample 2 , x, --- it 
is required to estimate G(@). The generalized method of maximum 
likelihood consists in using the estimates G,(@; 2, +++ , ,) in G for which 
[]t-: g*(@:) is a maximum. Under certain restrictions this method is 
consistent asn—o« . More generally, if the distribution function of x 
is of the form G*(a) = Jf F(x; 0)dG(6) and if this integral equation with 
kernel F(x; 6) defines a one-to-one continuous correspondence between 
G(@) and G*(x), then G(@) can be estimated by using the sample dis- 
tribution function of x, , +++ , x, to replace G*(x) and solving for the 
corresponding G(#). We can also apply the method of minimizing in G 
an appropriate measure of the deviation between the sample distribution 
function of 2 , ++: , % and G*(z). 

Any consistent method of estimating the mixing distribution G(é) 
from the sequence x, , #2 , «++ yields a solution of parametric statistical 
decision problems in the following manner: from past values 21, °** ,2no3 
we estimate G(@), and then use the corresponding Bayes solution of the 
decision problem to reach our decision for x, , even though the value 6, 
which produced «, is different from those which produced 2, , --- ceey 
In certain cases of long-term experimentation this approach seems more 
reasonable than the minimax method which decides on the course of 
action appropriate to 6, on the basis of x, only, and ignores the informa- 


tion about the prior distribution of @ which is contained in 2, pCR ae re 
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Officers for 1950. According to our constitution, the general officers 
are elected by Council each year, with an obligatory change in President 
every two years. In accord with these provisions the Council has 
elected Arthur Linder as President, J. W. Hopkins as Treasurer and 
C. I. Bliss as Secretary for 1950. Our new President is professor of 
mathematical statistics in the University of Geneva and in the Swiss 
Federal Institute of Technology at Zurich. Much of the success of the 
Second International Biometric Conference last summer in Geneva was 
due to the skill and tact with which he handled the many arrangements 
for the Conference. 

The following Council members were elected for the period 1950-52 
inclusive by mail ballot of the members of the Society: M. H. Belz, 
G. Darmois, R. A. Fisher, P. V. Sukhatme, O. Tedin and E. B. Wilson. 
J. W. Trevan was named to complete President Linder’s unexpired term 
as member of Council. Dr. Jane A. Russell and J. R. Wittenborn of 
Yale University served as tellers. 

Meetings in Stockholm. The International Union of Biological 
Sciences will meet this coming summer in Stockholm, Sweden on July 
7-11. Since the Biometric Society provides the secretariat of its Section 
on Biometry, we will be represented on the Executive Committee of the 
IUBS by President Linder and one other member of the Society. An 
International Congress of Botany on July 12-20, also in Stockholm, will 
follow the IUBS meetings. A biometrical program is planned in con- 
nection with these two meetings which we hope will aid in the develop- 
ment of a Scandinavian Region in the Society. 

BIOMETRICS. Since its formation, the Biometric Society has used 
BIOMETRICS as its journal, although ownership resided in the Amer- 
ican Statistical Association and its Biometrics Section. As an interna- 
tional organization, the Society has felt the need of having its own jour- 
nal. Because of the contributions to BIOMETRICS by the Society, 
both in content and subscriptions, the Council proposed to the American 
Statistical Association that the journal be transferred to the Society with 
appropriate safeguards for the interests of the members of the ASA who 
subscribe to BIOMETRICS but are not members of the Society. ‘Trans- 
fer of BIOMETRICS to the Society was approved successively by the 
ASA Biometrics Section, Board of Directors and Council at meetings in 
December 1949. President-elect Lowell J. Reed and the Chairman of the 
Biometrics Section, H. F. Dorn, were named as representatives of the 
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ASA to arrange the transfer. The Council of the Biometric Society has 
asked J. W. Hopkins, C. I. Bliss and Gertrude M. Cox to represent the 
Society in arranging the transfer. When negotiations have been com- 
pleted, details of the transfer will be published in BIOMETRICS. 

Proceedings of a biometrical-entomological clinic. A biometrical 
clinic for entomological problems was held in New York on December 
13, 1948, under the sponsorship of the Association of Economic Entomol- 
ogists and the Eastern North American Region of the Society. The 
proceedings were recorded electronically, transcribed and edited. The 
Council of the Society approved their being published through the 
Secretary’s office if a sufficient number of prepublication orders were 
received to cover the estimated expense. In response to notices sent to 
the members of both organizations, 98 orders were received from Society 
members and 489 orders from the economic entomologists. An edition 
of 700 copies was printed in February as a multilithed 64-page bulletin. 
As long as the supply lasts, copies can be obtained from the office of the 
Secretary at 50¢ for members of both sponsoring organizations and 75¢ 
for others. 

ENAR. The Eastern North American Region held its annual meet- 
ing in New York City on December 28-30 jointly with the Biometrics 
Section of the American Statistical Association and the American As- 
sociation for the Advancement of Science. At the Regional business 
meeting on December 30, the following officers were named for 1950: 
Vice-President, Joseph Berkson; Secretary-Treasurer, Walter T. Federer; 
Members of the Regional Committee from 1950 to 1952, Lila F. Knudsen 
and W. J. Youden. The scientific program consisted of three sessions. 
The first under the chairmanship of H. W. Norton concerned the use of 
rationally developed equations in biology, with papers by S. E. Luria 
and by E. V. Newman and Margaret Merrell. The second on long-time 
follow-up in morbidity studies was chairmaned by J. W. Fertig and 
included papers by P. M. Densen, by T. E. Harris, Paul Meyer and J. W. 
Tukey and by H. F. Dorn. The third session with Frederick Mosteller 
as chairman consisted of contributed papers by Joseph Berkson, by 
Jane Worcester and 8. 8. Stevenson, by W. J. Youden and by 8. W. 
Greenhouse and Nathan Mantel. 

The Biometric Society, E.N.A.R. and the Institute of Mathematical 
Statistics met jointly at Chapel Hill, North Carolina, March 17-18, 
1950. All sessions were joint meetings of the two organizations. Forty 
members of the Biometric Society were in attendance. 

Two papers due particular comment are the invited addresses on 
“Mathematical Models in Biology” by Dr. James A. Rafferty and on 
“Estimation of Genetic Parameters” by Professor C. R. Henderson. 
Dr. Rafferty stressed that much more cooperation among biologists, 
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mathematicians and statisticians is needed, if biological research is to go 
forward at the pace required by modern living. Professor Henderson 
outlined the estimation methods currently used by animal breeders, 
pointed out their difficulties and deficiencies, and asked the aid of 
mathematical statisticians in improving these methods. 

Highhghts of the meeting included a dinner on Friday evening at the 
Carolina Inn. Professor W. G. Cochran was toastmaster, and welcome 
was bid the two societies by Chancellor Robert B. House of the Uni- 
versity of North Carolina. Professor Gertrude M. Cox responded for 
the Biometric Society and Professor David F. Votaw for the Institute 
of Mathematical Statistics. 

At the close of the meetings, Saturday afternoon, the attendees were 
guests for tea at the home of Professor and Mrs. Harold Hotelling. 
Many visited the Morehead Building and enjoyed the spectacular show 
“Raster Awakening” given at the Planetarium. 

Professor Harold Hotelling was in charge of program for the Institute 
of Mathematical Statistics, and Professor H. L. Lucas in charge of 
program for the Biometric Society. Professor H. E. Robbins was in 
charge of arrangements and accommodations, and was assisted by 
Professor George E. Nicholson. These two should be commended for 
doing an heroic job. 

Two other regional meetings were held in April. The first was a joint 
session with the Society of Pharmacology and Experimental Thera- 
peutics at Atlantic City on April 19. Practical problems submitted by 
pharmacologists were discussed informally by a panel consisting of C. I. 
Bliss, W. G. Cochran, Frederick Mosteller and W. J. Youden, with 
E. G. deBeer as chairman. The second was a joint meeting with the 
American Mathematical Society at Oak Ridge, Tennessee on April 21, 
where the program included three papers on topics in bio-mathematics by 
J. Z. Hearon, by A. Rapoport and by C. W. Sheppard. 

Indian Region. The Indian Region held its annual meeting in Poona 
on January 5, 1950. The following were elected to serve on the governing 
Council for the year 1950: Vice-President, P. C. Mahalanobis; Secretary, 
C. Radhakrishna Rao; Treasurer, Anukul Chandra Das; Members, V. M. 
Dandekar, K. Kishen, K. R. Nair, U.S. Nair, V. G. Panse, P. B. Patnaik, 
B. Ramamurthy, F. N. Roy and P. V. Sukhatme. A scientific program 
included papers by V. M. Dandekar, by Basudeb Banerjee and Anukul 
Chandra Das, by S. Raja Rao, by 8. Janardan Poti and by C. Rad- 
hakrishna Rao. The meeting voted to invite the Biometric Society to 
hold an international session in India in the fall of 1951. 

Région Francaise. Région Frangaise held its annual meeting on 
February 28 at the Laboratory of Zoology of the Faculty of Sciences in 
Paris. The three ordinary members of the Regional Council were re- 
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elected, M. Lamotte, Mlle. Colette Rothschild and M. P. Schutzen- 
berger. The scientific sessions included three papers on problems of 
correlation, two of them by P. Rey and the third by M. P. Schutzen- 
berger. A second session of the Region late in April considered the 
probability of an all-or-none response as a function of a dose or other 
parameter. 

British Region. At the annual meeting of the British Region on 
March 14, the following officers were named: Vice-President, R. A. 
Fisher; Treasurer, A. R. G. Owen; Secretary, D. J. Finney; Regional 
Committee Members for 1950-1952, W. L. M. Perry and C. B. Williams. 
Abstracts of the papers read at the meeting appear in this issue of 
BIOMETRICS. 

The Interests of Members of the Biometric Society. We are indebted 
to Professor John W. Tukey of Princeton University for the following 
analysis of the widespread interests of our members. 

The directory published last summer lists the interests of 93% of 
the 900 who were members at that time. No interest was available for 
62 (830AR, 11BR, 7IR, 7 unattached, 3RF, 3ENAR, 3WNAR). A 
rough examination of the addresses and titles of these 62 suggests that 
their probable interests resemble those of the other 838, and so all 
members below are based on 838. Many members stated two, three 
or even four interests. ‘They have been prorated. Thus, Wilcoxon, 
who gave “insecticides, fungicides, herbicides, biometry”, was counted 
as 2/4 under ‘entomology, mycology, parasitology, nematology”’, 1/4 
under ‘‘agriculture’, and 1/4 under ‘‘biometry’’. There seemed to be 
no easily applied rule that would be fairer than this. (When information 
is gathered for another directory, it may be possible for each member to 
give better weighted information.) 

According to its letterhead, the Biometric Society is “An international 
society devoted to the mathematical and statistical aspects of biology”’. 
We should hope then, that most interests would be biological, and second 


would come mathematics and statistics. This is indeed the case, for 
we find 


Biological sciences Do, (441.4) 
Mathematics and statistics 22% (184 ) 
Bridge fields of biology 17% (139.8) 
Other areas 8% (672.3) 


Not only this broad breakdown, but also the more detailed distribution 
is of interest to members, officers and editors. Not all will agree with 


this classification, and so considerable detail is useful, since it permits 
ready rearrangement. 
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In the detailed breakdowns, the percentages are based on the totals 
of the broad classifications. 


BIOLOGICAL SCIENCES (441.4) 


22% Genetics (96 as follows: unspecified 42.3; plants 17.5; human 12; animals 10.5; 
poultry 4; microorganisms 1; quantitative 8.7). 

14% Bioassay and pharmacology (62.6 as follows: bioassay and antibiotics 35.8; 
pharmacology, therapeutics and toxicology 26.8). 

12% Medicine and public health (55 as follows: medicine, clinical research, ete. 14.3; 
surgery, pathology and specialties 16.7; hygiene and industrial health 8.5; 
public health 15.5). 

12% Physiology and nutrition (51.8 as follows: physiology 32; nutrition 19.8). 

11% Agriculture, forestry, fisheries (49.3 as follows: agriculture 22.8; forestry, horti- 
culture and pomology 16; fish and aquatic biology 10.5). 

8% Human specialties (34.3 as follows: psychology 22.8; human biology, anthro- 
pology and biotypology 11.5). 

15% Other specialties (66.7 as follows: evolution, population genetics and ecology 
17.5; entomology, mycology, parasitology and nematology 17.4; microbiology, 
bacteriology and virology 12.3; serology, hematology and immunology 7.3; 
cytology, embryology, anatomy and histology 6.2; herpetology and bird biology 
6). 

6% General biology (25.7 as follows: unspecified and applied 6.2; botany, zoology 
and systematics 10.3, quantitative, theoretical and mathematical 9.2). 


MATHEMATICS AND STATISTICS (184) 


92% Statistics (170 as follows: mathematical statistics and probability 71.8; unspeci- 
fied and applied 69.5; experimental design, inference and scientific method 
19.7; sampling 9). 

8% Mathematics (14). 


BRIDGE FIELDS OF BIOLOGY (139.8) 


56% Statistical (78.5 as follows: vital statistics, demography and actuarial 24.3; 
medical statistics 19.2; public health statistics 16.7; biostatistics 12.8; agri- 
cultural statistics and agricultural economics 5.5). 

24% Biometrics (33). 

20% Physical science (28.3 as follows: biochemistry 16; biophysics and radiation 
biology 12.3). 


OTHER AREAS (72.8) 


380% Engineering (21.3). 

28% Economics (20.5). 

16% Chemistry (11.5). 

11% Education ( 8.2). 

7% Physics and geophysics (5.3). 
8% Other sciences and technologies (6). 


In interpreting these tables, we must remember that a figure like “37 may repre- 
sent partial interests of 50, 60 or even 70 members. 


RECENT APPLICATIONS OF 
BIOMETRICAL METHODS IN GENETICS 


(1) EXPERIMENTAL TECHNIQUES IN PLANT IMPROVEMENT 


F. YatTrEs 


Rothamsted Experimental Station 
Harpenden, Herts, England 


WAS ONLY ASKED to contribute to this discussion a few weeks before 

the conference, and pressure of other work has prevented me from 
tackling certain investigations on which I hoped I might make some 
progress. The present paper, therefore, is more a statement of problems 
which require solution, than a presentation of any definite conclusions. 
The problems I want to consider are those which arise when planning a 
testing scheme which will be effective in testing the large number of 
new lines and varieties which are produced in the course of any pro- 
gramme of plant improvement. In the main I have followed the lines 
of a paper I gave to a Congress of Plant Breeders which was held at the 
John Innes Horticultural Institution last November. To those who were 
present at that meeting I must offer my apologies. I am, however, 
encouraged to repeat what I said there, because only a very brief 
summary of the meeting has been published (Lewis, 1949), and because 
the present audience is largely different. 

The reaction of the meeting itself was that the subject was of con- 
siderable importance and should be more fully discussed at the next 
Congress. Subsequent events showed, however, that the need for further 
discussion was not felt by all plant breeders, for another Congress was 
held this spring to which such troublesome characters as mathematical 
statisticians who might originate discussions of this type were not invited. 

The problems which I posed at the earlier meeting therefore remain. 
They have not been further discussed, and I personally have not made 
any progress with their solution. As far as I know no one else has either. 
But ignoring a problem does not get rid of it. I will therefore propound 
the problems again, in the hope that some here may be stimulated to 
making a contribution towards their solution. 
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DESIGN OF EFFICIENT TESTING MECHANISMS 


At all stages of the selection process we shall require to compare 
the various new lines and varieties that are produced in order to see 
which are sufficiently promising to be retained, either as a basis of further 
breeding work or as possible new varieties for commercial use. Quanti- 
tative characters such as yield must be tested in field trials, which may 
be regarded as the testing mechanism of plant breeding. As in all 
branches of science improvement of measuring instruments and testing 
mechanisms may be expected to result in progress in the science itself. 

The design of field trials for comparison of different varieties is a 
branch of experimental design and is the only aspect of the subject 
under discussion which has received much consideration by mathematical 
statisticians. Development during the inter-war period has been con- 
siderable. The introduction of the principle of randomization by Pro- 
fessor R. A. Fisher was the first and fundamental step—replication and 
local control (arrangement in blocks or other systematic patterns) had 
of course long been in use, though in Great Britain at least there had 
been a retrogression to large plots without replication. Randomization 
laid the foundation of further advances by providing an unequivocal 
means of estimating the experimental errors. This had a number of 
important consequences, one of which was that it became possible to 
judge the relative efficiency of different types of design in an objective 
manner. 

The second major development was the recognition of the importance 
of factorial design and the development of methods appropriate to 
elaborate factorial experiments—confounding, estimation of error from 
high-order interactions and fractional replication. The third major 
development, which was of particular interest to plant breeders, was the 
introduction of quasi-factorial and balanced incomplete block designs. 
These designs enable a large number of varieties (or treatments) to be 
compared in groups without the use of controls, thus permitting a greater 
degree of elimination of soil heterogeneity than would be possible with 
ordinary randomized blocks, while at the same time giving comparisons 
between pairs of varieties which are all of approximately equal accuracy, 
and computations which are reasonably simple. 

I do not wish to discuss the various alternative incomplete block 
designs at length here, but I would like to make a few general points. 
In the first place the type of design which is most appropriate depends 
on the nature of the tests we desire to make. It has been suggested, 
for instance, that since incomplete block designs are modern, anyone 
using designs involving controls is somehow behind the times. This, of 
course, is not necessarily,true. Under certain circumstances comparison 
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with standard varieties may be required; in such cases it may well be 
advantageous to use these as controls. Incomplete block designs are 
also unsuitable if a number of the strains to be tested are likely to fail 
completely, or if a number can be rejected by inspection without 
harvesting. 

Secondly there is an upper limit to the gain in efficiency that can 
be obtained by elaboration of the experimental layout. This limit is set 
by the inherent variability of the individual plants due to factors which 
operate substantially at random from plant to plant, and to the residual 
component of fertility variation which occurs with even the smallest 
plots. Probably the least variable type of element for a given size and 
shape of plot is the Latin square doublet: 


A B 
B A 


A comparison of the variability per plot of such doublets with the vari- 
ability of designs for testing large numbers of varieties can be used 
to assess their inherent efficiency. But we must not expect to obtain 
efficiencies of 100 per cent, since there will always be some loss of 
efficiency with a large number of varieties owing to the number of 
different comparisons on which information is required. Thus if doublets 
of the above type are used to compare v varieties in balanced pairs the 
efficiency factor, if inter-doublet information is not utilized, is v/2(v — 1) 
or, if v is large, 3. Other designs such as lattice designs in blocks and 
lattice squares may be expected to remove fertility differences some- 
what less effectively than Latin-square doublets, and will consequently 
have a variance per plot which is somewhat greater than that for a 
doublet, but they will have higher efficiency factors. The efficiency 
factor of a set of 7 X 7 lattice squares (49 varieties) for example, is 
#. If the variance per plot is the same as for a doublet the inherent 
efficiency is 75 per cent plus whatever information is derivable from 
inter-row and inter-column comparisons. This latter information in its 
turn cannot exceed 25 per cent, and if there are marked fertility irregu- 
larities which are eliminated by the rows and columns, will be very much 
less. 

I have not carried out any examination of the inherent efficiencies 
of lattice and lattice square designs, but I would hazard the guess that 
they will be found on the average to be not very far removed from 100 
per cent. If this is so it is idle to look for more precise types of design. 

There remains size and shape of plot. Up to a point the smaller the 
plots the greater the amount of information per unit area. Also up toa 
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point long narrow plots are more efficient than square ones of the same 
area. Limits are set in both cases, however, by the need for rejecting 
edge rows to avoid inter-plot competition; and even if edge rows do not 
have to be rejected inter-plant competition may result in the amount 
of information being maximum at some plot size well above the minimum 
of a single plant. If the greater part of the costs are proportional to the 
number of plots and not to total area maximum efficiency will be obtain- 
ed with considerably larger plots. 

One way of reducing the variability is to choose particularly uniform 
land for trials. This is good practice—indeed essential—as far as the 
central station is concerned, but the necessity of testing and to a certain 
extent selecting varieties under conditions similar to those in which they 
will be grown limits possibilities in this direction. 

The other major advance of experimental design referred to above, 
namely factorial design, has not been nearly as fully exploited as it 
should have been in plant selection. Factorial design is, of course, not 
of direct use, since the different varieties or strains cannot constitute 
more than a single factor. Far more rapid agricultural progress can, 
however, be made by including other factors such as fertilizers, time of 
planting and so forth in varietal trials. This is of value in the plant 
selection work itself—thus, varieties of cereals are urgently required 
which will stand up to heavy nitrogenous manuring without lodging, 
and the respective merits of different varieties in this respect can only 
be directly tested by a factorial design involving both varieties and levels 
of nitrogenous manuring. It is also of value in that it leads to indirect 
economies by providing information on the other factors, thus rendering 
independent trials on these factors unnecessary. Factorial designs of 
the conventional type are likely to be of greatest use at the later stages 
of the testing programme. At the earlier stages, when large numbers of 
lines or varieties are involved, the other factors can well be applied to 
whole blocks, e.g. those of incomplete block designs. 

Having evolved an efficient testing mechanism we must concern 
ourselves with the ways in which it can be applied in practice. It is this 
aspect of the matter which, it seems to me, has as yet been given very 
inadequate consideration, and it is here that we expect to see large 


increases in efficiency. 


BALANCE BETWEEN THE DIFFERENT STAGES OF THE SELECTION PROCESS: 
NUMBER OF STAGES, PROPORTION OF VARIETIES OR LINES RETAINED AT EACH 
STAGE, ACCURACY AT EACH STAGE 

The most efficient routine for plant breeding and selection will depend 
very much on the genetical situation, and consequently the optimal 
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methods for any given crop will have to be evolved by close cooperation 
between mathematical statisticians and geneticists engaged in the breed- 
ing of that crop. The maintenance of breeding stocks also presents 
somewhat different problems from the production of new commercial 
varieties. In the former the preservation of genetic potentialities is 
required, as well as progress in the desired direction. In a new com- 
mercial variety suitability to the particular needs of the moment is the 
governing factor. The essence of the situation in all crops, however, 
is that the production of new untested strains is a relatively simple job. 
It is the testing that involves the work. It may well, therefore, pay to 
produce a large number of lines, carrying out the initial tests with low 
accuracy and gradually refining the selection. This is what is in fact 
done intuitively by all plant breeders. 

Even with a single stage of selection at each generation the greater 
the number of lines, the greater the chance that some particularly good 
line is included. On the other hand, comparisons between the individual 
lines will be less accurate when a larger number are included, owing to 
the smaller number of plots that can be devoted to each line, and owing 
to the increasing error per plot as the number of varieties is increased. 

In 1939 at the 7th International Congress of Genetics, I pointed out 
(Yates, 1940) that neglecting the latter factor entirely, and assuming 
that the total number of plots is fixed, the average genetic advance due 
to the selection of the apparently best variety, instead of a random 
variety, where there is no retrogression in the absence of selection, will be 


Grae 
GaLN 


where n is the number of varieties, 


G is the genetic variance (distribution assumed normal), 

An is the experimental error variance, \ being independent of n, 

x, is the mean value of the greatest deviate of a sample of n from a 
normal population with unit standard deviation. (Tabulated in 
Fisher & Yates, Statistical Tables. 1938.) 


With A = 74G@ the optimum number of varieties will be 13, in 
which case the genetic variance will be somewhat less than the experi- 
mental error variance. With \ = ya 9@ the optimum number will 
be somewhat greater than 50, and the genetic variance will be about 
twice the experimental error variance. In terms of the ordinary analysis 
of variance the variance ratios between varieties and error will average 
about 1.8 and about 3 respectively. The former is less than the value 
required to give significance at the 5 per cent point. 
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These simple considerations serve to emphasize the value of testing 
a large number of varieties with moderate accuracy instead of only a few 
with very high accuracy. In any series of trials involving only a few 
varieties which give varietal differences that are large compared with 
their standard errors the question should always be asked: would not 
the work have been improved if the same experimental resources had 
been devoted to the comparison of a larger number of varieties? 

The above thoughts were prompted by the habit which had become 
common in certain quarters of judging the efficiency of a variety trial by 
the degree of significance obtained. I had hoped that others would 
continue the investigation, but so far as I know nothing further has 
been done. 

The number of stages introduced into the testing process at each 
generation presents similar problems. If we have a thousand new lines, 
for example, we can, if we wish, test all lines simultaneously, or we can 
test them in a series of stages. Thus we might use three stages, retaining 
yo at each of the first two stages and ending with 10 lines between 
which reasonably accurate comparisons are available; or we might use 
two stages, retaiming 75 at the first stage and ending with 25 lines; 
or we might use two stages, retaining 75 at the first stage and ending 
with 100 lines. Frankly, I have no idea of the relative merits of these 
various procedures, but I suspect that they may be very different. 

Testing in stages is analogous to sequential analysis in sampling 
procedure. One disadvantage, which must not be forgotten, is that it 
lengthens the testing process. This may be serious in certain circum- 
stances in agriculture, since each test normally occupies a season. 

The possibility of testing the suitability of parental lines by testing 
their progenies in a number of different crosses also introduces further 
interesting statistical problems which are in some ways analogous to 
factorial design. If in a cross-fertilized plant, for instance, we test 
all reciprocal crosses between 20 parents (380 lines) we shall obtain 
very accurate measures of the average merit of a parent, even though the 
tests of the individual lines are of low precision. The analysis of this 
type of data has already been discussed (Yates, 1947). A further problem 
that still requires consideration is what weight to give to the average 
parental scores, and what to the scores of the individual lines, when 


arranging the lines in order of merit. 


STANDARD VARIETIES 


In the introduction of new varieties of an established agricultural 
crop we are usually interested in comparing them with established 


varieties. These established varieties can often well serve as controls in 
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the trials. Control varieties provide a standard for measuring the relative 
merits of groups of lines produced in different years. I would emphasize, 
however, that if control varieties are used, three or more should be used 
and not, as is often the case, a single variety. Any one variety may have 
peculiarities such as an abnormal reaction to certain meteorological 
conditions, particular susceptibility to certain diseases, etc. Again in 
vegetatively propagated plants such as potatoes a standard variety 
may become infected by a virus and gradually deteriorate for this 
reason. Consequently, standard varieties, when they are used, should 
be gradually supplanted by the superior new varieties, which in turn 
become the standards. 

If the field trial testing scheme is not to become unmanageable 
varieties should not be permitted to enter the field trial stage until 
they have passed the necessary preliminary tests, but it is equally im- 
portant that potential new varieties should be introduced into the field 
trials at an early stage. This is a defect of the present system of testing 
varieties in the United Kingdom, where the official tests of the National 
Institute of Agricultural Botany always tend to deal with varieties which 
have, in fact, been in current use by the more progressive farmers for a 
number of years on the recommendations of their seed merchants, 
recommendations which are frequently, I suspect, based on somewhat 
inadequate evidence or on trials in other countries from which the new 
varieties have been imported. 

There must also be some orderly method of eliminating varieties from 
field trials, and, I trust, from general agricultural use, when they have 
got to the stage of being superseded by newer and better varieties. 


QUALITY FACTORS 


Testing for quality factors is of the utmost importance. It presents 
considerable practical difficulties, and has in consequence tended to be 
neglected. One trouble is that a fair bulk of produce is often required 
for quality testing. This may require the use of larger plots than would 
be the case if only yield was of importance. Certain quality factors such 
as flavour can only be tested by subjective means. This may have to be 
done by studying an order of preference, i.e. by ranking. C. I. Bliss has 
recently been doing some interesting work on these lines, using in- 
complete blocks. 


ENVIRONMENTAL FACTORS 


. In order to test the suitability of different varieties to different 
soils and climatic conditions, tests must be carried out at different 
centres. Consequently we must organise a chain of experimental sta- 
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tions, or alternatively a series of trials on widely dispersed commercial 
farms. Agronomic factors such as amounts of fertilizers, sowing date 
etc. can be, and should be, dealt with in the variety trials themselves 
by the methods of factorial design mentioned above. For certain factors 
such as lodging, disease response, etc. it may be possible to develop 
special laboratory or greenhouse tests. Such tests are of the greatest 
value as they lead to sharpened and more speedy criteria of selection, 
and enable much more rapid progress to be made. They must, however, 
be confirmed by correlation with field trials. 


CONDITIONS UNDER WHICH SELECTION IS MADE 


The conditions under which the earlier stages of selection are made 
are, I suspect, sometimes at fault. To take a specific example, the 
stiffness and shortness of straw are very important properties in cereals. 
Yet, so far as I know, the first stage of selection of new lines is generally 
made from plants sown in close proximity and judged on their individual 
merits, particularly by apparent vigour. Thus, there is at this stage, 
as there is in a wild population, considerable selective advantage in 
height. This seems to be a situation in which selection under deliberately 
unnatural conditions at the early stages might well result in considerable 
advances, but I leave this to the geneticists. 
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RECENT APPLICATIONS OF 
BIOMETRICAL METHODS IN GENETICS 


(2) THE ANALYSIS OF SELECTION CURVES 


Lurer L. Cava 
Dept. of Genetics, The University, Cambridge 


Qyeamenss DATA on spontaneous and artificial selection exist in 
the genetical literature, but no special attempts appear to have been 
made to produce statistical methods for their analysis. In this paper a 
method will be given for the analysis of a certain type of selection curves, 
namely those arising when two alleles in competition come to an equilib- 
rium after a sufficient number of generations, the frequency of each allele 
at equilibrium being independent of the initial frequencies. ‘The method 
can be easily extended to cover other types of selection processes. 

The usual genetical hypothesis in this case is that the heterozygous 
genotype Aa, is at an advantage over both homozygotes AA, aa. The 
equilibrium will then depend on the relative survival values of the two 
homozygotes (Fisher, 1922, Haldane, 1926). 


THE SELECTIVE PROCESS 
Let p and g = 1 — p be the gene frequencies of alleles A and a. 
Let the relative survival values of the three genotypes be 


genotype AA Aa aa 
survival value a B Y 
frequencies =p Png 


The last row gives the expected frequencies of the three genotypes if 
mating is at random. 

If generations do not overlap, then the change in p in the course of 
one generation is 
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This difference equation is not well adapted for mathematical treatment. 

And therefore, even in the case where generations are distinguishable, 

it seems convenient to replace it by a differential equation. In cases 

where the ratios a/8, y/8 do not greatly differ from unity, the differential 

equation will represent the course of the selection process to a good degree 
of approximation. 

There is however an additional reason, apart from mathematical 
convenience, which makes it preferable to represent the process by means 
of a differential equation. In Nature, as well as in most experimental 
arrangements, generations are often not separated in time; on the con- 
trary, they may overlap to such an extent, that no distinction is possible, 
and even the calculation of an average generation time is very difficult. 
The selective process is therefore best regarded as being a continuous 
one, and hence is most naturally represented by means of a differential 
equation. 

We shall therefore set up a differential equation which represents a 
certain continuous process. It is then assumed that the actual process 
of selection when the generations overlap approximates to this continuous 
one. 

The form of the differential equation is suggested by equation (1) 
and is obtained by various formal transformations of (1). Let dt be an 
element of time. Then in equation (1) substitute 


a — 
= A dt 
(2) 


ms 
abo COL 
B 


and set dp in place of Ap. The numbers A and C are essentially dif- 
ferences of Malthusian parameters. 
We obtain 


peep ier Adi ey pea 
= VA de Cai 


_ —pylCg — Ap) dt __ 
1 — Ap’ dt — Cg dt 
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which, on rejecting powers of dt higher than the first, becomes 

(P _ na(Cq — Ap) (3) 
dt j 

On integration of (3) we obtain 


i eters: 
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At equilibrium dp = 0, and the equilibrium gene frequencies pz and qz 
will be: 


| = i+ const. (4) 


jae “asa eae 
pa {re Gin = eae oy CT: 


We can therefore write equation (4) as follows: 


as log p + i log q — oid log | pz — p | = (A + C)(t + const.) (5) 
Dr dp Prde 


Equation (5) was suggested by Professor Fisher as providing a transfor- 
mation of gene frequencies which would give a function linear with time, 
if the equilibrium frequencies pz and gz are known. 

Let us call 


1 
Yi esa aia = ] = 6 
Py dea Pak caer 2) (6) 
the transformed gene frequencies, and 
Y=a-+bOt (7) 


the linear function connecting Y and time of selection. The slope of 
the linear function, b, which is equal to (A + C), in conjunction with the 
equilibrium frequency pz will allow the calculation of the survival values 
AandC. The position of the straight line (7) given by a, will depend on 
the initial gene frequency, being equal to Y, (for time ¢ = 0). The slope 
is independent of the initial conditions; therefore the problem of equality 
of A + C in independent experiments, started with different initial 
frequencies of the same genes, will be a problem of parallelism of the 
corresponding straight lines after transformation (6). The goodness of 
fit of the straight line to the transformed gene frequencies will provide a 
test of the hypothesis that the selective process can actually be described, 
in a satisfactory way, by means of the differential equation given above. 


THE PROCESS OF FITTING 


The problem of estimating the parameters A and C, and the other 
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related problems mentioned before, are essentially centered on the process 
of fitting the straight line (7) to experimental data. Table I gives values 
of Y as a function of p and py. When pz is exactly known, the observed 
gene frequencies will be transformed into Y by interpolation of the table, 
and the Y’s will be plotted against time. Best estimation of a and b in 
equation (7) involves, however, weighting of the Y values and, as the 
weight will depend on the expected values of Y, only a process of 
successive approximation will in practice give best estimates for a and b. 


TABLE I 
TABLE OF THE VALUES OF THE FUNCTION: 
1 1 1 
tee 10S e.g — 46 = 
Gee Tagen ica, PAD =p) 


Nez 0.50 | 0.45 | 0.40 | 0.35 | 0.30 | 0.25 0.20 0.15 0.10 0.05 


PN 


0.01 | —6.38| —6.93| —7.61] —8.43] —9.47] —10.82] —12.66] —15.29) —19.31] —24.35 
0.05 | —2.90| —3.05| —3.20] —3.35] —3.46| —3.47) —3.19]| —1.97] +3.27 +25.5210.95 
0.10 | —1.15| —1.07| —0.92| —0.65| —0.16] +0.77| +2.75) +-8.02 +16.91) +12.86]0,90 
0.15 | +0.08] +0.35] +0.76] +1.40] +2.48) +4.48] +9.03 +14.13] +10.36] +8.80)0.85 
0.20 | +1.15] +1.62] +2.31] +3.40| +5.28] +9.24 +12.50] +9.24| +7.52) +6.59]0.80 
0.25 | +2.20] +2.90] +3.96} +5.72] +-9.23 +11.43] +8.48] +6.90] +5.85) +5.23]0.75 
0.30 | +3.32| +4.34] +5.99] +9.18 +10.69| +7.93] +6.43] +5.45| +4.73] +4.28]0.70 
0.35 | +4.23] +6.19] +9.14 +10.15| +7.51] +6.07} +5.12| +4.43] +3.90} +3.56]0.65 
0.40 | +6.36] +9.14 +9.76} +7.18] +5.77| +4.84| +4.16] +3.65) +3.24] +2.97/0.60 
0.45 | +9.19 +9.49]+6.91| +5.52) +4.59] +3.92} +3.42} +3.02] +2.69] +2.47/0.55 
+9.30| +6.71| +5.29] +4.36] +3.70] +3.19] +2.80] +2.47| +2.22) +2.04/0.50 
+9,19] +6.52] +5.08] +4.14] +3.47] +2.97] +2.57| +2.26] +2.01) +1.80} +1.66/0.45 
+6.36| +4.86] +3.90] +3.22] +2.72] +2.33} +2.03] +1.78] +1.58) +1.40) +1.28]/0.40 
+4,63| +3.64] +2.95] +2.45) +2.06] +1.76] +1.52] +1.33] +1.17) +1.03} +0.94)0.35 
+3.32| +2.62] +2.12|+1.74| +1.45] +1.23] +1.04] +0.89] +0.77| +-0.67) +-0.60)0.30 
+2.20| +1.70] +1.34] +1.07] +0.86] +0.70} +0.57] +0.46] +0.37) +0.30) +-0.24]0.25 
+1,15] +0.82] +0.58]} ++0.40] +0.26] +0.15] +0.07] +0.00] —0.05) —0.10} —0.13}0.20 
+0.08| —0.11] —0.24] —0.34| —0.41] —0.45] —0.49] —0.52)| —0.54) —0.55) —0.55)0.15 
~1.15| —1.19] —1.21] —1.22) —1.21] —1.19] —1.18} —1.16] —1.13] —1.11) —1.10)0.10 
—2.90| —2.76| —2.63] —2.51] —2.40] —2.30] —2.20] —2.12] —2.04) —1.96) —2.91]0.05 
—6.38| —5.91| —5.50| —5.15| —4.85| —4.57| —4.33] —4.12! —3.92| —3.75| —3,63)/0.01 


0.50) 0.55} 0.60) 0.65 0.70 0.75 0.80 0.85 0.90 0.95 0.99 Ne 
PR 


In the more general case, however, the equilibrium frequencies are 
not known, or at least are not exactly known, and therefore must be 
estimated from the data. The same is true of the initial gene frequency, 
Po, Which will be known only within a sampling error. Therefore, the 
estimation of A and C will be obtained through the estimation of the 
three parameters a, b, of equation (7), and pz, from the data. 

Let us call 6; any of the three parameters; choose p so that pz > Dp; 
call n the number of observations on which each gene frequency is 
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determined; call p the observed, and P the expected (similarly g = | — p 
and Q = 1 — P) gene frequencies. The maximum likelihood equations 
are 


00; PQ 00; . 

there being one such equation for each of the three parameters. No 
direct solution of this system of equations is available, and therefore 
a method of successive approximation will be used, as is usual in other 
applications of maximum likelihood to related problems (see for instance, 
Probit Analysis, Fimney 1947). 

Equations (8) will be expanded according to Taylor-Maclaurin series, 
with consideration of first orders only, giving 


OL OL OL 
F) F a 
aL 00; 00; 00; 
— v A a 7 9 
06; me 06, sek. 00, Tages 005 e 


Trial values of the three parameters will be inserted in equations (9) 
and from this, the adjustments for each parameter, 60; , will be calcu- 
lated. The whole process can be repeated, using the adjusted values of 
the three parameters, and such cycles of computations will be repeated 
until no further adjustment is necessary (the adjustments being small in 
comparison of their standard error). All computations are greatly sim- 
plified by taking, in the second derivatives of L, p = P; in this case the 
system of equations will become, writing it in extenso, 


ips) OP ws ee _n_ oP dP 
2 ae OO ye » BG 00, meses PQ 06, 065 
n OP aP 
+ 9% 24 BO a6, a6, 
Oe _n OP oP (ey 
dX PQ 06. — obs Dd BQ 00, 065 002 » PG 060, 
(10) 
m Qe Oe 
+ 56 2) 36 38, a6, 
np = P) oP We OPSeP n OP dP 
eS es eee 50, pale Re tena ste eno Sa a 
» POs 086; » PQ 00, 005 2 PO 08 0b: 


ap \? 
50, re (S) 
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To obtain the three derivatives 0P/d6; which, in the present instance, 

are 0P/dpz , 0P/da, AP/db, let us consider the function 
¥ = gz log P + pz log Q — log (pe — P) — psqgsY = 0 (11) 


We shall then have 


Beer 
and similarly 
Sp = POPs — Pit, 5 = — POs — Pe 
where z stands for 
c= —-1 | tog gp - 1, - (os - pa | 
Peqe pe 


Inserting these values into equations (10), and setting 


op pg) 


i= POlpy— P) 
we shall obtain the following equations: 


dons = 6a >i nw + 6b di nwt + dpe D> nwz 
> nst = 5a >. nwt + bb D> nwt? + Spe D> nwtz (12) 


> nsz = ba D> nwz + 6b Do nwte + ope Dy nee? 


This method can be considered as a simple extension to the case of more 
than one parameter of the method of scoring, proposed for the analysis 
of linkage data by R. A. Fisher (1947), and which can be considered as 
the most direct form of solution of maximum likelihood equations by 
successive approximation. The quantities on the left in eq. (12) are 
then scores S,, in respect of the three parameters a, b, pz ; the matrix of 
the coefficients of the adjustments is the matrix of information J, in 
respect of the estimates of the three parameters, or the variance matrix 
in respect of the scores. System (12), in matrix notation will be 


S; = al 


I being the information matrix and 6 the vector of the adjustments 


da, 6b, 6pz . The adjustments are obtained from 
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6 == ‘Sale 


5 . < o =} 
and their variances from the elements of the principal diagonal in I. 
APPLICATION 


In an experiment of artificial selection on Drosophila pseudoobscura, 
Dobzhansky established populations with various initial gene frequencies 
of certain chromosome arrangements, which, for our purposes, can be 
considered as single genes, and observed that an equilibrium was reached 


80 
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after a certain number of generations. In one such experiment started 
with about one-third of chromosomes of the Standard type, and the rest 
of the Chiricahua type, and finished after nine months when about three 
fourths of the chromosomes were of the standard type, samples of about 
150 individuals were taken and their chromosomes examined at various 
times after the mixture of the two original types. Frequencies observed 
at given times (time units of 10 days) are given in the second column of 
table II. Numbers of chromosomes counted are given as n in the third 
column. The same data are shown graphically in fig. 1. 
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FIGURE II 


The equilibrium frequency is not known; from an inspection of fig. 1, 
assuming that equilibrium had not yet been fully reached at time when 
observation ceased, a first trial value of pz was taken at pz = 0.75, 
and hence values of the transformation (6) corresponding to the observed 
gene frequencies were interpolated from table I and plotted in fig. 2 
against time. These values are indicated as y in the fourth column of 
table II. 

In fig. 2, a straight line was fitted by eye to the data, being 


Y = 2.13 + 0.48¢ 


First trial values will be indicated with the suffix 1. They are 


ja 2-13 
b= 048 
pr, = 0.75 


Expected values, calculated from the straight line fitted graphically, are 
indicated with Y in the fifth column of table II, and from Y, expected 
values P are interpolated by means of table I (or an extension of it). 
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TABLE II 


FIRST CYCLE OF COMPUTATIONS ON DOBZHANSKY’S DATA 
t p n y y ‘B. nw ns Zz nwz nut 
0} .383) 1278 | 2.1 | 2.13] .3886 |40.138 |—1.3896 |+ 6.56)/263.2528) 0 
3] .530]} 300 |] 4.2 | 3.57] .491 | 5.04 |+8.030 |4+ 10.93) 55.0872)15.12 
6| .633/ 300] 6.8 | 5.01} .569 | 2.41 |+38.475 |+ 17.33} 41.7653)14.46 
10} .603) 300} 5.9 | 6.93) .635 | 0.920 |—1.104 |+ 30.88] 28.4096; 9.20 
13} .653} 300] 7.6 | 8.37] .667 | 0.459 |—0.349 |4+ 45.65) 20.9533] 5.967 
16} .653} 300] 7.6} 9.81] .689 | 0.2389 |—0.659 |+ 65.44] 15.6402) 3.824 
20} .704| 250 | 11.5 /11.73] .708 | 0.0912)—0.0420|+100.37) 9.1537) 1.824 
25} .720) 300 | 18.6 |14.13]-.722 | 0.0472)—0.0168/+157.81|} 7.4486) 1.180 
49 .3364|+2.9378 441.7107/51.575 


By means of the formulas given in the second section, from p, P 
and Y and from py = 0.75 the values of s, w and z are calculated, and 
hence the last five columns of table IT are obtained. 

From these, the following other sums of squares and products are 
easily calculated: 


S nwt’ = 428.85; S nwtz = 1591.88; S nwz” = 8004.36 
Snst = +2.559; S nsz = —15.8332 
With such values, the system of the adjustment equations is formed: 
scores information matrix X adjustments 
+ 2.938 = 49.34 6a+ 51.57 6b + 441.71 dp, 
+ 2.559 = 51.57 da + 428.55 6b + 1591.88 Spy 
— 15.833 = 441.71 6a + 1591.88 6b + 8004.36 dp, 
The inverse of the matrix of information is 
+0.0755123 +0.0244006 —0.0090198 
+0.0244006 +0.0167924 —0.0046861 


—0.0090198 —0.0046861 +0.0015546 
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and multiplying this by the matrix of scores, one secures the adjustments 
to the trial values of the three parameters 


6a, = 0.0755123 X 2.938 + 0.0244006 x 2.559 
+ 0.0090198 X 15.833 
= 0.4271 + 0.275 


0.1889 + 0.130 


6b; 


épe, = —0.06311 + 0.0394 


The standard errors of the adjustments are the square roots of the 
corresponding elements of the principal diagonal in the inverted matrix; 
thus, the standard error of 6a is ~/0.0755123 = 0.275. The adjusted 
values of the three parameters, after this first cycle, will then be 


d, = 2.56, b= 0.67, pp, = 0.687 


The adjustments are greater than their standard errors, and therefore 
it seems worthwhile repeating the computations, taking as estimates of 
the three parameters the adjusted values given above. Having per- 
formed this second cycle, the following adjustments were found: 


da, = —0.0222 + .255 


6b, = +0.0557 + .167 


I 


pe = 40.0122 4.023 


The adjustments are now well below their standard error; in normal 
cases, therefore, no further refinement would be necessary. In this case, 
however, a third cycle was carried out, and for reasons which will be 
apparent later, the inverse of the matrix of information is given in full 


below. 


Third cycle. Inverse of the matrix of information Scores 
+0.0592353 +0.0178369 —0.0035900 —0.3908 
+0.0178369 +0.0320025 —0.0031734 +1.1092 
—0.0035900 —0.0031734 +0.0004880 — 0.0825 


The solutions for the adjustments are: 
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da; = —0.00307 + .243 
6b; = +0.02879 + .179 
dpe, = —0.00216 + .022 


and the final estimates of the three parameters are: 


a,= 2.54, b= 0.76, pp, = 0.697 


In the final stage, the goodness of fit can be tested in the usual way by 
means of x’, comparing expected and observed gene frequencies. This 
has been here done at all stages, and the results are given below. x° was 
calculated with three figure accuracy only for the last fit. There are, in 
all, eight observations, and three parameters have been estimated from 
the data, so there remain five degrees of freedom. The probability 
corresponding to each x’ is given in the bottom line. 


Analytical fitting 
Graphical fitting 


after 1st cycle 


after 2nd cycle 


after 3rd cycle 


10.32 9.43 Th AMR 7.036 
<<a) <.10 <.30 <.30 
> .05 > .05 > .20 > .20 


One further point is the calculation of the values A and C, which are 
the ones that really matter from the genetical point of view. Since 


Pm Cc. do A = 
we shall have 
A= baz, C= bps 


To calculate the standard errors of A and C, an approximate formula is 
available, and we shall need some of the values contained in the inverted 
information matrix. Calling V, , V,, and W,,, the elements of the 
inverse information matrix of places ¢2. , ¢33 , C23 , We shall have: 


V(A) Ga Vs a b: Vee a= 2baeWoor 


V(C) = DrV> Sie bie + 2bpEWopr 
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Taking these values from the matrix obtained in the last cycle of compu- 
tations (where V, = 0.0320025; V,, = 0.0004880; W.,, = —0.0031734) 
we obtain numerically: 


A = 0.230 + 0.063 


C 


The values of A and C are expressed in the time scale used in the calcu- 
lations, which is here a time unit of 10 days approximately. The solid 
line in fig. I represents the theoretical curve they calculated. 


0.530 + 0.119 


DISCUSSION 


The method here developed allows the estimation of the parameters 
involved in a selection process, in which the heterozygous genotype has 
an advantage over both homozygotes, and therefore the frequency of the 
two alleles in competition comes to an equilibrium value, which depends 
on the relative survival values of the homozygous genotypes. When esti- 
mates of these parameters have been obtained, tests of goodness of 
fit become possible. It is however to be noticed that one necessary con- 
dition for the application, in the present problem, of the method of 
maximum likelihood and of the usual tests of goodness of fit, may not be 
fulfilled, namely the condition of the reciprocal independence of the 
observations. A selective process is a stochastic process, and a gene 
frequency at a given time depends on gene frequencies at earlier genera- 
tions unless the population is infinite. While such difficulty need not 
worry us for most populations of animals and plants in Nature, it seems 
possible that it may need to be taken into consideration for artificial 
populations, which are usually of small size. The total effect will be of 
a greater variance than that due to sampling only, i.e. a bad fit—which 
might however arise through a multitude of other causes. In the example 
considered above, the deviation from the theoretical course of selection 
is not significant, and therefore there is no detectable evidence of drift 
or other cause of departure. Where drift is suspected, the only method of 
treatment so far known is that which has been proposed by Fisher and 
Ford for the analysis of a natural population of Panaxia dominula. 

One consideration which applies to drift, however, is that when a 
small sample of the population is observed each time, the sampling 
variance will be many times greater than the variance of gene frequencies 
due to the finite size of the population. The smaller the ratio of the size 
of the sample to the size of the total population, the better the approxima- 
tion to the condition of independence necessary for an unbiased applica- 
tion of the maximum likelihood, and of tests of goodness of fit. In the 
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example of which use has been made here, the size of the sample was 
300, and the size of the population twenty times larger, a ratio that 
seems in fact to provide the conditions for the application of the test. 

This investigation was started with the aim of analyzing experimental 
data by Prof. A. Buzzati-Traverso; the analysis of such data will, how- 
ever, appear as an appendix to Buzzati’s paper. The whole investiga- 
tion would not have been possible without the generous advice and help 
of Prof. R. A. Fisher. My gratitude is also due to Dr. A. R. G. Owen 
for revision of the manuscript. 
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RECENT APPLICATIONS OF BIOMETRICAL 
METHODS IN GENETICS 


(3) SCORES FOR THE ESTIMATION OF PARAMETERS 


D. J. FINNEY 


Lecturer in the Design and Analysis of Scientific Experiment, 
University of Oxford 


HE BIOMETRIC STATISTICIAN should not rest content with developing 

sound methods for the statistical analysis of biological data: he should 
also put these methods into the form in which they can most easily be 
applied by a biologist with no extensive knowledge of mathematics. 
He will thereby frequently ensure that the labour of statistical analysis 
forms only a small part of the total effort in a particular research project. 
Experimenters often argue their preference for a non-efficient method of 
estimation or analysis on the grounds that it is more easily understood 
and performed. This may be a false economy, for a rearrangement of 
the calculations and the provision of suitable tables may much simplify 
the work and make fully efficient, or at least very highly efficient, 
methods available with the minimum of labour. 

In no branch of biology is the need for efficient statistical methods 
more evident than in genetics. The collection of genetical data may 
entail the breeding and recording of plants or animals for many years, 
or the patient search for informative human pedigrees. Yet records 
collected at great cost are often subjected to methods of analysis, 
chosen without regard to statistical efficiency, methods which may fail 
to extract as much as 50% of the information available. 

A very convenient computing scheme for the estimation of genetic 
parameters, such as gene frequencies and recombination frequencies, is a 
scoring system which assigns a numerical score to each individual or 
each family in the records and bases the estimate upon an average score 
for all the data. Such a method can be very quickly applied, especially 
when tables of scores have been prepared. 

Suppose that @ is a single unknown parameter to be estimated, say 
the frequency of recombinations between two loci. If we have data 
covering, say, two generations, all possible families may be classified 
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according to patterns, such that the relative frequencies of different 
patterns are independent of 6. For example, all records of two parents 
of specified phenotypes with a particular number of offspring would 
constitute one pattern, all records of a particular number of sibs for 
which no parents have been classified would constitute another, and so 
on. Within any one pattern, the probability of obtaining for the records 
individuals of particular phenotypes would be a function of 6, say P(@): 
using S( ) for summation over all families of the same pattern 


S(P) = 1. 


Take 6 as an approximate estimate of 6, obtained by any rapid 
method (by inspection, or by rough analysis of all or part of the data) 
and write 


6= %+ 4, 


where 6 is an adjustment to 6. Then to the first order in 6 


P(A) = PCH) + 8 


6=6o 


== IP ae IRGC say 
and 


S(P,) = 1, 


S(Po) = 0. 


Let z be a score calculated from the observed phenotype frequencies 
in a family: zg may be any function of the several frequencies, subject 
to the obvious restriction that its average value for all families of the 
same pattern must be dependent on @. Now 


E(z) = S(P2) 


= S(Poz) + dS(P%). 


Provided that S(P{z) is not zero, an average value of z from observations 
can be used in association with S(Poz) and S(P{z) to give an estimate 
of 6: this is then added to 4 to give a revised estimate of 6. For con- 
venience in handling, a score which leads directly to the revised estimate 
has advantages: write 


~ _ Sz) z 
Ln (t ee + ge)’ 
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then 
Ey) = 6 + 6 
= f. 


Consequently, simple averaging of y over all families of the same pattern 
gives an estimate of @ directly. Moreover, ignoring terms in 6, 


_ SP) = {SPA} 
(S(P2)}? 


Hence, if for any one pattern of family we write 
W = 1/Vqy), 
and then combine the scores from all the records into the weighted mean 


Wy 


saris 


(where >> denotes summation over all data), y will be the most precise 
average value of y that can be formed and will have expectation 0. 
Therefore, y may be taken as 6, , a revised estimate of 06. 

If y and W happen to be independent of 6 , this is the end of the 
calculation: y is the best estimate of @ that can be based upon the 
scores y. More usually, however, both y and W will be functions of 
6) , and iterative scoring must be practised. Values of y and W are 
obtained for an initial 4, ; these give y, which is taken as 6, ; new sets 
of values of y and W are then based on 6, , and the new y is taken as 6, . 
This continues until two successive 6’s are sufficiently close together 
for, their difference to be ignored. The latest value is then taken as 
the'estimate, 0, , and 


Vy) 


<I 


1 
PE AO ite green 
ea 
The procedure sounds laborious, but in fact is simple and rapid if tables 
of y and W as functions of # have been prepared for different patterns of 
family. 

The estimate reached depends upon the form of z adopted in each 
pattern of family. Indeed, the method is a general expression of many 
existing methods for the estimation of genetic parameters, which consist 
merely of equating some arbitrary function of phenotype frequencies to 
its expectation. A criterion is now needed for deciding between rival 
scoring systems. As is well known, an asymptotically efficient estimate 
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of 6, that is to say an estimate which in large samples has the least 
possible variance, may be obtained by the method of maximum likeli- 
hood; this requires that 


> (log P) 


be maximized for the data. This may be accomplished by scoring with 


d 
L = 19 (08 P) Z 


90 
Po 


The corresponding y score 1s 


with 


W 


lI 
= 
SI 

nw 


The weight W is a well-known expression for the total information 
provided by a single family of particular pattern if 6 = 6) , thus demon- 
strating the efficiency of the estimate based on Y. It is convenient to 
tabulate Y and 


P! 
A= WY = W — 
Won P, 
for different values of 6 , so that the revised estimate may easily be 
obtained as 


eae 
POSSE 

Often A and W are both linear functions of the phenotype frequencies 
amongst the progeny of a family, so that for specified parental types all 
family sizes may be covered by tabulation of weight and score per indi- 
vidual. The iteration goes as before. 

The efficiency of any other z-score may be assessed by comparing 
the information extracted by it with the total amount available: 


{ S(P%z)}? 


Le) : 
(FE \isee) Bae oh ce) a 
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This formula assists a decision as to whether a score other than the 
maximum likelihood Y involves a sacrifice of information unreasonably 
great in comparison with the saving of labour. It is sometimes possible 
to find a scoring system different from Y and more easily calculated 
which is nevertheless of full efficiency, but more usually other scores 
will sacrifice at least a little of the available information. 

Similar methods can be applied when two or more parameters must be 
estimated simultaneously, though the calculations are then necessarily 
more complicated. The maximum likelihood estimation of the param- 
eters of a normal tolerance distribution by means of the probit transfor- 
mation is an example outside the field of genetics (Finney, 1947; 1949a). 

Tables for certain applications of the method to gene frequency 
and recombination frequency estimation have been given t elsewhere 
(Finney, 1948a, b; 1949b), together with illustrations of their use. 
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(These papers contain bibliographies more complete than can be given here and show 

how the ideas were first developed by R. A. Fisher and others.) 


DISCUSSION ON RECENT APPLICATIONS OF 
BIOMETRICAL METHODS IN GENETICS 


(1) EXPERIMENTAL TECHNIQUES IN PLANT 
IMPROVEMENT 


F, Yates 


W. G. Cochran. Some time ago I presented the ideas in Dr. Yates’ paper 
tO a group of plant breeders in the United States. Using a number of plots that 
appeared typical and typical estimated values for the genetic and environmental 
variances, I showed that the optimum number of varieties would allow only one or 
two replications and that the F value had only a small chance of being significant. 
The audience was somewhat shocked, apparently because my results seemed to con- 
tradict the usual advice of the statistician that several replications must be used and 
that F tests must be made. This suggests that rather patient explanation will be 
needed in order to make Dr. Yates’ results clear and acceptable to plant breeders. 

With regard to selection in several stages, some unpublished work on two-stage 
selection has been conducted in the United States. The work suggests, as a specula- 
tion, that the optimum is fairly flat, i.e. that a considerable deviation from the 
optimum procedure does not greatly decrease the genetic advance. 

Rk. A, Fisher. The selectionist requires either significant, or at least subsignificant 
evidence of genetic variability at the first stage of quantitative testing. His program 
for the next 3-5. years may be decided on the strength of this evidence. So I can 
understand that American plant breeders should be a little shocked at the suggestion 
that the z test should be ignored at the first stage of selection. 

G. Pompilj. Je voudrai développer quelques considérations critiques tendant & 
éclaircir la signification logique du probléme du “‘testing”’. 

Je crois qu’on doit distinguer ce qu’on voudrait obtenir et ce qu’on peut obtenir. 
Ce qu’on voudrait obtenir est trés clair. Nous voudrions obtenir des procédés, tout 4 
fait objectifs et impersonnels, pour classifier, selon certains caractéres, les différentes 
lignes ou variétés qu’on a obtenues. 

Mais les données & notre disposition sont elles en général suffisantes pour ce but? 
Je crois que non. En effet le probléme que nous sommes posé, comme monsieur Gini 
la observé il y a plusieurs années, demande, pour étre résolu, deux ordres de con- 
naissances; & savoir: la distribution a priori de toutes les variétés en tenant compte 
des caractéres prefissées et la distribution des mémes variétés apres les expériences. 
Or, en général, nous ne connaissons pas la distribution a priori; il nous manque done 
un des deux éléments nécessaires pour resoudre le probléme. 
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(2) THE ANALYSIS OF SELECTION CURVES 
Lurer L. Cava 


J. B. S. Haldane. It is extremely gratifying that Dr. Cavalli has obtained such 
good agreement with a simple hypothesis. Insofar as there were no deviations from 
linearity, the data are consistent with the following hypotheses. 

1. Mating was at random, the influence of inbreeding and assortative mating 
being negligible. 

2. Selection was Darwinean, that is to say it proceeded as if the three genotypes 
survived to maturity in different proportions, so that only two parameters were 
needed. Often selection may depend on the interaction of two genotypes, as when a 
foetus heterozygous for Rh sensitivity has a decidedly different expectation of life 
according as its mother is or is not Rh-sensitive. If the q different matings, e.g 
AAY X AAG, Aa X aad have different fertilities, we should reach a five para- 
meter equation. 

3. Selection was of constant intensity, independent of the composition of the 
population, although this might, for example, have altered the culture medium. 

L. L. Cavalli. I entirely agree with Professor Haldane about the simplicity of the 
hypothesis considered. However, I believe that as a general policy introduction of 
new parameters may be dispensed with until goodness of fit with old parameters is 
satisfactory, and as far, of course, as old parameters make sense from a biological 
point of view. 


(3) SCORES FOR THE ESTIMATION OF PARAMETERS 
D. J. FINNEY 


M. J. R. Healy. When more than one parameter is to be estimated, several cycles 
of iteration are more frequently required. In this case the process leads to simul- 
taneous equations. During the iteration, it will save computation to adjust only the 
right hand side of these equations, exact values of the left hand side being only 
required at the final stage to obtain variances and covariances. 

G. Pompilj. Je voudrai donner un avertissement qui se rattache A ce que j’ai déja 
dit 4 propos de la premiére communication. II s’agit d’une faute de logique qui, par- 
tant de la statistique, va infecter les autres branches de la science et méme de la 
biométrie. 

Je parle des théories de signification et d’évaluation des paramétres, dont 
Vinconsistence logique a été démonstrée par monsieur Gini il y 9 dix ans. La chose 
est trop longue 4 développer ici. Je me limiterai done & vous recommander, si ga 
n’a pas trop de présomption, de vous méfier toujours de tous les “tests” de significa- 
tion et de toutes les théories qui promettent d’atteindre a Vévaluation des paramétres, 
car on peut démontrer que pour obtenir des résultats corrects il faut connaitre des 
éléments qui en général ne sont pas & notre connaissance. 

Other participants in the discussion included F. Chodat and F, Bernstein. 


INDUSTRIAL APPLICATIONS OF BIOMETRY 


BIOMETRIC METHODS IN THE CHEMICAL INDUSTRY 


O. L. Davies 


Biological Laboratories, Imperial Chemical Industries Ltd., 
Hexagon House, Blackley, Manchester 


Abstract 


N THE CHEMICAL INDUSTRY a great deal of research is directed to the 
I improvement of yield and quality of chemical products. This research 
usually consists of assessing experimentally either in the laboratory or on 
the plant, the effect of various changes in the process conditions, for 
instance, changes in the temperature and time of reaction, concentration 
and variety of the reactants. Several of these factors may be examined in 
the same experiment, and factorial arrangements lend themselves par- 
ticularly well to this kind of research. In these experiments the observa- 
tions are made sequentially, that is, one after another, and this gives 
a large measure of flexibility to the experimental design, because it is 
possible to modify the design during the course of the experiment in the 
light of the results obtained. Partial factorial designs in sequence are 
particularly adaptable to this type of experimental work, and two experi- 
ments relating to the problem of improving the persistence of penicillin 
in the body are described to illustrate these principles. 

Simple methods are introduced for the construction of partial factorial 
designs and these are illustrated by an investigation on the large scale 
manufacture of penicillin, in which a relatively large number of factors 
were investigated. In investigations on the plant scale it is usually 
necessary to allow for time trends and this may be done by using the 
well-known principles of confounding, with “blocks” spaced in time. 
Typical partial factorial designs for the 2" and 3” systems are given. 

In experimental work in industry we frequently know from past 
experience the magnitude of the experimental error. This information 
may be used to assess the significance of the effects in a partial factorial 
experiment when there are insufficient degrees of freedom in the experi- 
ment itself to obtain a reliable estimate of the experimental error. 
Another use is in estimating the size of a projected experiment and the 
considerations involved here are the control of the risks of errors of the 
first and second kind. The levels decided for the risks of these two 
types of errors are dictated by the economic consequence of making 
the errors of the first and second kind. 
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DISCUSSION FOLLOWING 
BIOMETRIC METHODS IN THE CHEMICAL INDUSTRY 


O. L. Daviss 


Besse B. Day. Dr. Davies has furnished us with an invaluable guide to the prac- 
tical use of partial factorials. Industrial statisticians are deeply indebted to bio- 
metricians for many modern statistical techniques which are equally applicable to 
engineering problems but often with a different emphasis. Thus the prime interest 
may be in the error term itseli—how reproducible are the results? 

Among the applications of interest in both biological and industrial fields is a 
war-time statistical development in ordnance testing which may have application to 
research for determining the 50% lethal dose. By this method of sensitivity testing 
the height from which a specimen of explosive is dropped is changed by a fixed incre- 
ment after each specimen has been tested, increased if the specimen did not explode 
and decreased if it did explode. This ensures many results near the 50% point from 
which the mean height for 50% and its standard deviation may be estimated. A tool 
with biological and industrial applications is spectrographic analysis. The evalua- 
tion of its reliability and of conversion factors or formulas are statistical problems 
which require experimental designs, regression analysis, etc. Bacteriological counts 
have an industrial parallel in estimating solid contaminants in used lubricating oils 
from crankcases. In general the problems in industrial experimentation requiring 
statistical tools are (1) instrumentation, more practically calibration, (2) development 
and standardization of test methods, (3) specifications—setting of tolerances, (4) im- 
provement of a product and (5) development of new products. 

E. C. Fieller. Dr. Davies’ paper was to be admired not only for its adaptations 
of modern experimental design to the needs of his industry, but also for the wise 
restraint with which these applications were made, and the care taken to ensure 
beforehand that the conditions under which fractional designs might be validly used 
were in fact justified. 

As an example of the difficulties that might arise without these preliminary as- 
surances, an experiment in the design of radio valves has been described, which 
started off on a highly fractional basis but had eventually to be made a complete 
factorial experiment. 

Many speakers at this and previous sessions have found difficulty in distin- 
guishing between Statistics and Biometry but perhaps the answer is suggested by 
A. GC. Bacharach’s remark that “Every chemical experiment involves some living 
organism, even if it is only the chemist.” 

A. Hald. I should like to stress two points. First, the need for statistical investi- 
gations of the testing technique. Very often the laboratory error is considerably 
larger than that shown by the usual duplicate analyses, because the conditions of the 
laboratory may vary in time and produce a trend in the measurements. This trend 
does not become apparent from duplicate analyses and is often more important than 
the errors of measurement proper. Therefore, it is necessary to keep the laboratory 
errors under statistical control. 
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The second point to be noted is the sequential nature of industrial experiments 
and the corresponding time trend. The fitting of trends with polynomials or with a 
moving average seem tome well suited to take this time factor into account and might 
be considered as useful supplements to the usual analysis of variance. Especially as 
industrial experiments are often very costly, it may pay to use regression analysis 
instead of the analysis of variance. 

H. Astrand. The discussion of Dr. Davies’ paper makes me doubt that ‘‘bio- 
metrics” describes it adequately. Would it not be better to speak of “statistical 
methods”, in spite of the fact that to a large extent they are developed by bio- 
metricians? 

It seems that the field for applying statistical methods in industry is immense. 
I feel sure that many more statisticians from industry will take part in future meet- 
ings of our Society and with that expectation I close the session. 


NEWS AND NOTES 


The Survey Research Center of the University of Michigan will hold 
its Third Annual Summer Institute in Survey Research Techniques from 
July 24 to August 18, 1950. The following courses will be offered: 
Introduction to Survey Research, Survey Research Methods, Sampling 
Methods in Survey Research (Introductory), Sampling Methods in Sur- 
vey Research (Advanced), Mathematics of Sampling, Statistical Methods 
in Survey Research, and Techniques of Scaling. In addition, the intro- 
ductory courses will be given from June 26 to July 21. This will permit 
students who are attending the full eight-week summer session of the 
University (June 26 to August 16) to register for the introductory 
courses during the first four weeks. It is expected that this special 
session will attract men and women employed in business or govern- 
mental research or other statistical work and university instructors and 
graduate students with a particular interest in this area of social science 
research. All courses are offered for graduate credit and students must 
be admitted by the Graduate School... Yale University’s Committee 
on Statistics offered a special course in applied statistics during recent 
months. The following speakers with their topics took part in this 
course. J. W. Tukey (Princeton) Order statistics, signed ranks, W. J. 
Youden (Bureau of Standards) Measurements in the physical sciences, 
D. F. Votaw, Jr. (Yale) Matrix theory for statisticians, P. J. Rulon 
(Harvard) Discriminant function in psychology and education, H. 
Hotelling (North Carolina) Multivariate analysis, C. I. Bliss (Yale) 
Discriminant function in biological assay, F. Mosteller (Harvard) Sta- 
tistical problems in social psychology, T. W. Anderson (Columbia) Use 
of stochastic equation models in economics, R. A. Fisher(Cambridge) 
Enumeration and recombination in genetics, W. G. Cochran (Johns 
Hopkins) Estimation of components of variance and their use, W. Allen 
Wallis (Chicago) Prediction of future observations, sequential analysis, 
statistics of the Kinsey Report, S. S. Wilks (Princeton) Industrial 
statistics, H. F. Dorn (National Institute of Health) Statistical problems 
in the measurement of morbidity, and Donald Mainland (Dalhousie) 
Statistics in clinical medicine. During his stay at Yale, each visiting 
lecturer was available for a limited number of individual or group con- 
ferences. ... L. Otis Emik is in the Commissioned Corps of the United 
States Public Health Service. He is Scientist with the Statistical Branch, 
Communicable Disease Center at Atlanta. His duties are to develop a 
program of Veterinary Statistics. Both census and research programs 
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are being organized, dealing with those diseases which affect both man 
and animals. Mr. Emik before 1950 was with the United States Sheep 
Experiment Station and Western Sheep Breeding Laboratory at Dubois, 
Idaho. ... E. L. LeClerg spent the month of January at Mayaguez, 
Puerto Rico, giving daily lectures on experimental design to the technical 
staff of the Federal Experiment Station. He is Research Coordinator 
in Field Crops for the Agricultural Research Administration of the 
United States Department of Agriculture at Washington, D.C... . 
F. E. Satterthwaite is quality control engineer of the Plastics Division 
in the Chemical Department of the General Electric Company in Pitts- 
field, Massachusetts, according to an announcement by F. W. Warner, 
engineering manager of the division. Mr. Satterthwaite has been quality 
control engineer for the product service division of the Company in 
Bridgeport, Connecticut. ... B. Schneider, formerly with the University 
of West Virginia, Morgantown, is now with the Department of Animal 
Husbandry, The State College of Washington, Pullman. 
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